The 511 keV emission from positron annihilation in the Galaxy 
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The first gamma-ray line originating from outside the solar system that was ever detected is the 
511 keV emission from positron annihilation in the Galaxy. Despite 30 years of intense theoretical 
and observational investigation, the main sources of positrons have not been identified up to now. 
Observations in the 1990's with OSSE/CGRO showed that the emission is strongly concentrated 
towards the Galactic bulge. In the 2000's, the SPI instrument aboard ESA's INTEGRAL 7-ray 
observatory allowed scientists to measure that emission across the entire Galaxy, revealing that 
the bulge/disk luminosity ratio is larger than observed in any other wavelength. This mapping 
prompted a number of novel explanations, including rather "exotic" ones (e.g. dark matter an- 
nihilation). However, conventional astrophysical sources, like type la supernovae, microquasars 
or X-ray binaries, are still plausible candidates for a large fraction of the observed total 511 
keV emission of the bulge. A closer study of the subject reveals new layers of complexity, since 
positrons may propagate far away from their production sites, making it difficult to infer the 
underlying source distribution from the observed map of 511 keV emission. However, contrary to 
the rather well understood propagation of high energy (>GeV) particles of Galactic cosmic rays, 
understanding the propagation of low energy (~MeV) positrons in the turbulent, magnetized 
interstellar medium, still remains a formidable challenge. We review the spectral and imaging 
properties of the observed 511 keV emission and we critically discuss candidate positron sources 
and models of positron propagation in the Galaxy. 
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The existence of a particle with equal mass but oppo- 
site c harge to that of the electron was predicted by iDirad 
(ll93lh . who named it " anti-electron" . Unaware of Dirac's 
prediction, lAndersonl (|l932[ ) found the first experimental 
hints for such a particle in cloud-chamber photographs of 
cosmic rays, and he named it p ositron. His finding was 
confir med the following year by iBlackett and Occhialinil 
(1933), w ho identified it with Dirac's anti- electron. One 
year later. iKlemperer and Chadwickl ([l934) detected the 
characteristic 7-ray line at 511 keV resulting from e~-e + 
annihilation, a convincing proof that positrons are in- 
deed elec tron's antipa r ticles. That same year, Croatian 
physicist iMohorovicid (|1934f ) predicted the existence of 
a bound system composed of an electron and a positron 
(analogous to the hydrogen atom, but with the proton re- 
placed by a positron), which he calle d "electr u m". T his 
state was experimentally found by iDeutschl (|l95ll ) at 
MIT and became known as positronium. 



For about 30 years after their discovery, all detected 
po sitrons we r e o f terres trial origin. Those d e tected 
by lAndersonl (|l932h and IBlackett and Occhialinil (|l933l ) 
were created by cosmic ray interactions with m olecules 
in Earth's atmosphere. Ijoliot and Curid (|l934l ) identi- 
fied another positron producing process: /3 + radioac- 
tivity of artificially created unstable nuclei. The first 
positrons of extrater restrial origin were reported by 
Ide Shong et all ([l964) , who loaded a spark chamber on a 
stratosphe ric balloon to de tect positrons within the cos- 
mic rays. iGinzburj ( 1956h had already suggested that 
high energy p-p interactions in cosmic rays would pro- 
duce pions 7r + , which would decay to positrons (via \i 
decays) . The production r ate of those pions was evalu- 
ated by iPollack and Faziol (| 19631 ) who predicted a 7-ray 
flux from the Galaxy at 511 keV of ~10 -3 cm -2 s _1 . 



The properties of e~-e + annihilation were explored in 
the 1940's. Direct e~-e + annihilation produces a single 
7-ray line at 511 keV, while the annihilation of positro- 
nium produces a composite spect rum with a lower en- 
ergy continuum and a 511 keV line ( Ore and PowelllJl949l 
and Sec. III. A)) . ISteckerl ([l969) was the first to point out 
that in the conditions of the interstellar medium, most 
positrons would annihilate after positronium formation; 
this would reduce the 511 ke V flux from cosmic rays t o 
values lower than evaluated bv lPollack and Faziol ([1963). 

The 511 keV emission of e + annihilation was first de- 
tected from the general direction of the Galactic center 
in the early 197Q's . by balloon-borne in struments of low 
energy resolution ([Johnson et a/I Il972h . It was unam- 
biguously identifie d a few years la t er wi th high resolu- 
tion Ge detectors ( Leventhal et a/.Lll978l ). It is the first 
and most intense 7-ray line originating from outside the 
solar system that was ever detected. Its flux on Earth 
(~10 -3 cm -2 s _1 ), combined with the distance to the 
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Galactic center (~8 kpc 1 ), implies the annihilation of 
-2 10 43 e+ s" 1 fSec. IILB.3p . releasing a power of -10 37 
erg s _1 or ~10 4 L in 7-rays. Assuming a steady state, 
i.e. equality between production and annihilation rates of 
positrons, one should then look for a source (or sources) 
able to provide ~2 10 43 e + s _1 . If the activity of that site 
were maintained to the same level during the ~10 10 yr of 
Galaxy's lifetime, a total amount of positrons equivalent 
to ~3 M would have been annihilated. 

Several years earlier, the Sun had already become the 
first astrop hysical laboratory for the study of positron an- 
nihilation (Crannell et al ., 1976). The solar annihilation 
7-ray line had been detected with a simple Nal instru- 
ment aboard the OS 0-7 satellite ([Chupp et all Il975l ). 
The Solar Maximum Mission (SMM), designed for solar 
flare observations and launched in 1980, featured a 7-ray 
spectrometer with exceptional stability. Based on de- 
tailed measurements with SMM, positrons in solar flares 
were found to originate from flare-accelerated particles 
when hitting the upper photosphere. Nuclear interac- 
tions of flare-accelerated protons and ions with atomic 
nuclei of the photosphere produce radioactive nuclei and 
pions that d ecay by emission of positrons, which annih i- 
late locally ([Murphy et all l2QQ5bl : iRamatv et aE 1 19831 ). 

Imaging the Galaxy in annihilation 7-rays was con- 
sidered to be the exclusive way to identify the cosmic 
e + sources (assuming that the spatial morphology of 
the 7-ray emission reflects the spatial distribution of 
the sources, i.e. that positrons annihilate close to their 
production sites). Because of the difficulties of imag- 
ing in the MeV region, progress was extremely slow in 
that field: only in the 1990s were the first constraints 
on the spatial distribution of the 511 keV emission in 
the inner Galaxy obtained by the OSSE instrument 
aboard the Compt o n Gamma Ray Obse rvatory (CGRO, 
ICheng et all Il997l iPurcell et a/1 Il997l ). The most reli- 
able imaging of the 511 keV emission was obtained by 
the SPI instrument aboard ESA's INTEGRAL Gamma 
Ray Observatory: the emission is strongly concentrated 
in the inner Galaxy (the bulge, iKnodlseder et all l2005f ) 
and a weaker emission i s detect ed from the Galactic disk 
( Weidenspointner et all l2008al ), unlike the situation in 
any other wavelength. 

Several candidate sources of positrons were proposed 
over the years: radioactivity from j3 + decay of unsta- 
ble nuclei produced in stellar explosions, high energy in- 
teractions occurring in cosmic rays or near compact ob- 
jects (like pulsars and X-ray binaries) or the supermassive 
black hole in the Galactic center, etc. For a long time, 
radioactivity from 56 Co produced in thermonuclear su- 
pernovae (SNIa) appeared as the most promising candi- 
date, provided that just a few per cent of the released 
positrons could escape the supernova remnant and anni- 
hilate in the interstellar medium. However, none of the 
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candidate sources has a spatial pattern resembling the 
one of the detected 7-ray emission. In particular, the 
release of the first year of SPI data, revealing the bulge 
but not yet the disk, prompted a series of "exotic" expla- 
nations involving dark matter particles, superconducting 
cosmic strings, etc. The confirmation of disk emission 
a few years later, made such explanations lose much of 
their interest, without completely eliminating them up 
to now. 

The spectral analysis of the 511 keV emission es- 
tablished already at the late 1970's that most of 
the positrons annih ilate after positronium formation 
(Bussard et 

aHE979). 

This result constitutes an impor- 
tant diagnostic tool for the physical properties of the 
annihilation medium, as analyzed in Guessoum et al. 
(1991). Only recently, i.e. in the 2000's, was it realized 
that the spectral analysis may also provide important 
hints on the e + source(s). In particular, positrons appear 
to annihilate at low energies, while in most candidate 
sources they are produced at relativistic energies; dur- 
ing the long period of their slowing down, positrons may 
travel far away from their sources, making the detected 7- 
ray emission useless as a tracer of their production sites. 
Unfortunately, propagation of low energy positrons in the 
turbulent, magnetized interstellar plasma of the Galaxy, 
is poorly understood at present. 

In this paper we present a synthetic view of the various 
facets of this complex issue, concerning the production, 
propagation and annihilation of positrons in the Galaxy, 
in relation with the characteristic signature of that an- 
nihilation, namely the 511 keV emission. The paper is 
structured as follows: 

In Sec. [II] we first present a historical account of ob- 
servations of the 511 keV emission, which illustrates the 
difficulties of 7-ray line astronomy in the MeV range. We 
also provide a summary of the latest SPI/INTEGRAL 
data analysis, of relevant observations in other wave- 
lengths and of the constraints imposed on the e + sources. 

Sec. Ill provides astrophysical background material 
concerning the stellar and supernova populations of the 
Milky Way, as well as the properties of the various phases 
of the interstellar medium (ISM) and of the Galactic 
magnetic fields in which positrons propagate. This mate- 
rial can be skipped at a first reading by astronomers, but 
it contains important and updated information, which is 
used in all other sections. 

In Sec. IV we discuss the physical processes and candi- 
date sources of positron production (radioactivity from 
stars and supernovae, high energy processes in cosmic 
rays, compact objects and the central supermassive black 
hole, dark matter, other "exotica", etc.) and, in some 
cases, we present new estimates of their e + yields. We 
discuss critically the properties of those sources, in the 
light of the observational constraints presented in Sec. II. 

Sec. V summarizes the various physical processes of e + 
slowing down and annihilation, taking into account the 
properties of the ISM, as well as the corresponding 7- 
ray spectral signature. The spectral analysis of the SPI 
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data is then used to constrain the energy of the emit- 
ted positrons (thus eliminating some of the candidate 
sources) and to constrain the properties of the annihila- 
tion medium. The results of this analysis offer some hints 
for e + propagation away from the sources. 

The intricacies of low energy positron propagation are 
discussed in some depth in Sec. VI, in the light of recent 
work. The implications of e + propagation for the 511 keV 
emission are also discussed. This is not a mature topic 
yet, and our poor knowledge of the plasma properties 
in the inner Galaxy prevents any definitive conclusions. 
A synthetic summary of the subject and directions for 
future research are presented in the last section. 



Spin and momentum conservation constrain the release 
of annihilation energy in the form of photons. Para- 
positronium annihilation releases two photons of 511 keV 
each in opposite directions (as in the case of direct e~-e + 
annihilation). Ortho-positronium annihilation requires a 
final state with more than two photons from spin con- 
servation; momentum conservation distributes the total 
energy of 1022 keV among three photons 3 producing a 
continuum of energies up to 511 keV (Fig. [I]). The cor- 
responding lifetimes before annihilation (in vacuum) are 
1.2 10 -10 s for p-Ps and 1.4 10 -7 s for o-Ps. 

If a fraction fp s of the positrons annihilate via positro- 
nium formation, then the 3-photon 7-ray continuum of 
ortho-positronium will have an integrated intensity of 



II. OBSERVATIONS 

Gamma-ray observations in the MeV domain (from 
~100 keV up to a few MeV) provide access to Galactic 
positrons through three main windows: 

i) Emission from e + annihilation at sub-relativistic en- 
ergies, with its prominent 511 keV line and the associated 
three-photon continuum from positronium annihilation; 

ii) Continuum emission at energies E>0.5 MeV 
from energetic positrons propagating through interstel- 
lar space and annihilating " in-flight" ; 

hi) Emission of characteristic gamma-ray lines from 
radioactive nuclei, such as 26 Al and 44 Ti, which also pro- 
duce positrons by /3 + -decay. 

In this section, we describe the relevant observations 
and the constraints that they impose on our understand- 
ing o f Galactic positrons (see also iDiehl and Leisingl 
2009). We start with a brief description of the radiative 
signatures of positron annihilation. 



A. Radiative signatures of positron annihilation 

The annihilation of a positron with an electron releases 
a total (rest-mass) energy of 1022 keV in the form of two 
or more photons. Direct annihilation of a e~-e + pair at 
rest produces two photons of 511 keV each. The situa- 
tion is more complex in the case of positronium (Ps). The 
ground state of positronium has two total spin states, de- 
pending on the relative orientations of the spins of the 
electron and the positron. The singlet state has antipar- 
allel spins, total spin S=0, is denoted as 1 S'o and is known 
as para-positronium (p-Ps). The triplet state has parallel 
spins, total spin S=l, is denoted as 3 Si and is known as 
ortho-positronium (o-Ps). From the (2S+1) spin degen- 
eracy, it follows that Ps will be formed 1/4 of the time 
in the p-Ps state and 3/4 of the time in the o-Ps state. 2 



2 The energy difference between the two spin states ("hyperfine 
splitting") is 8.4 10 -4 eV. Transitions between these states sim- 
ilar to the spin-flip transition in hydrogen, which produces the 
astrophysically-important 21 cm line, are unimportant due to the 
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(i) 



The remaining fraction 1 — fp s will annihilate directly 
to 2 photons of 511 keV each, to which should be added 
the 2-photon contribution of the para-Positronium state; 
thus, the 2-photon (511 eV line) intensity will be: 



727 oc 2(1 -f Pa ) 



1.5 h 



(2) 



By measuring the intensities of the 511 keV line and of 
the Ps continuum one can then derive the positronium 
fraction 



fp s = 



8 /37A27 
9 + 6/ 37 // : 



27 



(3) 



which offers a valuable diagnostic of the physical condi- 
tions of the ISM where positrons annihilate, (see discus- 
sion in Sec. V.E). 

The state of Ps can be formally treated as the one of 
hydrogen atom. The corresponding Schroedinger equa- 
tions are identical, with the reduced mass being half of 
the electron mass in the case of Ps. Because of that, 
the frequencies of the de-excitation spectral lines of ar e 
roughly half of those of the H atom ( Canter et aZ.L fl975| ) . 
Radiative recombination of ortho-Ps could, under certain 
circumstances, be observed i n the near infra-red by next 
gener ation instruments (see lEllis and Bland- Hawthorn! 
l2QQ9h , 



short Ps lifetimes (see text). 
3 Annihilation into a larger number of photons (an even num- 
ber for para-positronium, an odd number for ortho-positronium) 
is possible, but the corresponding branching ratios are negligi- 
ble (~10 -6 for four photons in the case of para-Positronium or 
five photons in the case of ortho-positronium). Even lower are 
the branching ratios for annihilation into neutrino-antineutrino 
pairs. Finally, a single-photon annihilation is also possible, pro- 
vided momentum conservation is obtained through positronium 
being bound to another particle, such as a dust grain. None of 
those cases is important in astrophysical settings. 
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FIG. 1 Spectrum of ortho-positr onium annihilation w ith the 
three-photon continuum (from lOre and Powell 119491 ) . 



B. Observations of ~ MeV emission from e + annihilation 

1. Early balloon and satellite observations 

The first evidence for Galactic 511 keV emission was 
obtained in the early 1970's, through a s eries of balloon 
flights by teams from Rice University. IJohnson et all 
announced the first detection of a celestial gamma- 
ray line originating from outside the solar system. Using 
a sodium iodine (Nal) scintillation detector they observed 
a spectral excess with a flux of 1.8 x 1(T 3 ph cm 2 s 1 
at an energy of 473 ± 30 keV during a balloon flight in 
1970. A second balloo n flight of the same team in 1971 
confirmed this signal ( Johnson and Havmesl 0~973). Al- 
though the authors mentioned e + annihilation as a pos- 
sible origin of the observed feature, the significant offset 
between observed (473 keV) and expected (511 keV) line 
centroid led them to conclude that the feature was, per- 
haps, due to radioactive decay of unknown origin. 

ILeventhal (|l973h proposed an interesting alternative 
explanation that could reconcile observations with a 
positron origin of the feature. As e + annihilation in 
the interstellar medium occurs mainly via positronium, 
a superposition of a narrow line at 511 keV and a con- 
tinuum emission below 511 keV is expected (Fig. [1]). 
Because of the poor spectral resolution of the Nal detec- 
tor and the low signal/noise ratio of these balloon flights 
these two contributions could not be disentangled, and 
were reported as a single emission peak at an energy of 
~ 490 keV, in reasonable agreement with the observa- 
tions. 

It took five y e ars before a group from Bell-Sandia 
( Leventhal et all Il978l ) could confirm this conjecture. 
With their high-resolution balloon-borne germanium 
(Ge) detector they could separat e the line and c ontin- 
uum components of the emission. ILeventhal et all (|l978f ) 
located the narrow component (FWHM=3.2 keV) at an 



energy of510.7±0.5 keV, consistent with the expecta- 
tions for e + annihilation at rest. The observed line flux 
of (1.2 ± .2) 10~ 3 ph cm~ 2 s~ 1 was below the value re- 
ported by I Johnson et all (|l972[ ), as expected if the ear- 
lier m easurements wer e a su perposition of two compo- 
nents. ILeventhal et all ( 1978h also detected the positro- 
nium continuum component, and the comparison of line 
and continuum intensities implied that 92% of the anni- 
hilations occurred after the formation of a positronium 
atom. 

Subsequent observations of the Galactic center 
by different balloon-borne instruments in the 1970's 
found a surprising v ar iabilit y of the 511 k eV line 
flux (lAlberrihe et all Il98ll: [Gardner et all Il982l : 
Havmes et all 1 19751; ILeventhal et all 198Qh . But 



Albernhe et all ((1981) recognized that the flux measured 



by the various balloon experiments increased with 
increasing size of the detector's field of view, which 
could mean that the anni hilation emission wa s extended 
along the Galactic plane. Rie gler et all ( 198lh proposed 
a different scenario, based on analysis of their HEAO-3 
satellite data. These data showed a decline by almost 
a factor three of the 511 keV flux between fall 1979 
and spring 1980, suggesting that positron annihilation 
was variable in time. From the At ~6 months interval 
between the observations they inferred a maximum size 
of Ar rsj cAt ~ 0.3 pc of the annihilation site, which 
implies gas densities for the annihilation medium of 10 4 
- 10 6 cm -3 . These extreme conditions suggest that the 
positrons were produced by a compact source such as 
a massive black hole wi thin 4° of the Galactic center 
([Lingenfelter et aZ.L fl98lh . 

While balloon-borne experiments see med to establish 
the variability of the 511 k e V em ission ([Leventhal et all 
1 19821 1 19861 ; iPaciesas et all Il982[ ) , contemporaneous ob- 
servations by the Solar Maximum Mission (SMM) satel- 
lite did not confirm such a trend. SMM carried a 
Nal-detector with a large field of view (130°) and pro- 
vided the first long-term monitoring of the inner Galaxy 
(1980-1987). The variability of the 511 keV emission 
was c onstrained to be less than 30% ([Share et all 1 19881 , 
1990). The apparent disagreement between the balloon 
and SMM observations could still be understood by as- 
suming an extended distribution of 511 keV emission 
along the Galactic plane; however, in order to recon- 
cile the observations with a time variable source one 
had to adopt rather complex scenarios. For example, 
iLingenfelter and Ramatvl ( 19891 ) suggested the combina- 
tion of a steady, extended 511 keV emission along the 
Galactic plane and a compact variable source at the 
Galactic center (assumed to be active from 1974 through 
1979) in order to explain all data available at that time. 
Such a scenario could not be ruled out, since no imaging 
of the 511 keV emission had been achieved and the mor- 
phology and spatial extent of the positron annihilation 
emission were essentially unconstrained. 

The hypothesis of a time variable central Galactic 
positron source was revitalized in the early 1990s by 
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the observation of transient 7-ray line features with 
the SIGMA telescope. The french coded mask im- 
ager SIGMA was launched in 1989 on board the soviet 
GRAN AT satellite. It was the first imaging 7-ray in- 
strument, with an angular resolution of 15 arcmin and 
it used Nal detectors covering the energy range 35-1300 
keV. In October 1991, an unusual spectrum was observed 
from IE 1740.7-2942, during an outbreak of this hard 
X-ray source which lasted ~17 hours ( Bouchet et all 
Il991 : ISunvaev et all Il99lf ). Superimposed on a candi- 
date typical black hole continuum spectrum, there ap- 
peared a strong (flux F ~ 10 -2 ph cm -2 s _1 ) and 
broad (FWHM ~ 200 keV) emission line centered at 
about 440 keV. If interpreted as broadened and red- 
shifted annihilation line, this observation seemed to make 
IE 1740.7-2942 be the long sought compact and vari- 
able source of positrons. Follow-up observations led to 
the classification of I E 174 0.7-2942 as the first micro- 



quasar ([Mirabel et all |1992[ ) : a binary system involving 
a compact object (neutron star or black hole, see sec. 
IV. B. 3) accreting material from its companion and emit- 
ting part of the accreted energy in the form of jets. It 
was therefore proposed that IE 1740.7-2942 would occa- 
sionally emit jets of positrons (produced in e~-e + pairs), 
some of which would annihilate in the inner edge of 
the accretion disk as presumably observed by SIGMA; 
the remaining positrons would eventually lose their en- 
ergy and give rise to time- variable narrow 511 keV line 
emission. Different SIGMA teams also reported narrow 
and/or broad 7-ray lines near 511 keV, lasting for a day 
or so, from the t r ansient X-ray source "Nova Muscae" 
([Goldwurm et all Il992fc ISunvaev et all Il992h and the 
Crab nebula (iGilfanov et all Il994). Another transient 
7-ray li ne source was disco vered from archival HEAO 1 
data bv lBriggs eTall (|l995h . 

However, the line feature seen by SIGMA was not 
seen in simultaneous obser vations of IE 174 0.7-2942 per- 
formed with t he OSS E ([Jung et all Il995h and BATSE 
(Sm ith et all [l996a) instruments aboard the NASA 
Compton Gamma Ray Observatory (CGRO, launched 
in 1991). Besides, BATSE data did not confirm the 
transient event seen by SIGMA from the Crab nebula 
(Smi th et all Il996al ). Moreover, the search on 6 years 
of BATSE data did not reveal any transient line fea- 
ture from a n y direc tion of the sky ([Cheng et all 1 19981 : 
ISmith et all Il996bh . Similarly, 9 years of SMM data 
did not show any transient event f rom the Galactic cen- 
ter direction ( Harris et all [l994a) or the Crab nebula 
([Harris et a/lll994bl ). A reanalvsis of HEAO 3 data then 
revealed that the d rop in 511 keV flux re ported earlier by 
Ried er et all (|l98lh was not significant ( Mahonev et all 
1994). Thus, the idea of a steady 511 keV Galactic emis- 
sion was gradually established. 

The contradictory results obtained during the 1980s 
and early 1990s provide a dramatic illustration of the 
difficulties affecting the analysis of 7-ray line data. In 
this domain, astrophysical signals rarely exceed the in- 
strumental background by more than a few percent and 



any systematic uncertainty in the treatment of the back- 
ground immediately disturbs the analysis. In particular, 
the time variability of the instrumental background (due 
to changing radiation environments along the orbital tra- 
jectory, or due to solar activity) can easily fake time vari- 
able signals. In addition, hard X-ray sources often exhibit 
highly variable continuum emission components that may 
further affect the data analysis and require their proper 
modeling; this concerns, in particular, the densely pop- 
ulated regions towards the Galactic center, which were 
not spatially resolved by older instruments. 



2. Early mapping of the spatial distribution of e + annihilation 

Before the launch of CGRO in 1991 with its OSSE 
collimated (11.4° x 3.8°) spectrometer, the spatial dis- 
tribution of 511 keV line emission was only poorly 
constrained. Hypotheses on a possible extent of the 
emission were mainly b ased on theoretical expectations 
( Kozlovskv et a/.Ul987l ), on the differen t fluxes received 
by detectors with different fields of view ([Albernhe et all 
119811 : iDunphv et all . [l983), and on a marginal detection 
of the 511 keV line near ~ 25° Ga lactic longitude with 
the balloon-borne GR IS telescope ([Gehrels et all Il99ll : 
iLeventhal eTall . [l989) . Nine years of OSSE observations 
drastically improved this situation. 

OSSE data could clearly exclude a single point 
source as the origin of observed 511 keV line emission 
(jPurcell et all Il994f ). The data were best understood 
in terms of an extended source consisting of a symmet- 
rical bulge (centered on the G alactic center ) and emis- 
sion from the Gala ctic plane. ICheng et all ([l997) and 
iPurcell et all (|l997h established the first 511 keV line 
emission map of the central Galactic ridge (Fig. [2]). Be- 
yond the aforementioned components, there was hint of a 
third component located at Galactic coordinates of lon- 
gitude / ~ —2° and latitude b ~ 12°, dubbed the Posi- 
tive Latitude Enhancement (PLE). However, the inten- 
sity and morphology of t his feature were on ly weakly 
constrained by the data (Mil ne et all l2001bl ), and the 



non- u niform exposure of th e sky m a y have biased the sk y 
maps (IVon Ballmoos efall. l2003l)jKinzer et all (|l996h . 



IKinzer et all ([200 lh and lMilnd (|l998h studied the spatial 
distribution of the continuum emission from positronium 
annihilation and concluded that it follows closely the dis- 
tribution of the 511 keV line. However, no PLE was 
visible in the continuum emission image ([Milne et all 
l2001al) . 

At this point one should mention that images in the 
hard X-ray and soft 7-ray domains are obtained through 
complex non-linear iterative deconvolution techniques, 
and they generally represent only a family of solutions, 
which explains the observed data within the given sta- 
tistical and convergence constraints. The reader should 
be aware of this particularly important point when in- 
specting all images in this paper. For instance, other 
OSSE images of the Galactic 511 keV line emission are 
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FIG. 2 OSSE 511 keV line map of the galactic center re- 
gion (top) and cor responding exposure map (bottom) (from 
IPurcell et a/lll997h . 




FIG. 3 511 keV line map (top) and positronium continuum 
map ( bottom) derived from one year o f INTEGRAL/SPI data 
(from IKnodlseder et all 120051 and IWeidenspointner et all 
l2006l respectively). 



presented in ([Milne et a/ll2001al ) or ([Milne et aZ.I . l2QQ2h . 

Several models have been proposed to describe the 
spatial dis tribution of the ann i hilation emission observed 
by OSSE ([Kinzer et all l200ll . llg^Milng. et all l2001bt 
iPrantzosl . Il993fc IPurcell et all Il997l Il994t ?). They 



all had in common a two-component emission from a 
spheroid located in the inner Galaxy and from the ex- 
tended Galactic disk (see Sec. III. A for a detailed dis- 
cussion of the Galaxy's morphology). However, both 
morphology and relative intensity of these two compo- 
nents were only poorly constrained by the data; depend- 
ing on the adopted model, the spheroidal component was 
claimed to be dominant or sub-dominant, i.e the Galactic 
spheroidal-to-disk flux ratio was only constrained in the 
broad interval F sp h er /Fdisk ~ 0.2 — 3.3. The uncertainty 
on the total Galactic 511 keV line flux was also rather 
large, spanning the range 1-3 10 -3 ph cm _2 s _1 . 

Despite the considerable progress achieved by OSSE 
observations, the origin of Galactic positrons remained 
unclear. The data did not constrain the morphology 
of the 511 keV emission enough to clarify the under- 
lying source population(s). Yet, the strong concentra- 
tion of the 511 keV emission towards the Galactic bulge 
led several authors to suggest that the j3 + decay of ra- 
dioactive 56 Co, produced by Galactic Type la supernovae 
(SNIa) , should be the dominant Galactic pos itron source 
(Kinzer et al, 2001, 1996; Milne et al, 2002). The emis- 



sion from the Galactic disk was generally attributed 
to radioactive j3 + decays of 26 Al, 56 Co and 44 Ti pro- 



duced by a variety of stel lar sources ([Kinzer et ^D,[l996; 
IPurcell et aH . ll997Lll99l . In fact, 26 Al had already been 
detected from the inner region of the Galaxy through its 
characteristic 7 -ray line at 1809 keV in the 1979/1980 
HEAO-C data ([Mahonev et a/lll982f ). and its contribu- 
tion to Galactic-disk e + production was established (Sec. 
II.B.2 and IV.A.2). 



3. Imaging with INTEGRAL/SPI 

With the la u nch o f ESA's INTEGRAL observatory 
( Winkler et al 1 l20Q3h in 2002 for a multi-year mission, 
a new opportunity became available for the study of 
Galactic e + annihilati on. The SPI imaging spectrometer 
( Vedrenne et all [20031 ) combines for the first time imag- 
ing with high-resolution spectroscopy. The spatial resolu- 
tion of 3° (FWHM) of SPI, though inferior to telescopes 
optimized for slightly lower energies (SIGMA, IBIS), is 
superior to that of OSSE; its energy coverage and sen- 
sitivity around the annihilation line and its large field 
of view allow an improved study of the 511 keV emis- 
sion morphology. The spectral resolution of ~ 2.1 keV 
(FWHM, at 0.5 MeV) is comparable to that of other Ge 
detectors employed on balloons or the HEAO 3 satellite, 
allowing for a spatially resolved fine spectroscopy of the 
signal (including the underlying continuum emission). 

The first 511 keV line and positr onium continuum 
all- sk y maps have been presented by IKnodlseder et all 
( 2005 ) and IWeidenspointner et all (|2006l ), respectively. 
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FIG. 4 511 keV line map derived from 5 years of INTE- 
GRAL/SPI data (from Wei denspointner et all l2008aT ). 



based on approximately one year of SPI data (Fig. [3]). 
The two maps are compatible with each other (within 
their uncertainties), suggesting that the positronium 
fraction does not vary over the sky. The images illustrate 
the remarkable predominance of the spheroidal compo- 
nent. In contrast to OSSE data, which suggested a rela- 
tively strong disk component, the Galactic disk seemed to 
be completely absent in the first year SPI images. Model 
fitting indicated only a marginal signal from the Galac- 
tic disk, corresponding t o a bulge-to-disk flux ratio > 1 
( Knodlseder et all l2QQ5h . This strong predominance of 
the Galactic bulge, unseen in any other wavelength, stim- 
ulated " unconvention al" models invol ving dark matter 
(Sec. IV.C). However. [Prantzosl (|2006l ) pointed out that 
the data could not exclude the presence of disk emission 
of a larger latitudinal extent (resulting from positrons 
propagating far away from their sources), which could be 
rather luminous and still undetectable by SPI, because 
of its low surface brightness. 

After accumulating 5 years of INTEGRAL/SPI data 
the 511 keV line emission all-sky image revealed also 
fainter emission extending along the Galactic plane 
(Fig. [4]). With a much improved exposure with respect 
to the first year (in particular along the Galactic plane), 
511 keV emission from the Gala c tic dis k is now clearly 
detected ( Wei denspointner et all 2008a). However, the 
detailed quantitative characterization of components of 
511 keV emission requires parameterizing these in the 
form of (necessarily idealized) spatial emission models 
fitted to the data. No unique description emerges at 
present, since both the spheroid and the disk may have 
faint extensions contributing substantially to their total 
7-ray emissivities. It turns out that the bulge emission 
is best described by combining a narrow and a broad 
Gaussian, with widths (FWHM, projected onto the sky) 
of 3° and 11°, respectively. Another, more extended com- 
ponent is needed to fit the data, a rather thick disk of 
vertical extent 7° (FWHM projected on the sky). The 
model implies a total e^ 



s _1 and a spheroid/disk ratio of 1.4 (Table UJ. It should 
be noted, however, that alternative models, involving ex- 
tended components of low surface brightness (thus far 
undetected by SPI) are also possible. One such alterna- 



TAB LE I Two model fits of the G alactic 511 keV emission 
(from Weiden spointner et al. I l2008bh : fluxes, photon emissiv- 
ities and e + annihilation rates (computed for a positronium 
fraction of / ps =0.967, see Sec. II. B. 4). Notice that "thin" 
and "thick" disks have not the same meaning as in Sec. III. 
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tive (Weidenspointner et all l2QQ8bh involves a centrally 
condensed but very extended halo and a thinner disk 
(projected vertical extent of 4°), with a spheroid/disk 
ratio of 6 (Table U). 

With more SPI data, it was possible to proceed to 
more detailed constraints on the morphology of the disk 
emission. The flux in the disk component remains con- 
centrated to longitudes |/| < 50°; no significant 511 keV 
line emission has been detected from beyond this interval 
so far. The accumulated SPI data yield a flux from nega- 
tive longitudes of the Galactic disk that is twice as large 
as the flux from an equivalent region at positive longi- 
tudes. The significance of this asymmetry is still rather 
low, about ~ Act. Indications for such an asymmetry 
were already noticed in the OSSE data (M. Leising, pri- 
vate communication). It should be noted, however, that 
a different analysis of the s ame SPI data finds no ev i- 
dence for a disk asymmetry ([Bouchet et all 120081 l2Q10h , 
although it cannot exclude it, either. Clearly, clarifying 
the asymmetric or symmetric nature of the disk profile 
should be a major aim of the 511 keV studies in the years 
to come 4 . 



4. Spectroscopy with INTEGRAL/SPI 

Before INTEGRAL, the spectral shape of the positron 
annihilation emission was only poorly constrained by ob- 
servations. All high-resolution observation s suggested a 
mode s t line broadening of FWHM^ 2 keV ([Harris et al. , 
19981 : iLeventhal eTaU . Il993l : iMahonev et all Il994 : 



Smith et~a/l Il993h . The excellent spectral resolution of 



4 INTEGRAL will continue operations until 2012, at least. 
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TABLE II Results of spectral analysis of Galactic 511 keV 
emission from the region within 8°from the Galactic center. 
I n , r n , h and F b are the flux and width (FWHM) of the 
narrow and broad lines, respectively, 1^ is the flux of the 
ortho-positronium continuum and A c is the amplitude of the 
Galactic continuum at 511 keV. The first set of error bars 
refers to la statistical errors a nd the second set to systematic 
errors ffrom ljean et aZ.L l20Q6h . 



Parameters 




Measured values 


In 


(10 -3 s _1 cm" 


- 2 ) 


0.72 ± 0.12 ± 0.02 


r„ 


(keV) 




1.32 ± 0.35 ± 0.05 


h 


(10~ 3 s" 1 cm" 


2 ) 


0.35 ± 0.11 ± 0.02 


T b 


(keV) 




5.36 ± 1.22 ± 0.06 


I37 


(10~ 3 s" 1 cm 


- 2 ) 


4.23 ± 0.32 ± 0.03 




(10 -6 s _1 cm 


" 2 keV" 1 ) 


7.17 ± 0.80 ± 0.06 



SPI allows now for the first time to study the spectrum 
of this emission in great detail and for different regions. 

Spectral results for the Galact i c sph eroid al emission 
were p resented bv lChurazov et all ( 20051 ) and I Jean et al\ 
([2006f h based on the first year of SPI data. The 

line displays no s pectral shift, i.e . it has an energy 

£=511±0.08 keV (jChurazov etall l2QQ5h and it is com- 



posed of two spectral components (assumed to be rep- 
resented by Gaussians): a narrow line with a width of 
FWHM=1.3 ±0.4 keV and a broad component with a 
width of FWHM=5.4±1.2 keV (Fig.©. The width of the 
broad line is in agreement with the broadening expected 
from positronium annihilation via charge exchange with 
hydrogen atoms (see section V.B.2). The narrow line 
component contains ~ 2/3 of the total annihilation line 
flux while the broad one makes up the remaining ~ 1/3 
of the flux. Table [TT] summarizes the results of the spec- 
tral analysis of the Galactic 511 keV emission after the 
first year of SPI data. 

SPI also clearly detected the ortho-positronium contin- 
uum with an intensity tha t corresponds to a positronium 
fraction of f Ps =97±2 % (IJean et all 120061 : see Eq. [3|). 
This value is in good agreement w ith earlier measure - 
ments obtained by O SSE (97±3 %, lKinzer_ela/J , Il996h 
and TGRS (94±4 %. lHarris et~a/l Il998l ). 

The shape of the annihilation line and the relative in- 
tensity of the ortho-positronium continuum are closely 
related to the physical conditions such as temperature, 
ionisation stage and chemical abundances of the interstel- 
lar medium in which positrons annihilate. These condi- 
tions, obtained from the analysis of the measured spec- 
trum, are presented and discussed in Sec. V.E. Important 
complementary information on the energies of the anni- 
hilating positrons is obtained from the analysis of the ob- 
served continuum emission at somewhat higher energies 
(above 511 keV and into the MeV region), as discussed 
in the next section and Sec. V.B. 
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FIG. 5 Fit of the spectrum of the annihilation emission mea- 
sured by SPI with narrow and broad Gaussian lines and an 
ortho-positronium continuum. The p ower- law account for the 
Galactic diffuse continuum emission ( Je an et al .1 .12006). 



C. Relevant observations at MeV energies 

1. The MeV continuum 

Positrons are typically emitted at relativist ic energies, 
in some cases even far above 1 MeV (Sec. IV). They be- 
have essentially like relativist ic electrons of cosmic rays, 
by producing bremsstrahlung and inverse- Compton emis- 
sion while slowing down to the thermal energies (eV) 
of the interstellar medium, where they eventually an- 
nihilate. But positrons may also annihilate "in flight" 
while still having relativist ic energies, giving rise to a 
unique 7-ray continuum signature at energies above 511 
keV (as the center-of-mass energy is transferred to anni- 
hilation photons; Sec. V). The shape and amplitude of 
this 7-ray emission depend on the injection spectrum of 
positrons and the corresponding total annihilation rate. 
For positrons injected at low energies (of the order of 
~MeV, such as those released by radioactivity), the am- 
plitude of the in-flight annihilation continuum above 1 
MeV is quite small, while for sources injecting positrons 
at much higher energy (such as cosmic-ray positrons from 
pion decay), the annihilation 7-ray spectrum would ex- 
tend up to GeV energies and include a considerable 7- 
ray flux. The high energy 7-ray continuum above 1 MeV 
therefore constrains the energy and the annihilation rate 
of relativistic positrons, when all other sources of such 
high energy emission are properly accounted for. 

Diffuse Galactic continuum emission has been well- 
measured at least in the inner part of the Galactic 
disk (longitudes -30° < I < 30°) in the hard-X-ray 
through 7-ray regi me by INTEGRAL, OSSE, COMP- 
TER, and EGRET dBouchet et all I2QQ8I : iKinzer et all . 
Il999l : IStrong et all . Il99j). The spectrum of the underly- 
ing continuum emission in the 511 keV regi on is best rep- 
resented as a power-law with index 1.55 ([Bouchet et all 
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FIG. 6 Spectrum of the inner Galaxy as measured by various 
instruments, compared to various theoretical estimates made 
under the assumption that positrons are injected at high en- 
ergy: the four pairs of curves result from positrons injected at 
100, 30, 10 and 3 MeV (from top to bottom) and correspond 
to positrons propagati ng in neu t ral (s olid) or 50% ionized 
(dotted) media ffrom lSizun et al. , 2006). This constrains the 
injected positron energy (or, equivalent ly, the mass of decay- 
ing/annihilating dark matter particles) to a few MeV. 



2008) and is mostly due to emission from cosmic-ray elec- 
trons and positrons (Sec. IIV.B.2|) . The corresponding 
emission process es are modelled in detail in e.g. the 
GALPROP code (IStrong et q/1 l2Q07h , which includes 3D 
distributions of interstellar gas and photon fields, as well 
as all relevant interaction cross sections and constraints 
from near-Earth observations of cosmic-ray fluxes and 
spectra. This model reproduces well the entire range of 
7-ray observations, from keV to GeV energies; however, 
there exist tantalizing hints of residual emission in the 
MeV region, when comparison is made to COMPTEL 
measurements 5 : the data points appear to lie on the 
high side of model predictions. In view of the possible 
systematic uncertainties of such measurements, but also 
of the parameters of GALPROP code and the possible 
contributions of unresolved X-ray binaries, some (though 
little) room is still left for a contribution of in-flight e + 
annihilation to the MeV continuum. 

The physics of the in-flight annihilation of positrons 
will be analyzed in Sec. V. Here we simply note 
that the corresponding constraints on the injection en- 



5 The CGRO/COMPTEL data points have an uncertainty of up 
to a factor 2 due to the difficulty of producing skymaps with 
a Compton telescope with high background. The most reliable 
COMPTEL values come from a maximum-entropy imaging anal- 
ysis and are model-independent, but the zero level is uncertain 
and contributes to the systematic error. The enormous gap in 
sensitivity between the current Fermi mission (>30 MeV) and 
the 1-30 MeV range (sensitivity factor ~100 !) highlights the 
urgent need for new experiments in the MeV range. 



ergy of positrons wee p ointed out many years ago by 
lAgaronvan and Atovanl ( 198lf ). They showed that the 
positrons which are responsible for the Galactic 511 keV 
line cannot be produced in a steady state by the decay of 
the 7r + created in proton-proton collisions or else the in- 
flight annihilation emission s hould have been detected. A 
simil ar argumen t was u sed by lBeacom and Yiiksell ( 20061 ) 
and lSizun et al\ ([2006) to constrain the mass of the dark 
matter particle which could be the source of positrons in 
the Galactic spheroid (see section IV. C). If such parti- 
cles produce positrons (in their decay or annihilation) at 
a rate which corresponds to the observed 511 keV emis- 
sion, then their mass should be less than a few MeV, oth- 
erwise the kinetic energy of the created positrons would 
have been sufficiently high to produce a measurable 7- 
ray continuum emission in the 1-30 MeV range (Fig. [6]). 
The same argument allows one to constrain the initial 
kinetic energy of positrons and thus to eliminate several 
classes of candidate sources, like e.g. pulsars, millisecond 
pulsars, magnetars, cosmic rays etc., as major positron 
producers (Sec. IV. D). 



2. Gamma-rays and positrons from radioactive 26 AI 

26 Al is a long-lived (half-life t 1 / 2 =7.4 10 5 yr) radioac- 
tive isotope. It decays by emitting a positron, while the 
de-excitation of daughter nucleus 26 Mg emits a charac- 
teristic 7-ray line at 1808.63 keV. Based o n predict i ons o f 
nucl eosynthesis calculations in the 19 70's jArnettl (|l977[ ) 
and iRamatv and Lingenfelterl (|l977h suggested that its 
7-ray emission should be detectable by forthcoming space 
instruments. The detection of the 1809 keV line from the 
i nner Galaxy with the HEAO-C germanium spe ctromet er 
dMahonev et a/I ll 982) came as a small surprise ([Clavtonl 
Il984f ) because of its unexpectedly high flux (~4 10 -4 
cm -2 s _1 ). Being the first radioactivity ever detected 
through its 7-ray line signature, it provided direct proof 
of ongoing nucleosynthesi s in our Galaxy (see review by 
iPrantzos and DiehlL fl996l ). 

Several balloon experiments, and in particular the 
Gamma- Ray Spectrometer (GRS) aboard the Solar Max- 
imum Mission (SMM) rap idly confirmed the HEAO-C 
finding ([Share et a/| [l985). Early experiments had large 
fields of view (130° for SMM, 42° for HEAO-C) with no 
or modest imaging capabilities. The first map of Galactic 
1.8 MeV emission was obt ained with the CO MPTEL in- 
strument aboard CGRO ( Diehl et a/.U 19951 ), which had 
a spatial resolution of 3.8° (FWHM) within a field of 
view of 30°. The sky map derived from the 9-year sur- 
vey of COMPTEL is shown in Fig. U\ Unlike the 511 
keV maps of Fig. [3] and 3] the 26 Al emission is concen- 
trated along the Galactic plane (brightest within the in- 
ner Galactic radian) and is irregular, with emission max- 
ima aligned with spiral- arm tangents. The Cygnus region 
stands out as a significant and bright emission region. 
The "patchy" nature of the 26 Al7-ray emission suggests 
that massive stars are the most important contributors 



11 



% Cyqms 
AurigoM Per? 




FIG. 7 Map of Galactic 26 Al 7-ray emissio n after 9-year 
obser vations with COMPTEL/CGRO ffrom iPKischke et all 
l200lh . 



to Galactic 26 Al, as suggested at a time when the mor- 
phology of 26 Al emission was unknown ( Prantzosl Il99ll 
and Sec. IV. A. 2). It is consistent with the (statistically 
significant) similarity to the Galactic free-free emission 
map, which reflects electron radiation from HII regions 
ionized fr om the same mass ive stars that eventually re- 
lease 26 Al (|KnodlsederUl999f ). 

The total flux of 26 Al 7-rays depends slightly on the 
measuring instrument. In terms of statistical precision, 



the SMM result of 4.0±0.4 10" 4 ph cm^s^rad- 1 has 
been considered the canonical value. Imaging instru- 
ments, however, have consistently reported lower flux 
values of 2.6±0.8 1(T 4 ph cm-V^ad -1 (COMPTEL) 
and 3.1±0.4 10 -4 ph cm _2 s _1 rad _1 (SPI), respectively. 
The latest SPI value is compatible with the full range 
of measured values by other instruments (within statis- 
tical uncertainties), and we adopt it here. The detected 
flux translates into a decay rate of 26 Al which depends 
slightly on the adopted 3 D distribution of 26 Al in the 
Galaxy ([Diehl et all [2006). The most recent analysis of 
SPI data results in a r ate o f N 2 e= 4.3 10 42 s -1 or 2.7 
M /Myr ([Wang et a/ll2009l ). Assuming a steady state, 
i.e. equality between production and decay rates, this is 
also the present production rate of 26 Al in the Galaxy; 
recent models of massive star nuc leosynthesis can read- 
ily explain such a production rate ( Diehl et a/ll2006l and 
Sec. IV.A.2). 

Being predominantly a /3 + -emitter (with a branching 
ratio of / e + ,26=82%, see Table lVIIj) 26 Al is itself a source 
of positrons. The corresponding Galactic e + production 
rate is iV e + 5 26 = f e +, 26^26 ~ 3.5 10 42 s _1 . This consti- 
tutes a significant contribution to the total Galactic e + 
production rate (Sec. II. A. 3 and Table H]): 17% of the 
total e + annihilation rate and almost half of the (thick) 
disk in the double bulge+thick disk model, or 10% of 
the total and 70% of the thin disk in the Halo+thin disk 
model. We shall see in Sec. IV that positrons from other 
/3 + -decaying nuclei can readily explain the remaining disk 
emissivity, while the bulge emissivity remains hard to ex- 
plain. 



D. Summary of observational constraints 

The results of the analysis of Galactic 7-ray emissions 
in the MeV range can be summarized as follows: 

1) Intensity: The total rate of positron annihilation 
observed in 7-rays is at least L e +=2 10 43 s _1 , depending 
on the adopted source configuration. Most of it comes 
from the bulge (unless there is important emission from 
an extended, low surface brightness, disk). 

2) Morphology: The bulge/disk ratio of e + annihilation 
rates is B/D ~1.4; however, substantially different ratios 
cannot be excluded if there is important emission of low 
surface brightness (currently undetectable by SPI) either 
from the disk or the spheroid. About half of the disk 
emission can be explained by the observed radioactivity 
of 26 Al (provided its positrons annihilate in the disk). 
There are hints for an asymmetric disk emission with 
flux ratio F(l <0°)/F(l >0°)~1.8, which has yet to be 
confirmed. 

3) Spectroscopy: The ratio of the 511 keV line to the 
E<511 keV continuum suggests a positronium fraction 
of 97±2 % and constrain the physical conditions in the 
annihilation region. The observed continuum at ~MeV 
energies can be mostly explained with standard inverse 
Compton emission from cosmic ray electrons. A con- 
tribution from unresolved compact sources is possible, 
while a (small) contribution from high-energy (>MeV) 
positrons annihilating in flight cannot be excluded. 

These are the key observational constraints that should 
be satisfied by the source (s) and annihilation site(s) of 
Galactic positrons. We shall reassess them in the light of 
theoretical analysis in the end of Sec. IV and V. 



III. THE GALAXY 

The expected spatial distribution and intensity of the 
positron annihilation emission obviously depends on the 
corresponding distribution of the potential e + sources, as 
well as on the properties of the ISM in which positrons 
first slow down and then annihilate. One may distin- 
guish two types of e + sources, depending on whether 
their lifetimes (75) are shorter or longer than the lifetime 
of positrons in the ISM (r e +). Calculation of the total e + 
production rate requires in the former case (rs < r e + ) an 
estimate of (i) the Galactic birthrate Rs of the sources 
and (ii) the individual e + yields n e + (i.e. the average 
amount of positrons released by each source). In the lat- 
ter case (rs > T e +), the total number of such sources 
in the Galaxy Ns is required, as well as the individual 
e + production rate h e + of each source. In the former 
class belong supernovae or novae and the corresponding 
positron production rate is N e + = Rsn e +; in the lat- 
ter class belong e.g. low mass XRBs or millisecond pul- 
sars, and the corresponding positron production rate is 
N e+ =N s h e+ . 

The galactic distribution of any kind of stellar source of 
positrons is somewhat related to the distribution of stars 
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in the Milky Way. Similarly, the birthrate of any kind 
of positron source is somewhat related to the Galactic 
star formation rate. In this section we present a sum- 
mary of current knowledge about the stellar populations 
of the Milky Way and their spatial distribution and we 
discuss the birthrates of stars and supernovae. Since the 
slowing down and annihilation of positrons depend on 
the properties of the ISM, we present a brief overview 
of the ISM in the bulge and the disk of the Milky Way. 
Positron propagation depends also on the properties of 
the Galactic magnetic field, which are reviewed in Sec. 
III.D. Finally, the main properties of the Milky Way's 
dark matter halo are presented in Sec. III.E. 



A. Stellar populations 

The Milky Way is a typical spiral galaxy, with a total 
baryonic mass of ~5 10 10 M , of which more than 80% 
is in the form of stars. Stars are found in three main 
components: the central bulge, the disk and the halo, 
while the gas is found essentially in the plane of the disk. 
Because of its low m ass, estimated to 4 10 8 M i.e. less 
than 1% of the total, (jBell et a/lkoO/l ). the Galactic halo 
plays no significant role in the positron production. The 
bulge contains ~l/3 of the total mass and an old stellar 
population (age>10 Gyr). The dominant component of 
the Milky Way is the so-called thin disk, a rotationally 
supported structure composed of stars of all ages (0-10 
Gyr). A non negligible, contribution is brought by the 
thick disk, an old ( >10 Gyr) and kin ematically distinct 
entity identified by iGilmore and Reidl (|l983h . 

To a first approximation, and by analogy with external 
galaxies, the bulge of the Milky Way can be considered as 
spherical, with a density profile either exponentially de- 
creasing with radius r or of Einasto-type (p(r) ocexp(-A 
r a ). Measurements in the near infrared (NIR), concern- 
ing either integrated starlight observations or star counts 
revealed that the bulge is not spherical, but elongated. 
Recent models suggest a tri-axi al ellipsoid, but its exact 
shape is difficult to determine ( lLopez-Corredoira et all 
20051 iRattenburv et all l2QQ7h because of the presence of 
a Galactic bar. The mass of the b u lge lies in th e rang e 
1-2 10 10 M . (|Dwek et all Il995fc iRobin et all l2QQ3h . 
By comparing colour-magnitude diagrams of stars in the 
bulge and in metal-rich globular clusters, iZoccali et al\ 
( 20031 ) find that the populations of the two systems are 
co-eval, with an age of ^10 Gyr. 

The innermost regions of the bulge, within a few hun- 
dred pc, are dominated by a distinct, disk-like compo- 
nent, called the Nuclear Bulge which contains about 10% 
of the bulge stellar population (~1.5 10 9 M ) within 
a flattene d region of radius 230±20 pc and scaleheight 
45±5 pc ([Launhardt et all l2QQ2h . It is dominated by 
three massive stellar clusters (Nuclear Stellar Cluster or 
NSC in the innermost 5 pc, Arches and Quintuplex), 
which have a mass distribution substantially flatter than 



the classical Salpeter IMF 6 (|Fkerl . 120081 ) . Finally, in the 
center of the Milky Way, at the position of SgrA* source, 
lies the supermassive Galact ic black hole (SMBH) with 
a total mass of - 4 10 6 M (iGillessen et all 120081) . 

The Sun is located in the thin disk of the Milky 
way, at a distance of Rq ~8 kpc from the Galac- 
tic center; a recent evaluation, based on Cepheids, 
gives i^=7.94zb0.37(stat istical)±0.26(systematic) kpc 
( Groenewegen et all 120081 and references therein). Fur- 
thermore, the Sun is not located exactly on the plane, but 
at a distance from it z Q ~25 pc, as evaluated from the 
recen t analysis of the S loan Digital Sky Survey (SDSS) 
data (jJuric et aZ.L l2QQ8h . 

In studies of the Milky Way the solar neighborhood 
plays a pivotal role, since local properties can, in gen- 
eral, be measured with greater accuracy than global 
ones. The total baryonic surface density of the solar 
cylin der 7 is estimated to Ht=48.8 M pc~ 2 (|Flvnn et all 
2006), with ~13 M pc 2 belonging to the gas (see 
Sec. III.B). This falls on the lower end of the dynam- 
ical mass surface density estimates ( from kinematics of 
stars perpendicularly to the plane) which amo unt to 
S D =50-62 M ^ pc~ 2 (iHolmberg and Flvnnl . l200l ) or 57- 
66 M pc~ 2 (|Bienavme et all 120061 ). Thus, the values 
for the baryon content of the solar cylinder, summarized 
in Table IIII[ should b e considered rather as lower lim- 
its ( Flvnn et a/ll2006l ): the total stellar surface density 
could be as high as 40 M pc -2 . 

The density profiles of the stellar thin and thick disks 
can be satisfactorily fit with exponential functions, both 
in the radial direction and perpendicularly to the Galac- 
tic plane. The recent SDSS data analysis of star counts, 
with no a priori assumptions as to the functional form 
of the density profiles finds exponenti al disks with scale - 
lengths as displayed in Table HTJ (from ljuric et aD,|2008). 
The thin and thick disks cannot extend all the way to 
the Galactic center, since dynamical arguments constrain 
the spatial co-existence of such rotationally supported 
structures with the pressure-supported bulge. The ex- 
act shape of th e "central hole" of the disks is poorly 
known (see, e.g. iFreudenreichl 1 19981 : IRobin et~al\ , 120031 
for parametrizations), but for most practical purposes 
(i.e. estimate of the total disk mass) the hole can be 
considered as trully void of stars for disk radius R <2 
kpc. 

The data presented in this section (as summarized 
in Table IIII|h allow one to estimate the total mass of 



6 Stars are born with a mass distribution called Initial mass func- 
tion (IMF). Observed IMFs of young stellar clusters in the Milky 
Way and other galaxies have similar IMFs, with the upper part 
(M>1 M ) described by a power-law (dN/dM oc M~^ 1+x \ 
where X is th e slope of the IMF; in most cases, X=1.35 (as 
determined by Salpeter, 1955, for the local IMF. 

7 The solar cylinder is denned as a cylinder of radius 500 pc cen- 
tered on Sun's position and extending perpendicularly to the 
Galactic plane up to several kpc. 
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TABLE III Properties of the stellar populations of the thin 
and thick disk a . 



TABLE IV Physical properties (typical temperatures, hydro- 
gen densities and ionization fractions) of the different ISM 
phases in the Galactic disk. 









Thin 


Thick 


Mass density 


Po,o 


(M pc" 


" 3 ) 4.5 1(T 2 


5.3 10~ 


Surface density 


S 


(M pc" 


" 2 ) 28.5 


7 


Scaleheight 


H Q 


(pc) 


300 


900 


Scalelength 


L 


(pc) 


2600 


3600 


Star mass 


M D 


(10 10 M 


3 ) 2.3 


0.53 


(Age) 


<A> 


(Gyr) 


5 


10 


(Metallicity)© 


<[Fe/H]> 


(dex) 


-0.1 


-0.7 



a: The indice here denotes quantities measured at 
Galactocentric distance i?©=8 kpc. Average quantities are 
given within ( ). 



the thin and thick disks as M^ ? t^ n =2.3 10 10 M and 
MD y thick=0.5 10 10 M , respectively, in the galactocen- 
tric distance range 2-15 kpc. Overall, the disk of the 
MW is twice as massive as the bulge. 



B. Interstellar matter 

Interstellar matter is primarily composed of hydrogen, 
but it also contains helium (~ 9 % by number or 28 % by 
mass) and heavier elements, called "metals" (~ 0.12 % 
by number or 1.5 % by mass in the solar neighborhood). 
All the hydrogen, all the helium, and approximately half 
the metals exist in the form of gas; the other half of the 
metals is locked up in small solid grains of dust. In- 
terstellar dust manifests itself through its selective ab- 
sorption of starlight (leading to extinction, reddening, 
and polarization of starlight) and through its thermal 
infrared emission. Dust grains cover a whole range of 
radii, from a m i n <; 100 A to a max J> 0.25 /im, as implied 
from the overall shape of extinction curves which can 
be reproduced with a power-law dist r ibution in radius , 
N(a) da oc q~ 3 5 da (|Boulanger et 120001 : iKim et all 
1994: iMathis" et al. L 119771 1. Overa ll, gas and dust a ppear 
to be spatially well correlat ed ( Boulanger et all Il996t 
iBoulanger and PeraultL Il988l ). 

Interstellar gas can be found in molecular, atomic (cold 
or warm) and ionized (warm or hot) forms. The physical 
properties of the different gas component s in the Galac- 
tic disk were reviewed by IFerrierd (|200ll ) and are sum- 



marized in Table IIVI The gas properties in the Galactic 
bulge are less well established, but on the whole, all gas 
components appe ar to be hotte r and denser in the bulge 
than in the disk ([Ferriere et aZ.L l2QQ7h . 

Spatially, the molecular gas is confined to discrete 
clouds, which are roundish, gravitationally bound, and 
organized hierarchically from large complexes (size ~ 
20-80 pc, mass ~ 10 5 -2 10 6 M©) do wn to small clump s 
(size £ 0.5 pc, mass £ 10 3 M©) ([GoldsmithL Il987l ). 
The cold atomic gas is confined to more diffuse clouds, 
which often appear sheet-like or filamentary, cover a wide 



Phase 



T (K) n H (cm- 



Molecular (MM) 10 - 20 10 2 - 10 6 

Cold neutral (CNM) 20 - 100 20 - 100 

Warm neutral (WNM) 10 3 - 10 4 0.2 - 2 

Warm ionized (WIM) - 8000 0.1 - 0.3 
Hot ionized (HIM) - 10 6 0.003 - 0.01 



< 10 " 4 
4 10~ 4 - 10" 3 
0.007 - 0.05 
0.6-0.9 
1 



range of sizes (from a few pc up to ~ 2 kpc), and have 
random motions with typic al velocities of a few km s _1 
(Kulkarni and Heiles, 1987). The warm and hot compo- 
nents are more widespread and they form the intercloud 
medium. 



The different gas components also differ by their spa- 
tial distributions at larg e scales. The o bservational sit- 
uation was revie wed bv IFerrierd ( 200lh for the Galac- 
tic disk and by IFerriere et all ([20071 ) for the Galactic 
bulge. Fig. Ogives the radial variation of the azimuthally- 
averaged surface densities of H2, Hi, Hll and the to- 
tal gas (accounting for a 28% contribution from He), 
while Fig. Ogives the vertical variation of their respective 
space-averaged volume densities averaged along the solar 
circle (at R = R®). The total interstellar masses of the 
three gas components in the Galactic disk are uncertain 
and their sum in the Galactic disk, i.e. between 2 and 20 
kpc, is probably comprised between ~ 0.9 10 10 M and 
- 1.5 10 10 M©. 

Since most of the transport and annihilation of 
positrons takes place in the Galactic bulge, its gas distri- 
bution deserves a more detailed description. In the bulge, 
the interstellar gas is roughly equally divided between the 
neutral (molecular + atomic) and ionized components, 
and the neutral component is ~ 90% molecular. The 
molecular gas tends to co ncentrate in the so-called c en- 
tral molecular zone (CMZ; iMorris and Serabvnll 19961 ), a 
thin sheet parallel to the Galactic plane, which, in the 
plane of the sky, extends out to R ~ 250 pc at longitudes 
I > 0° and R ~ 150 pc at I < 0° and has a FWHM thick- 
ness ~ 30 pc. Projected onto the Galactic plane, the 
CMZ would appear as a 500 pc x 200 pc ellipse inclined 
(cloc kwise) by ~ 70° to the line of sight ( Sawada et all 
120041 ). Outside the CMZ, the mole cular gas is con- 
tained in a significantly tilted disk ([Burton and Lisztl 
11992b ILiszt and Burto 

2 H978), extending in the plane 
of the sky out to R ~ 1.3 kpc on each side of the GC 
and having a FWHM t hickness ~ 70 pc. According to 
ILiszt and Burtonl ( 1980f ), the tilted disk has the shape of 
a 3.2 kpc x 1.0 kpc ellipse, which is tilted (counterclock- 
wise) by ~ 13? 5 out of the Galactic plane and inclined 
(near-side down) by ~ 70° to the plane of the sky. The 
tilted disk is also believed to feature an elliptical hole in 
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FIG. 8 Azimuthally- averaged surface densities of interstel- 
lar atomic, molecular and ionized hydrogen as functions of 
Galactic radius. The total gas (bottom) includes a 40% con- 
tribution by helium. Notice the change of scal e at R—2 kpc 
For R <2 kpc (bulge) data are derived by Ferrier e et al 
(2007), based on a compilation of earlier works: ISawada et al 



(2004 ) for the molecula r gas in the Central Molecular Zone, 
iLiszt and Burton] (Il980h for the neutral gas in the tilted disk, 
and lCordes and Lazid (2002) for the i onized gas. In all pan- 
els, disk data (R >2 kpc) are from: iDamd (Il993|) (solid): 
lOlling and Merrifieldl fo OQl) (do tted) ; lNakanishi and Sofud 
(|2006l ) for HI and fNakanishi and Sofuel (120031) fo r H 2 , respec- 
tive ly (dashed) ; and iKa lberla and Dedes (2008) for HI and 
and lPohl et all (|2008T ) for H?< respectively (dot-dashed). The 
curve in the HI I panel is from the N E2001 free-electron den- 
sity model of lCordes and Laziol ([2002 ) (for simplicity, we iden- 
tified the H density with the free-electron density, i.e., we ne- 
glected the contribution of free electrons originating from He 
in the HIM). 



the middle, just large enough to enclose the CMZ. The 
spatial distribution of the ato mic gas is arguably simila r 
to that of t he molec u lar ga s ( Burton and Lisztl (1992), 
but see also ICombesI (| 199 if ) for another point of view) , 
with this difference that the atomic layer is about three 
times thicker than the molecular layer, both in the CMZ 
and in the tilted disk. The ionized gas, for its part, is not 
confined to either the CMZ or the tilted disk; it appears 
to fill the entire bulge and to connect with the ionized 
gas present in the disk. 
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FIG. 9 Space-averaged volume densities of interstellar H2, 
Hi and HII, averaged along the solar circle (R = R@), 
as functi ons of distance f r om t he Galactic plane Z. Data: 
H2 from Bronfma n et all (Il988fh rescaled to Xco = 2.3 x 
10 20 cm " 2 K" 1 km" 1 s, as in foiling and Merrifieldl (|200ll ): 
Hi from iDickev and Lockmanl (|l990f h scaled up by a fac- 
tor 1.58 such as t o mat ch the Hi surface density of 
Oiling and Merrifield (2001). All the other curves are from 
the same sources as in Fig. [8] 



tween the different ISM phases as well as the supersonic 
random motions observed in all of them bear witness to 
a highly turbulent state. Mainly responsible for this tur- 
bulence are the powerful winds and terminal supernova 
explosions of the most massive stars. Interstellar turbu- 
lence manifests itself over a huge range of spatial scales, 
from <; 10 10 cm up to J> 10 20 cm; throughout this range, 
the power spectrum of the free-electron density in the lo- 
cal ISM is consist ent with a Kolmogorov-like power law 
([Armstrong et aZ.L fl995h . 



The dramatic density and temperature contrasts be- 



C. Star formation and supernova rates 

Determination of absolute values of star formation 
rates (SFR, in M yr _1 ) constitutes on e of the most cha l- 
lenging tasks in modern astrophysics ([Kennicuttl Il998). 
In the case of the Milky Way, methods based on counts 
of various short-lived objects (with lifetimes less than a 
few Myr, like e.g. pulsars, SN remnants or OB associa- 
tions) are used. Those methods establish in fact the rela- 
tive star formation rate across the Galactic disk. Surface 
density profiles of various stellar tracers appear in Fig. 
[TUl For the calibration of the SFR profile one needs to 
know either the total SFR of the MW disk or the local 
one in the solar neighborhood. A "ball-park" estimate 
of the former value is obtained by noting that the late 
spectral type (Sbc) of the MW suggests a slow formation 
at a relatively steady rate (SFR) over the past AT —10 
Gyr, leading to (SFR)=M D , thin / AT -2.3 M yr" 1 . 
Most empirical estimates of the present-day total Galac- 
tic SFR, based on the aforementioned tracers (and as- 
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FIG. 10 Surface densities of stars+gas, SFR, SN rates and 
scaleheights of gas and stars as a function of Galactocen- 
tric distance. Star profiles are f r om da ta of Table II and 
gas profile is the one of iDamd (|l993f ) (bottom pane of 
Fig. 6). D ata for SFR are from : iLvne et all (|1985T ) 
(open deles): [Case and Bhattacharval (| 19981 ) (filled circles)', 
iMcKee and Williams! (| 19971 ) (open squares): iGuibert et ali 
(| 19781 ) (filled squares). The solid curve is an approximate 
fit, normalized to 2 M©yr -1 for the whole Galaxy. The same 
curve is used for the CCSN rate profile (third panel), normal- 
ized to 2 CCSN/century; the SNIa rate profile is calculated 
by Eq. [H and normalized to 0.5 SNIa/century (Table |V|) . 



sumptions on the IMF) pr oduce values within a factor 

IMcKee and WilliamsL Il997l 



of two of the (SFR) (e.g. „ „„ 

iRobitaille and Whitnevi l2010l and references therein) . 

In the context of Galactic positrons, special attention 
should be paid to the star formation activity in the cen- 
tral regions of the bulge. The massive star population 
of the three major star clusters inside the Nuclear Bulge 
clearly indicate important recent star formation, obvi- 
ously fed from the gas of the Central Molecular Zone 
(CMZ). Deep field observations of l ate- ty pe stars with the 
NICMOS /H UBBLE riFiger et all I2QQ4D and with SIN- 
FONI/VLT (jManess et all l2QQ7h suggest that the star 
forming activity in that region has proceeded at a rela- 
tively steady rate, of the order of a few 10 ~ 2 M /yr, over 
the past ~10 Gyr. 



From the theoretical point of view, SN are now classi- 
fied mainly in thermonuclear supernovae (the explosion 
energy being due to the thermonuclear disruption of a 
white dwarf accreting matter in a binary system) and 
core collapse supernovae (CCSN, where the energy origi- 
nates from the gravitational collapse of the iron core of a 
massive star having exhausted all its nuclear fuel). Ther- 
monuclear supernovae are identified with SNIa (lacking 
hydrogen in their spectra) and are observed in all types 
of galaxies: old ellipticals with no current star forma- 
tion, but also young, star forming, spiral and irregular 
galaxies. All other supernova types (SNII, SNIb, SNIc) 
are exclusively observed in the star forming regions of 
spirals (i.e. inside spiral arms) and irregulars 8 . 

No supernovae have been observed in the Galaxy in 
the past four centuries, and the handful of so called 
"historical supernovae" offe rs a very biased estima te of 
the Galactic SN frequency ( Tammann et q/| [l994). All 
methods used to determine the Galactic SN rate which 
are based exclusively on Galactic data suffer from vari- 
ous systematic uncertainties and converge to a value of 
Rccsn= a few per century (|Diehl et all I2QQ6I and ref- 
erences therein). The most accurate way to evaluate the 
Galactic SN rate is, probably, through statistics of SN 
rates in external galaxies. The work of iMannucci et all 
(2005), corrected for various observational biases, offers a 
valuable database for such an estimate and can be used, 
along with the stellar masses of the various Galactic com- 
ponents (Sec. III. A) , to derive the Galactic rate of the 
main SN types (Table EJ). 

The spatial distribution of core collapse SN in the MW 
should obviously follow the one of the SFR (Fig. [TO)) . 
Such an azimuthally averaged surface density masks the 
fact that CCSN are exclusively concentrated inside spiral 
arms. The scaleheight of core collapse SN should be com- 
parable to the scaleheight of the molecular gas, i.e. less 
than 100 pc, and little varying with Galactocentric dis- 
tance. More difficult is the evaluation of the radial profile 
of SNIa, since the progenitor white dwarfs may originate 
from stars of a wide variety of stellar masses (1-8 M ) 
and corresponding lifetimes (10-0.05 Gyr). Various mod- 
els have b een deve l oped in order to calculate the SNIa 



rate (e.g. iGreggiol 12005 ). A useful empirical approac h 
is the one adopted in IScannapieco and Bildste nT (|2005h . 



where the SNIa rate is calculated as the sum of two terms: 
one depending on the stellar mass M* and one on the 
SFR of the system, i.e. 



RsNIa 

century 



A 



10 10 M ft 



B 



SFR 

M yr- 



(4) 



with parameters A and B empirically determine d 
([Scannapieco and Bildstenl 120051 : ISullivan et all [2006). 



8 The degree of mass loss suffered by the massive star prior to the 
explosion determines the appearance of the core collapse super- 
nova as SNII (little H lost), SNIb (all H and little He lost) or 
SNIc (all H and most He lost). 
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TABLE V Supernova rates in the Milky Way 









SNIa 


Core collapse SN 




— — — 

Stellar mass 


Spectral type Specific rate 6 


Rate 


Specific rate 6 Rate 




1 nio i\/r 
1U M0 




SNuM 


century 1 


QTVTnl\/r /-./-vn-i-n v»ir — 1 

olNuM century 


Bulge 


1.4 


E0 


0.044 


0.062 




Nuclear Bulge 


0.15 


Sbc/d-Irr c 


0.17-0.77 


0.025-0.115 


0.86-2.24 0.13-0.33 


Thin disk 


2.3 


Sbc 


0.17 


0.4 


0.86 2 


Thick disk 


0.5 


E0 


0.044 


0.022 




Total bulge 


1.5 






0.087-0.18 


0.13-0.33 


Total disk 


2.8 






0.42 


2 


Total Milky Way 


4.3 






0.5-0.6 


2.13-2.33 


Bulge/Disk ratio 


<0.5 






0.21-0.43 


0.06-0.15 



SNuM: = 1 SN per 10 10 M© per century; a: See Sec. III. A for references; b Mannucci et al. (2005); c Very uncertain, in view 
of uncertainties in star formation efficiency and slope of IMF (see text). 



The parametrized SNIa profile in the Milky Way disk ap- 
pears in the third panel of Fig. [101 Taking into account 
the nature of the SNIa progenitors, it is expected that 
the distribution of SNIa vertically to the disk plane will 
follow the corresponding distribution of the thin disk, i.e. 
with a scaleheight of 300 pc (an insignificant contribution 
from the thick disk is also expected). 



D. Interstellar magnetic fields 

The magnetic field strength, £>, in cold, dense regions 
of interstellar space can be inferred from the Zeeman 
splitting of the 21-cm line of Hi (in atomic clouds) and 
centimeter lines of OH and other molecules (in molec- 
ular clouds). In practice, it is the line-of-sight compo- 
nent of the magnetic field, Bu, that is measured. With 
appropriate statistical corrections for projection effects, 
it is found that in atomic clouds, B is typically a few 
/iG, with a slight tendency t o increase with increasing 
density (Heiles and TrolandL 120051 : iTroland and HeilesL 
1986), while in molecular clouds, B increases approxi- 
mately as t he square root of dens ity, from ~ 10 /iG to 
- 3 000 ^G (ICrutcher] . [l999l . l2QQ7D . 

The interstellar magnetic field, B, in the ionized 
medium is generally probed with Faraday rotation mea- 
sures of Galactic pulsars and extragalactic radio sources. 
Here, too, the quantity that is actually measured is Bu. 
Faraday rotation studies have provided the following in- 
teresting pieces of information. 

1) B has a uniform (or regular) component, B u , and 
a random (or turbulent) com ponent, B r . Near the Sun , 
B u ~ 1.5 /iG and B r ~ 5 /iG ([Rand and Kulkarnil . fl989l ) . 
Away from the Sun, B u increases towar d the GC, to 
;> 3 /iG at R = 3 kpc (|Han et all . l2QQ6h . i.e., with an 
exponential scale length <; 7.2 kpc. In addition, B u de- 
creases away from the midplane, albeit at a very uncer- 
tain rate; for reference, the exponential scale height in- 
ferred from t he rotation measures of ex tragalactic sources 
is ~ 1.4 kpc (jlnoue and TabaraL fl98ll ). 



2) In the Galactic disk, B u is nearly horizontal and 
generally dominated by its azimuthal component. It is 
now widely accepted that B u reverses several times with 
decreasing radius, but the number and radi al locations of 
the r e versals are s till highly co n troversial (|Brown et al , 
20071: lHan et all . 120061 Il999fc iRand and Lvnel Il994 : 
Vailed . l2QQ5h . These reversals have often been inter- 
preted as evidence that B u is bisymmetric (azimuthal 
wavenumber m = 1), while an axisymmetric (rn = 0) 
field would be ex pected from dynamo theory. In reality, 
M en et all (|2008l ) showed that neither the axisymmetric 
nor the bisymmetric picture is consistent with the exist- 
ing pulsar rotation measures, and they concluded that 
B u must have a more complex pattern. 

3) In the Galactic halo, B u could h ave a significant ver - 
tical component. For the local halo, iTavlor et all {2009) 
obtained (B u )z ~ —0.14 /iG above the midplane (Z > 0) 
and (B u )z ~ +0.30 /iG below the midplane (Z < 0), 
whereas [Mao et all (|2010h obtained (B u ) z ~ 0.00 /iG to- 
ward the north Galactic pole and (B u )z ~ +0.31 /iG 
toward the south Galactic pole. In contrast to the situ- 
ation in the Galactic disk, the azimuthal component of 
B u shows no sign of reversal with decreasing radius. 

4) At low latitudes (basically, in the disk), B 



appe a rs to be rough l y sym metric in Z (|Frick et al 



l200ll : IRand and Lvnel Il994f ). while at high latitudes 
(in the halo), B u appears to be roughly antisym- 
metric/ symmetric in Z inside/outside the solar circle 
(|Han et all . Il997l Il999h . Finding B u to be symmet- 
ric in the disk and antisymmetric in the inner halo is 
consistent with the predictions of dynamo theory and 
with the results o f galactic dynamo calculation s (e.g., 
iMoss and Sokoloffl [20081 : iRuzmaikin et all Il988f ). How- 
ever, there is no reason to believe that B u is simply a 
pure quadrupole (in the disk) and a pure dipole (in the 
inner halo), sheared out in the azimuthal direction by the 
Galactic differential rotation. In this respect, one should 
emphasize that the picture of an azim u thally -sheared 
pure dipole, originally proposed by iHanl ([20021 ) and of- 
ten used in the cosmic-ray propagation community (e.g., 
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lAlvarez-Muniz et all l2QQ2t |p rouza and SmidaL l2QQ3). is 
supported neither by numerical simulations of galactic 
dynamos nor by observations of external edge-on galax- 
ies, w hich generally reveal X-shaped field patterns ([Becki 

1x1). 

A more global method to map out the spatial dis- 
tribution of B rests on the observed Galactic syn- 
chrotron emission . R elying on the synchrotron map of 
iBeuermann et al\ ( 1985 ) and assuming e quipartit i on be- 
tween magnetic fields and cosmic rays, iFerrierd (|l998h 
found that the total magnetic field has a local value 
^5.1 /iG, a radial scale length ~ 12 kpc, and a local 
vertical scale height ~ 4.5 kpc. Besides, synchrotron po- 
larimetry indicates that the local ratio of ordered (regu- 
lar + anisotropic random) to total fields is ~ 0.6 (|Beckl 
l200lf ), implying an ordered field ~ 3 /iG near the Sun. 

In the vicinity of the GC, the interstellar magnetic 
field has completely different properties from those pre- 
vailing in the Galactic disk. In that region, systems 
of nonthermal radio filaments were discovered, which 
run nearly perpendicular to the Galactic plan e and 
pass through it wi t h litt le or no distortion ([Lisztl Il985b 
lYusef-Zadeh et all Il984f >, The morphology of the fila- 
ments strongly suggests that they follow magnetic field 
lines, and radio polarization measurements (assuming 
synchrotron emission) confirm that the magnet ic field 
in th e filaments is oriente d along their long axis (|Reichl 
I1994J : iTsuboi et all Il985h . From this, it has naturally 
been concluded that the interstellar magnetic field near 
the GC is approximately vertical, at least close to the 
midplane. Farther from the midplane, the filaments tend 
to lean somewhat outwards, consistent with the inter- 
stellar magnet ic field having an overall poloidal geometry 
(lMorrid . [l99Qh . 

The radio filaments have equipartition or 
minimum- energy field strengths ~ 50 — 200 /iG 
( Anantharamaiah et all Il99ll ; lLaRosa et all [2004, and 
references therein). On the other hand, the fact that the 
filaments remain nearly straight all along their length 
suggests that their magnetic pressure is stronger than 
the ram pressure of the ambient interstellar clouds, 
or, equivalent ly, that the i r field strength is B J> 1 mG 
(lYusef-Zadeh and MorrisL Il987l ) . 

Low-frequency radio observations of diffuse nonther- 
mal (supposedly synchrotron) emission from a 6° x 2° re- 
gion centered on the GC imply that the diffuse ISM near 
the GC has a minimum-energy field strength J> 6 /iG 
- possibly up to ~ 80 /iG if the cosmic-ray proton-to- 
electron energy ratio is as high as 100 and the filling 
factor of the synchro tron-emitting gas is as low as 0.01 
(lLaRosa et aZ.L l2QQ5h . A more reliable estimation of the 
general field stren gth in the GC r egion emerges from the 
recent analysis of 1 Crocker et al\ {2010), which combines 
radio and 7-ray data and comes to the conclusion that 
B % 50 M G. 

Far-infrared/submillimeter polarization studies of dust 
thermal emission from the GC region indicate that the 
magnetic field inside GC molecular clouds is roughly par- 



allel to the Galactic plane ( Novak et all 120031 ). More 
precisely, the field direction appears to depend on the 
molecular gas density, being nearly parallel to the plane 
in high-density regions and nearly perpe ndicular to it in 
low-density regions (|Chuss et all l2QQ3h . Near-infrared 
polarization studies of starlight absorption by dust also 
find the magnetic field inside GC molecular clouds 
to be roughly horizontal, although without any obvi- 
ous correlation between field direction and gas density 
(iNishivama et aZ.L l2QQ9h . 



Zeeman splitting measurements have yielded mixed re- 
sults. In the circumnuclear dis k, the innerm o st mo lecu- 
lar region wi t h ra dius <; 7 pc, iKilleen et all ([l992) and 
iPlante et al\ (|l995h derived line-of-sight magnetic fields 
|£||| ~ 2 mG and \Bu \ ~ 0.6 - 4.7 mG, r espectively. 
Farther from the GC, ICrutcher et all ( 19961 ) measured 
values of \Bu\ ranging between ~ 0.1 and 0.8 mG to- 
ward the Main and North cores of Sgr B2. In con- 
trast, lUchida and Guestenl (|l995) reported only non- 
detections, with 3<j upper limits to \Bu\ ~ 0.1 — 1 mG, 
toward 13 selected positions within a few degrees of the 
GC (including Sgr B2). 

Faraday rotation measures have also yielded somewhat 
disparate results. The disparity lies not so much in the 
absolute value of B\\ , which is generally estima t ed at 
a few /iG (e.g., \Qi^_et all Il99ll ; ITsuboi et all Il985l : 



lYusef-Zadeh and MorrisL Il987h. bu t mor e in the (/, b)- 
dependence of its sign. " iNovak et alJ(|2003h . who collected 
all the available rotation measures toward synchrotron 
sources within 1° of the GC, found that Bn reverses sign 
both across the rotation axis and ac ross the midplane . 
A different pattern was uncovered by lR,ov et al. (2005), 
who derived the rotation measures of 60 background ex- 
tragalactic sources through the region < 6°, |6| < 2°) 
and obtained mostly positive values, with no evidence for 
a sign reversal either with / or with b. 

The properties of the turbulent ma g netic field are not 
well established. [Rand and Kulkarnil (|l989h provided a 
first rough estimate for the typical spatial scale of mag- 
netic fluctuations, ~ 55 pc, although they recognized 
that the turbulen t field cannot be characterize d by a sin- 
gle scale. Later, iMinter and Spanglerl (|l996h presented 
a careful derivation of the power spectrum of magnetic 
fluctuations over the spatial range ~ (0.01 — 100) pc; 
they obtained a Kolmogorov spectrum below ~ 4 pc and 
a flatter spectrum consistent with 2D turbulence above 
this scale. In a complementary study, iHan et all ( 20041 ) 
examined magnetic fluctuations at larger scales, rang- 
ing from ~ 0.5 to 15 kpc; at these scales, they found 
a nearly flat magnetic spectrum, with a ID power-law 
index ~ —0.37. 

The, poorly understood at present, properties of the 
turbulent Galactic magnetic field as well as its overall 
configuration, are extremely important for understanding 
positron propagation in the Milky Way (Sec. VI). 
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FIG. 11 Composite magnetic energy spectrum in our Galaxy. 
The thick solid line is the large-scale spectrum. The thin 
solid and dashed/dotted lines give the Kolmogorov and two- 
dimensional turbulence spectra, respectively, inferred from 
the Minter and Spangler (1996) study. The two-dimensional 
turbulence spectrum is uncertain; it probably lies between 
the d ashed \E B (k) oc k ~ 2/s ] and dotted [E B {k) oc k~ 5/s ] lines 
(from lHan etall\2004 \. 
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FIG. 12 Dark matter density profile (top) and rotational ve- 
locity (bottom) of the Milky Way; the various components 
(bulge, stellar disk, gas and dark halo) contributing to the 
latter are also indicated. In both panels thick and thin 
solid curves correspond to NFW and isothermal ( "ISO") dark 
halo profiles, respectively. Data points are from ISofu7 et al. 
(|2QQ8h . 



E. The dark matter halo 



A large body of observational data on the extragalac- 
tic Universe suggests that its mass is dominated by non- 
baryonic dark matter. In the presently widely accepted 
"standard" cosmological model (ACDM, for Cold Dark 
Matter with cosmological constant A) dark matter ac- 
counts for a fraction VLdm ~24% of the overall mat- 
ter/energy budget of the Universe, baryons for ~4 % and 
dark e nergy - or cosmologi cal constant - for the remaining 
~72% (|Bartelmannl . l2QQ9l ). 

The presence of dark matter in spiral galaxies is de- 
duced from the fact that their rotation curves beyond 
a radius of ~ 3 scalelengths do not fall off as rapidly 
as expected from their baryonic content. In the case of 
the Milky Way, the rotation curve is poorly determined 
beyond the Sun's location (R =8 kpc). It is then as- 
sumed, rather than directly inferred from observations, 
that the MW is found inside a dark matter halo with 
a density profile Pdm(t) similar to those found in nu- 
merical simula tions of structure for mation in a ACDM 
universe (e.g. iNavarro et all Il997l ). In the absence of 
baryons such simulations predict approximately univer- 



sal density profiles Pdm(t) oc 



with k being itself 



a positive function of radius r: k(r) oc r s ("Einasto 
profile"). Because of finite numerical resolution, val- 
ues of k cannot yet be reliably determined in the inner 
halo. Some simulations find fc=1.5 in the inner galaxy 



( Moore et a/.Ul999l ) , but the analys i s of on e of the largest 
simulations so far (|Navarro et all 120081 ) suggests that 
&=0.9±0.1, i.e. a value compati ble with the value o f 
k=l in the classical NFW profile (iNavarro et a/| [l997). 
For values of k >1 mass diverges as r — >0 (cuspy pro- 
f iles). Including interactio ns of dark matter with baryons 
Blu menthal et aL 1 19861) or with a central black hole 
Gondolo and SilkL 1 19991 ) generically tend to enhance the 



(IGondolo and Silk. 199 9]) generically tend 
cusp (e.g. ISellwood and McGaughL 120051 ). 



The shape of the dark matter density profile in the 
inner Galaxy is obviously crucial for the corresponding 
profile of the putative positrons released from dark mat- 
ter decay, annihilation or de-excitation (see Sec. IIV.C|) . 
Since dark matter is sub-dominant in the inner Galaxy 
(see Fig. [T2]) , observations of the rotation curve can- 
not help to determine its density profile. Analyzing 
observations of the optical depth of the inner Galaxy 
to microlensing events (which are a ffecte d only by the 
baryonic mater) iBinnev and Evansl ( 200lf ) find k ~0.3. 
On the other hand, rotation curves of dwarf galaxies 
(which are dominate d by dark matter) s ystematically 
sugge st flat profiles ([Gentile et all 120071 ; ISpano et all 
2008) with k ~0, such as those obtained in the ca se of 
cored isothermal dark halo (see also iMerrittl l2010f ). A 
useful parametrization of the density profiles is 
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TABLE VI Models for the Milky Way dark halo profile. 



DM profile: 


ISO 


BE 


NFW 


M99 


a 


2 


1 


1 


1.5 


P 


2 


3 


3 


3 


7 





0.3 


1 


1.5 


r (kpc) 


5 


10 


20 


30 


po (M pc- 3 ) 


5. 10" 2 


7. 10" 2 


1.0 10" 2 


1.7 10" 3 


po (GeV cm- 3 ) 


1.89 


2.65 


0.38 


0.065 



ISO: Isothermal; BE: Binney and Evans (2001); NFW 
(Navarro et al. 1997); M99 (Moore et al. 1999) 



p(r) = 



(5) 



where po and ro are, respectively, the characteristic 
mass/energy density and radius of the halo and a, /3, 7 
are parameters with values (found, either from simula- 
tions or from observations), reported in Table fVTl 

The shape of the dark halo profile may deviate from 
spherical symmetry. A triaxial shape arises naturally 
from the fact that gravitational collapse of the halo starts 
first (and proceeds more rapidly) in one direction. How- 
ever, ot her processes may subsequently erase it (e.g. gas 
cooling, iKazantzidis et all l2004f ). Various observations 
in the Milky W ay have been interpret e d as suggest- 
i ng a spherical jM artinez-Delgad o et al. oblate 
([Martinez-Delgado et all . l2004f ) or prolate ([Helmil . l2004f ) 
dark halo, but in any case, deviations from spherical sym- 
metry appear to be small. 

Structure formation in the ACDM model leads to a hi- 
erarchy of dark haloes embedded within the main halo 
of a galaxy. Since smaller galaxies are more dark matter 
dominated than larger ones, the strongest signal from 
dark matter annihilation may not arise from the main 
halo, but from satellite galaxies. This important is- 
sue has been extensively studied recently. Analyzing 
one of the largest "M ilky Way size" simulations so far, 
ISpringel et all (|2008l ) find that the most intense emission 
is expected to arise from the main halo. We shall further 
discuss this point in Sec. IV. D. 3. 



IV. POSITRON PRODUCTION: PROCESSES AND 
SOURCES 

A. Radioactivity from stellar nucleosynthesis 

1. Radioactivity 

Positrons are emitted by the /3 + -decay of unstable nu- 
clei which turns a proton into a neutron, provided the 
mass difference between parent and daughter nucleus is 
AM > m e c 2 (where m e is electron's mass and c the light 
velocity). /3 + -decay of unstable nuclei produced in stellar 



explosions was one of the earliest candi dates proposed to 
explain the Galactic 511 keV emission (|Clavtonl . ll973h . 

Astrophysically important e + emitters are produced 
in proton-rich environments, either hydrostatically (e.g. 
in massive star cores) or explosively (in novae or SN 
explosions); in both cases, proton captures occur on 
shorter timescales than the corresponding lifetimes of 
decaying nuclei along the nucleosynthesis path. Impor- 
tant e + emitters are also found in the Fe-peak region and 
are produced in the so-called Nuclear Statistical Equilib- 
rium (NSE) regime, at temperatures T > 4 10 9 K. In the 
short timescale of the explosion (r ~1 s in SNIa and in 
the inner layers of CCSN) weak interactions can hardly 
operate and material is nuclearly processed under the 
effect of strong interactions alone, moving in the neu- 
tron vs. proton (N — Z) plane along a N/Z ~ const 
trajectory. The original stellar material is essentially 
composed either of 28 Si (N = Z = 14) in the case of 
CCSN or of 12 C (TV = Z = 6) and 16 (N = Z = 8) 
in the case of the white dwarf progenitors of SNIa; this 
N/Z ~ 1 ratio is mostly preserved during the explo- 



sion. Since the last stable nucleus with N 
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Ca 



(N = Z = 20), NSE reactions produce mostly unsta- 
ble Fe-peak nuclei, which decay later back to the nuclear 
stability valley by electron captures (EC) or e + emission. 
This is typically the case of the most abundant Fe-peak 
nucleus 56 Fe (Z = 26, N = 30), which is produced as 
56 Ni (Z = N = 28) through the decay chain 56 Ni — > 
56 Co — > 56 Fe; the first decay proceeds by EC and the 
second one by both EC and e + emiss ion, with branchin g 
ratios of 81% and 19%, respectively (jNadvozhinl . Il994t ). 

Other important astrophysical e + emitters are dis- 
played in Table IVII1 along with various relevant data. 
An important feature of f3 + -decay is that positrons are 
released with energies in the MeV range, i.e. they nat- 
urally satisfy the constraint imposed by the continuum 
observations of the inner Galaxy in that energy range 
(see Sec. II.B.l). 

Contrary to all other e + sources presented in this sec- 
tion, it is well established that stellar radioactivities con- 
tribute at a non-negligible level to the e + production rate, 
because of the observed presence of 26 Al in the disk (see 
Sec. II. B. 2). The uncertainties related to their overall 
contribution stem from two factors: 

i) In the case of short-lived radioactivities (i.e. with 
lifetimes short compared to the characteristic timescales 
of SN expansion) positrons are released in high density 
environments and in magnetic fields of unknown config- 
uration. Those conditions render difficult the evaluation 
of the fraction of e + escaping to environments of suffi- 
ciently low density for their annihilation photons to be 
detectable. This is the case of 56 Co in SNIa. 

ii) In some cases of long-live radioactivities, the cor- 
responding stellar yields and/or the frequencies of the 
nucleosynthesis sites are quite uncertain. Indirect meth- 
ods should then be used to evaluate their contribution to 
the Galactic e + production. This is the case of 44 Ti. 
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TABLE VII Astrophysically important positron-emitting radioactivities 



Nuclide 


Decay chain 


Decay mode Lifetime 
and e + BR a 


Associated 7-ray lines 
Energy in keV (BIT) 


Endpoint e + Mean e + 
energy (keV) energy (keV) 


Sources 


56 Ni 


56 Ni - 
56 Co- 


-> 56 Co* 
-+ 56 Fe* 


EC b 
e + (0.19) 


6.073 d 
77.2 d 


158(0.99), 812(0.86) 
9^<Wf) 1 7"! 1 771 (0 1 ^ 


1458.9 


610 


SNIa 


22 Na 


22 Na - 


22 Ne* 


e+ (0.90) 


2.61 y 


1275(1) 


1820.2 


215.9 


Novae 


44 Ti 


44 Ti _ 

44 Sc- 


-+ 44 Sc* 
-> 44 Ca* 


EC b 
e + (0.94) 


59.0 y 
3.97 h 


68(0.94), 78(0.96) 
1157(1) 


1474.2 


632. 


Super novae 


26 A1 


26 A1 _ 


-> 26 Mg* 


e + (0.82) 


7.4 10 5 y 


1809(1) 


1117.35 


543.3 


Massive stars 



(a) BRiBranching Ratio (in parenthesis); (b) EC: Electron capture 



2. Massive stars: 26 AI and 44 Ti 

The observed irregularities in the 26 Al 7-ray emis- 
sion along the pl ane of the Galax y as shown in the 
COMPTEL map (jDiehl et all Il995h suggest that mas- 
sive stars are the dominant source, as these are the 
only candidate 26 Al sour c es c lustered along spiral arms 
(see iPrantzos and DiehlL I1996L and Sec. II.B). 26 Al 
is produced in such stars both hydrostatically (during 
H-burning) and explosively (in the C-Ne-O) layers; it 
is ejected by the Wolf-Rayet stellar winds in the for- 
mer case and by th e super nova explosion in the latter. 
iLimongi and Chieffil (|2006l ) find that in their Z =Z 
models explosive nucleosynthesis is always dominant; 
however, models with rotation and at Z >Z (appro- 
priate for the inner Galaxy) may modify this conclusion 
somewhat. Stellar nucleosynthesis models find typical 
yields of ~10 -4 M of 26 Al per star, which combined 
with the derived CCSN frequency in the Galaxy (Table 
|V]) results in a production rate comparable to the ob- 
served one of ~2.7M /Myr; thus, the nucleosynthesis of 
26 Al is considered to be rather well understood quan- 
titatively (within a factor of 2). Independently of the- 
oretical considerations, however, the observed Galactic 
decay rate of 26 Al corresponds to a production rate of 



e+,26 



4 10 
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in the Galactic disk. 



44 Ti is produced in the innermost layers of th e super- 
nova , in the "q-rich freeze-ou t" regime of NSE (|Meverl 



1 19931 [Thielemann et a/| [l996). Its yields are much more 
uncertain than the ones of 26 Al because of uncertain- 
ties either in the nucl ear reaction rates (which affect its 
yield s by a factor of 2 , iThe et a/1 120061 : iMagkotsios et all 
120081 ) or, most importantly, in the explosion mechanism 
itself ([Wooslev and Weaver! Il995fc iTimmes et a/1 Il996h . 
Moreover, asphericity effects (due e.g. to rotation) ap- 
pear to be critical, leading to the production of substan- 
tially higher 44 Ti yields (and 44 Ti/ 56 Ni r atios) than in the 
case of spherically symmetric models ( Nagataki et all 
Il998h , 

Observations offer little help in this case. 44 Ti 



has been directly detected in the Cassiopeia A 
(CasA) SN remnant, th rough its 7-ray li nes, both 
with COMPTEL / CGRO (llvudin et all \l994 ) and with 



SPI/INTEGRAL (jRenaud et all 120061 ). Its presence is 
also indirectly derived in SN1987A, the closest observed 
supernova in the past four cent uries, since it is required 
to ex plain the late lightcurve ([Motizuki and Kumagaii 
\2004 . In both cases the derived 44 Ti yield is F 4 4 ~ 
2 10 -4 M , substantially larger than predictions of spher- 
ically symmetric models, but comparable to predictions 
of aspherical models. Asphericity is also favoured for 
CasA and SN1987A on the basis of other observables 
([Prantzosl 12004 and references therein). Does this mean 
that typical sources of 44 Ti are aspherical and have the 
aforementioned yield? 

CasA is found at a distance of ~3 kpc from the Earth 
in the outer Galaxy (outside the active star forming re- 
gions of the inner Galaxy) and its age is estimated to 300 
yrs (much larger than the 44 Ti lifetime). That a super- 
nova with such properties is the only one detected so far 
through its 44 Ti lines, despite the sensitivity of COMP- 
TEL/CGRO and SPI/INTEGRAL Gal actic surveys, ap- 
pears to be s t atisti cally improbable ([The et all i2006; 
Ren aud et al. , 2006). It may imply that typical 44 Ti 
sources are rare, i.e. with frequencies much lower than 
the CCSN frequencies of Table [Vj and, consequently, 
much larger yields. Sub-Chandrasekhar mass SNIa (i.e. 
thermonuclear SN induced by surface He-detonation, see 
Sec. IV. 4) are potential candidates, since they produce 
10-20 ti mes more 44 Ti than a typ ical massive star ex- 
plosion ([Wooslev and Weaverl Il994l ); but, provided that 
such objects exist and have the required yields, their fre- 
quencies are totally unknown. 

In those conditions, the only way to evaluate the Galac- 
tic 44 Ti production rate is through a nucleosynthesis ar- 
gume nt, based on i) the solar ( 44 Ca/ 56 Fe) = 1.2 10~ 3 
ratio ([Loddersl |2003) , i.e the ratio of the stable prod- 
ucts of 44 Ti and 56 Ni decays and ii) the knowledge of 
the current production rate of 56 Fe, based on disk SN 
frequencies of Table [V] and on presumably well-known 
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typical yields of 56 Fe: Y 5 s 6 NIa -0.7 M (see Sec. IV.A.4) 
and Y£j CS N —0.07 (from the observed lightcurve of 
SN1987A. lArnett et all Il989h . The production rate of 
44 Ti is then: 



^44 = (^56^j ( R SNIaY56 NIa 



RccsnY 5 c 6 csn ) (6) 



and the corresponding e + production rate is iV e + 5 44 — 3 x 
10 42 e + s _1 , i.e. comparable to the one of 26 Al. Thus, the 
two long-lived radioactivities together may account for 
most, if not all, of the disk production rate of positrons, 
as revealed by the SPI/INTEGRAL analysis. The same 
analysis, applied to the bulge (and assuming the bulge 
44 Ca/ 56 Fe ratio to be solar) leads to a e + production 
rate three times smaller, i.e. an insignificant fraction of 
the obervationally required rate for that region. 



3. Hypernovae and 7-ray bursts 

Hypernovae are very energetic supernova explosions, 
with typical observed kinetic energies >10 52 ergs (i.e. 
about ten times larger than normal su pernovae) and 
ejecte d 56 Ni masses of —0.5 M (e.g. INomoto et all 
l2007f ). Their properties are usually interpreted in terms 
of aspherical explosions of rotating massive stars (with 
mass >30 M ). The rotating Fe core implodes to a 
black hole, around which the surrounding material forms 
a short-lived (—0.1 s) accretion disk. The gravitational 
energy of accretion is partially transferred (by some still 
unclear mechanism) to two jets along the rotation axis, 
which launch the supernova explosion. Heavy nuclei 
(among which 56 Ni) are formed in the hot basis of the 
jet and ejected in the ISM. This model was originally 
proposed to account for the phenomenon of Gamma- Ray 
Bursts (GRB), the most powerful electromagnetic bea- 
cons in the Universe, releasing — 10 51 erg in short flashes 
of 7-rays bea med along the jet direction (the "collap- 
sar" model of IWooslevl 1 19931 ). Observed metallicities of 
GRB host galaxies are ty pically few times lower than so- 
lar ([Savaglio et a/ll2008f ): such low metallicities prevent 
substantial losses of mass and angular momentum and 
allow for a rapid rotation of the core at the moment of 
the explosion, a crucial ingredient of the collapsar model. 

Hypernovae/GRBs have been suggested as potential 
sources of the Galactic positrons, produced either from 
the 56 Ni decay ([Nomoto et all l2QQlfc ICasse et all \2004 
or from pair creation, as photons backscattered from the 
ionized mediu m ahead of the jet interact with the GRB 
7-ray photons ([Parizot et all [20051 : iBertone et a/l l2006). 
Because of the complex (and still very uncertain) nature 
of those objects, the corresponding positron yield is vir- 
tually unknown. In the light of the observational (and 
theoretically motivated) constraint of low metallicity for 
the progenitor stars, the existence of such objects in the 
metal-rich bulge (see Sec. III.A.l) should be excluded. 
Besides, a small bulge/disk ratio would be logically ex- 
pected in that case, contrary to observations. 



4. Thermonuclear supernovae (SNIa) 

SNIa display a remarkable uniformity in their proper- 
ties, like e.g. the peak luminosity, which is attri buted to 
the p ower input of —0.7 M of radioactive 56 Ni ([ArnettL 
1982) 9 . There is general agreement that SNIa result from 
the thermonuclear disruption of a white dwarf, igniting 
explosively its carbon. The thermonuclear flame may 
propagate either subsonically (deflagration) or superson- 
i cally (detonation) inside the white dwarf; Maz zali et al\ 
( 20071 ) show that the SNIa variety can be understood 
within a single, combined, model, involving both defla- 
gration and detonation. There are two main scenarios for 
the precursors of SNIa: the single degenerate (SD) model, 
in which accretio n is made from a main s equence or red 
giant companion ( Wh elan and IbenL Il973h ; and the dou- 
ble degenerate (DD) model, which involves th e merging 
of tw o white dwarfs in a c lose b ina ry system ([ Webbink, 
Il984t llben and Tutukovl Il984l ). iParthasarathv et ali 
(|2007l ) discuss all available observational evidence and 
find that the SD channel i s by f ar the dominant one (but 
see iGilfanov and Bogdanl l2010l for a different view) . 

Most studies of SNIa were made in the framework of 
the SD scenario and, up to the late 1990s with one- 
dimensional (ID) models. Detailed ID models exploring 
the various possibilities (and the corresponding parame- 
ter space) have been developed over the years. Perhaps 
the most successfu l ID model developed so far is the so- 



called W7 model (INomoto et all 1 19841), the phy sics of 



which has been updated in llwamoto et all (|l999h ; it is 
a deflagration model producing in its inner layers —0.7 
M of 56 Ni and negligible amounts of other positron 
emitters. A more accurate description of reality is pur- 
sued by the upcoming; generation of multi-dimensional 



models dTravadio et all 12004 


, Bravo and Garcfa-Senz, 


2006, 1 Schmidt and Niemever, 


2006L iRouke et all 12007. 


Rouke and Niemever. l2007h. 


Preliminary results show 



interesting features for the stratification of radioac- 
tivities, in particular the presence of substantial 56 Ni 
amounts within outer, high velocity, layers (Fig. [T3|) . 

The fate of the /3 + -decay products (7-rays and 
positrons) in the expandin g SNIa eject a has b e en ex- 
tensively studied in ID ([Gomez- Gomar et 

ID, [1998, 



Milne et al. 



( Isern et 



. ]2004ft . and more recently, in 3D models 
. l2007l ISim and Mazzali l2QQ8h . Generically, 
before peak luminosity, the SNIa envelope is opaque, and 
both the energy of the explosion and of /3 + -decay are de- 
posited in and diffuse outwards through the ejecta. After 



9 In fact, the 56 Ni mass m ay vary by a factor of ~10, as shown by 
Stritzinge r et al\ |2006), who find values in the range of 0.1-0.9 
for a sample of seventeen well observed SNIa. However, ober- 
vations indicate that the shape of the SNIa light curves is asso- 
ciated to the 56 Ni mass (with brighter SNI a fading more s lowly) 
and after correction is made for that effect (Phillips, 1993) SNIa 
can indeed be used as "standar d candles" for th e determination 
of cosmological distances (see Leibundgut, 2001 for a review). 
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v [cm/s] 

FIG. 13 Mass fraction of 56 Ni as a function of eject a veloc- 
ity after a SNIa explosi on. The sp ace- averaged profile of 3-D 
mode (solid curve from iRopke et all (|20Q7h ) is compared to 
observational data (dotted curve). Both theory and observa- 
tions find non-negligible am ounts of 56 Ni in the outer (high 
velocity) ejecta. Figure from lRopke et al] (|2007[ ). 
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FIG. 14 Distribution of emitted positron kinetic energies 
as estimated by Segre (1977, dashed curve) compared to the 
spectrum of escaping positrons from W7, as estimated by 
Chan and Lingenfelter (1993, solid curve) and Milne et al. 
(1999, solid histogram). The slowing of the positrons leads to 
the mean energy shifting from 632 to 494 keV . Figure from 
Milne et al. (1999). 



the peak of the bolometric lightcurve (~20 days after the 
explosion) the luminosity evolves from radioactive energy 
deposits and increasing energy leakage in a way follow- 
ing (surprisingly closely, given the interplay of these pro- 
cesses) the decay rate of 56 Co. About 6 months later, 
the ejecta are completely transparent to 7-rays and the 
SNIa luminosity results almost exclusively from energy 
deposited by positrons from ongoing radioactive decays. 
If positrons are trapped (escape) a flattening (steepening) 
of the light curve results. How many positrons ultimately 
escape to the ISM depends on the distribution of the par- 
ent radioactivities within the supernova, the evolution of 
its density, temperature and ionization profiles and, most 
importantly for the late phases, on the unknown config- 
uration of its magnetic field. P rogenitor white dwarfs 
have field strengths of 10 5 -10 9 G. IChan and Lin genfelter 
( 19931 ) found that, in the case of radially combed mag- 
netic fields and fully mixed ejecta a substantial frac- 
tion of 56 Co po sitrons (>1Q%) ma y escape. Building on 
the same ideas, iMilne et all (1999) compared SNIa mod- 
els to observations of late lightcurves of a dozen SNIa 
(mostly in B and V bands) and concluded that, typi- 
cally, a few % of positrons finally escape the ejecta; the 
average positronic yield of a SNIa is n e +(SNIa)~8 10 52 
(corresponding to an escape fraction of f esc ~0.03). They 
also found that the mean energy of escaping positrons is 
-0.5 MeV (Fig. EH}. 

The corresponding Galactic positron yield is then es- 
timated as 



N e +,SNIa = n e +(SNIa) RsNIa 



1.6 io 4 ^- 1 (7) 



but the bulge/disk positron emissivity ratio is B/D— 0.4, 
considerably less than derived from observations. 

This simplified picture may not apply to SNIa in gen- 
eral, though. Bolometric observations (including the near 
IR) of the late lightcurves of SN 2000cx (a rather peculiar 
at early times SNIa) and of SN 2001el and SN2003hv (two 
typical SNIa), interpreted in the fr amework of ID models, 
suggest that no positrons e scape (Sollerma n et a/ll2004 
IStritzinger and Sollermanl 120071 Leloudas et all l2QQ9h : 
in that case, despite their large /3 + -decay yields, the SNIa 
would be insignificant e + producers. 

However, 3D effects may considerably alter the stratifi- 
cation of radioactivities inside the SNIa ( Blinnikov et all 
120061 ), allowing for substantial amounts of 56 Ni to be 
mixed out to the surface (Fig. [13]) and for the re- 
leased positrons to escape at early times (i.e. when 
the lightcurve is dominated by 7-rays, and not yet e + 
deposition) without being noticed. Studying the very 



where RsNia is the SNIa frequency in Table [Vj The 
total e + yield is comparable to the observed Galactic one, 



deposition; witnout being noticed, stud ying tne very 
early optical spectra of six SNIa, iTanaka et al\ (2008) 
find indeed indications for asphericity and substantial 
amounts of 56 Ni present in the high velocity ejecta (v ~10 
000 - 15 000 km/s). Positrons produced by the subse- 
quent decay of 56 Co may escape the ejecta if the mag- 
netic field of the supernova is radially combed. Obser- 
vations indicate that this may be the dominant config- 
uration of magnetic fields in young SN r emnants (e.g. 
iMilne et all . Il993fc iKothes and Reich! l200l[ ) although the 
origin of such a configuration is not yet clearly u nder- 
stood (|jun and Normanl Il996l ISchure et all [2008). On 
the other hand, the late lightcurve of SNIa may also 
be (at least partially) powered by internal conversion 
and Auger electrons re leased from the decay of 57 Co 
(Seitenzahl et al I l2QQ9h . thus allowing for some Co 
positrons to escape. Thus, the issue of the positron yield 
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of SNIa is no t settled yet : n P + may well be as high as 
envisioned by Milne et all (|l999h (albeit for different rea- 
sons), but also much lower. 

When the SN becomes sufficiently diluted, the an- 
nihilation 7-ray ph otons may be directly observed. 
Kale mci et al. I d2QQ6h find no such signal in observations 
with SPI/INTEGRAL of the SN remnant SN1006. They 
exclude then SNIa as major e + producers in the Galaxy 



under the assumption that the e + lifetime is r e + <10 5 yr. l__ 



However, even the low energy positrons of /3 + -decay may 
live much longer before annihilation and then escape the 
SN remnant, especially in the case of a radially combed 
magnetic field. 



5. Novae 

Novae result from explosive H-burning on the surfaces 
of white dwarfs in binary systems. Accretion of material 
from the companion star increases the density, pressure 
and temperature at the base of the white dwarf envelope, 
up to the point where hydrogen ignites in degenerate 
conditions and burns explosively at peak temperatures 
of several 10 8 K. The ejected mass is ~10 -4 M and 
is substantially enriched with the material of the white 
dwar f, leading to CO or to ONe novae (see, e.g. llHernanzl 
2005 for a review). 

Major positron producers are 13 N and 18 F (produced 
in the hot-CNO cycle) and, in the case of ONe novae, 
22 Na (produced in the hot-NeNa cycle). The short life- 
times of 13 N and 18 F (r=862 s and 158 min, respectively) 
make unlikely a substantial escape of their positrons 
from the nova eject a. Positrons from 22 Na decay cer- 
tainly escape and recent calculations suggest that up to 
10 -8 M^ of this nucleus may be produced in ONe novae 
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( Hernanz and Josel l2006h , releasing up to n e + 

e+ 

The novae freque ncy in the Galaxy is estim a ted t o 
Rnova -35 yr" 1 ( (|Shafterl Il997l iDarnlev et oil l2006f ). 
About 1 /3 of those ma y originate from ONe white dwarfs 
(iGil-Pons et all 12003). leading to a Galactic e + produc- 
tion rate of N e +,Novae = Rnov a n e + j7lova ~1.5 10 41 s" 1 , 
ie. smaller by two orders of magnitude than the observed 
rate in the bulge or in the disk. It should be noticed that 
ONe novae appear mostly close to the Galactic plane 
(jdella Valle and Livid . fl998l ). 



B. High energy processes in cosmic rays and compact 
objects 

1. High energy processes 

la. Inelastic p-p collisions 

Relativistic protons and heavier nuclei are present in 
many astrophysical environments in the Galaxy. Their 
inelastic interactions with interstellar gas produce sec- 
ondary particles including numerous neutral and charged 
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FIG. 15 The energy spectra of positrons from the decay of 
7r + produced in collisions of isotropic monoenergetic protons 
with protons at rest for various proton kinetic energies (from 
bottom to top): 0.316 . 0.383, 0.464 0.562 , 0.681, 1.0, 1.78, 
3.16, 10.0, 100.0 GeV ((Murphy et aZlll987h . 



pions and kaons pp — » tt + X, pp — » K + X. In 
turn, decay of positively charged mesons produces sec- 
ondary positrons. The dominant channel is pion decay 



/i + — >• Vpif e t ? , though a non- negligible con- 
tribution comes from the charged kaon decays. The two 
main kaon decay modes contributing to the secondary 
e ± spectrum are K ± — >> fiu^ (63.5%) and if ± — >• 7r°7r ± 
(21.2%). The processes as the source of secondary cos- 
mic ray positrons and dif fuse 7 ray emission have been 
thoroughly studied (e.g., [Strong et all 120071 and refer- 
ences therein). A review of the experimental data for 
pion production in proton-proton collisions and relevant 
cr oss section parameteri zations <50 GeV were presented 
bv lBlattnig et al\ ( 20001 ). New parameterizations of neu- 
tral and charged pion cross sections which provide an ac- 
curate description of the experimental data in a wide en- 
ergy range from the p ion production thre s hold up to 10 5 
TeV are dis cussed in iKelner et al\ (|2006l ) , iKamae et all 
(I20Q6L l2QQ7h . 

The energy spectra of positrons from the decay of 7r + 
mesons produced in collisions of isotropic monoenergetic 
protons with protons at rest are shown in Fig. [15J they 
typically present a maximum at E ~ 30 — 40 MeV. 

lb. 7 — 7 pair production 

Positrons can also be produced in photon-photon 
interactions when the product of their energies is > 
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FIG. 16 The total cross section of 7 + 7 — » e + + e~ reac- 
tion as a function of the Lorentz factor of created positron in 
the center of mass frame (solid). The total cross section a a 
for two-photon pair annihilation reaction e + + e~ —^7 + 7 
(unpolarized) is indicated as a dotted curve. 



2m 2 c 4 /(l — cos#) where is the angle between the pho- 
ton directions. The total unpolarized cross section for 
the creation of e ± by two photons 77 — ^ e + e~ can be 
expressed as a function of a dimensionless velocity j3 of 
the produced particles in the center of mass frame: 



(3-/3 2 )lnl±|-2/?(2-/? 2 ) 



(8) 

where <jt ~ 6.65 10 25 cm 2 is the Thomson scatteri ng 
cross section (Breit and Wheeled . Il934t iGreinerL l2Q()l . 



In Fig. [16] the cross section is presented as a func- 
tion of the Lorentz factor of created positron in the 
center of mass frame. The positron production due to 
photon-photon collisions is suppressed at the threshold 
and reaches a maximum at j3 ~ 0.7. 

The production of a e + e~-pair by a single photon 
is possible in magnetic fields B J>10 12 G observed in 
highly magnetized object s such a s pulsars and ma gnetars 
(jDaughertv and Hardingl . 11983c iKlepikovl Il954f ). This 
occurs with significant probability when the photon en- 
ergy is ~ 3/ (#12 sin 0) MeV, where B12 is the external 
magnetic field strength in units 10 12 Gauss and is the 
angle between the photon direction and the magnetic 
field. 



2. Galactic cosmic rays 

The majority of positrons in cosmic rays (CR) are 
believed to be secondaries produced by interactions of 
relativist ic particles with interstellar gas; however re- 
cent measurements of posi tron fraction i n cos mic rays 
e+/(e-+e+) by PAMELA (jAdriani et alll2009h indicate 
that there may be another component at high energies. If 
produced by CR interactions, t he positron fraction is ex- 
pecte d to decrease with energy (jMoskalenko and Strong) . 
Il998h while the PAMELA data show it rises above ~10 
GeV. The origin of this additional component is inten- 
sively debated. The ideas proposed can be roughly di- 
vided into two broad classes: con ventional so urces, such 
as SNR or nearby pulsars (e.g. iBlasil . 120091 and refer- 
ences therein) and exotic sources such as W IMP annihila- 
tion or decay (e.g. lArvanitaki et all 120091 and references 
therein). The predicted behaviour of the positron frac- 
tion at energies higher than currently measured (~100 
GeV) depends on the model and can be used to distin- 
guish between different possibilities. 

In this section we will discuss positrons produced by 
conventional CR interactions with interstellar gas. The 
most important factors are the energetics of cosmic rays 
and their diffusi on in the int e rstell ar medium (for more 
details see, e.g.. IStrong et aZJ , l2Q07h . 

The major CR sources are believed to be supernovae 
(SN) and their remnants (SNR) with some fraction com- 
ing from pulsars, compact objects in close binary sys- 
tems, and stellar winds. Recen t observations of X-ray 
and 7-ray emission from SNRs ([Aharonian et all 120061 : 
Pan nuti et al. L l20Q3h reveal the presence of energetic elec- 
trons, thus testifying to efficient acceleration processes 
near these objects. The total power of Galactic CR 
sources needed to sustain the observed CR density is es- 
timated at ~10 41 erg s _1 which corresponds to approxi- 
mately 10 50 erg per SN, if the SN rate in the Galaxy is 
1 every 30 years (Table [Vj. This value is ~10% of the 
corresponding total kinetic power of the SN ejecta, an ef- 
ficiency which is in agreement with the predictions of dif- 
fusive shock acceleration theor y ( Blandford and Eichlerl 
Il987l : IJones and Ellison! . Il99lf ). After injection into the 
interstellar medium (ISM), cosmic rays remain contained 
in the gaseous disk for ~15 myr and in the Galaxy 
for ^100 Myr be f ore es caping into intergalactic space 
( Berezinskii et al\ Il990f ). Note that the latter value 
is much larger t han estimates based on the so-called 
leaky- box model ([Yanasak et aZ.L [2Q0lh : see lStrong et all 
(2007) for a full discussion of this point. 

Propagation of cosmic rays in the ISM is usually mod- 
elled as diffusion, where the energetic particles scatter on 
iregularities (fluctuations) of the turbulent Galactic mag- 
netic field (see Sec. III.D). The diffusion equation may 
include stochastic reacceleration in the ISM, convection 
by the Galactic wind, continuous and catastrophic energy 
losses, nuclei fragmentation, radioactive decay, and pro- 
duction of secondary particles and i sotopes (for a recen t 
review of cosmic ray transport, see lStrong et all 120071 ). 
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Isotopes of light elements (Li, Be, B) in cosmic rays are 
almost all secondaries produced in spallations of heavier 
(CNO) nuclei during CR propagation. If the diffusion is 
fast (slow), the secondary nuclei are present, after prop- 
agation, in small (large) amounts; therefore, the relative 
abundances of secondary and primary nuclei can be used 
to determine the propagation parameters 10 . The derived 
propagation parameters (timescale of CR confinement, 
diffusion coefficient, et c.) are model depen dent and can 
vary significantly (e.g.. iPtuskin et a/] . [2006). 

The production spectra of secondary particles are de- 
termined by the kinematics of the collision and depend 
on the ambient spectrum of cosmic rays while their prop- 
agation is governed by the same propagation equation as 
for other cosmic ray species. The production rate of sec- 
ondary positrons slightly depends on the as sumed propa- 
gatio n model and is about (1 — 2)10 42 s _1 (iPorter et all 
2008), i.e. ~5-10% of the Galactic e + annihilation rate. 
A cosmic ray origin of the positrons annihilating at the 
Galactic center can still be reconciled with the produc- 
tion rate if cosmic ray intensities were significantly higher 
in the past (Sec. IV. B. 5 for an analogous situation for 
the Galactic supermassive black hole). 

Heliospheric influence (modulation) changes the spec- 
tra of cosmic ray particles below ~10-20 GeV/nucleon 
as they propagate from the boundaries of the so- 
lar system toward t he o r bits of the inner planets 
(jGleeson and AxfordL Il968k iParkerl Il965[ ). The helio- 
spheric modulation is a combination of effects of con- 
vection by the solar wind, diffusion, adiabatic cooling , 
drifts, and diffusive acceleration (e.g., IPotgieterl 1 19981 ). 
Though the e + fraction inside the heliosphere is small 
(-0.1), the e+ flux in the ISM below -1 GeV is esti- 
mated to be comparable to C R electron flux at the same 
energies ( Strong et a/ll200"3 ). 

Direct information about the CR fluxes and spectra 
in distant locations is provided by the Galactic diffuse 
7-rays. Continuum diffuse emission is expected in the 
hard X-ray and 7-ray regime from the physical processes 
of positron annihilation (through formation of positro- 
nium), inverse Compton scattering and bremsstrahlung 
from CR electrons and positrons, and via decay of neu- 
tral pions produc ed by interactions of CR nuclei w i th the 
interstellar gas ([Porter et all 120081 : IStrong et all . l2000l 
I2004D . Positron annihilation in flight (continuum) may 
contribute in the MeV range (fAharonian and Atovanl 
200(| iBeacom and YiikselL l2006f ) . That contribution can 
be determined from a comparison of model predictions 



10 The stable secondary /primary ratio does not allow one to derive 
a unique set of propagation parameters. The radioactive isotope 
abundances are then used to break the degeneracy. Four radioac- 
tive isotopes, 10 Be, 26 Al, 36 C1, and 54 Mn, are commonly used to 
probe the effective Galactic volume filled with cosmic rays and 
derive the CR confinement time in the Galaxy. Their half-lives 
range from 3.07 10 5 yr ( 36 C1) to 1.60 10 6 yr ( 10 Be) with the 
shortest half-life being most sensitive to the local structure. 



to the data obtaied by INTEGRAL, COMPTEL, and 
EGRET (and now by Fermi/L AT launched in June 2008, 
lAtwood et al. I, l2QQ9h . The analysis of the INTEGRAL 
data shows th at most of the emis sion between 50 keV 
and ~1 MeV (Bouchet et al ., 2008) is produced via in- 
verse Co mpton scattering of b ackground photons off CR 
electrons ( Porter et all l2QQ8h and that CR positrons in 
distant regions of the Galaxy (including the direction of 
the Galactic center) are mostly secondary (see also Sec. 
II.C.l). 



3. Pulsars, milli-second pulsars and magnetars 

A large fraction of high energy sources in the Milky 
Way consists of rapidly rotating magnetized neutron 
stars, which belong to several sub-classes: rotation pow- 
ered pulsars (including Crab-like and Vela-like pulsars), 
accretion powered pulsars (including milli-second pul- 
sars or ms pulsars 11 ) and strongly magnetized rotating 
objects (including soft 7-ray repeaters and magnetars). 
All these objects have been the subject of int ense ob- 
serva t ional and theoretical work (see reviews by iRudak , 
| 200ll: ISeiradakis and Wielebinskil . 120041 : iLorimerl 12005 : 
lHarding and Lail 120061 ; and references therein) . Typical 
values of their main properties (magnetic field intensity, 
rotation period, activity lifetime, estimated birthrate and 
total number in the Milky Way) are provided in Ta- 
ble Enn 

The high energy radiation and/or a high magnetic field 
(B* > 10 12 Gauss) of those objects are associated with 
intense e~-e + pair creation. The pairs are further accel- 
erated in parallel electric fields in the polar caps or in 
the outer gaps close to the light cylinder 12 . The 7 — 7 in- 
teraction of secondary photons produced by the primary 
particles yields a pair cascade; its particles can eventually 
escape into the pulsar wind. Pair creation is accompanied 
by different high-energy photon production channels that 
directly contribute to the pair cascade: curvature radia- 
tion, magnetic inverse Compton scattering, synchrotron 
radiation, and photon splitting. 

It suffices of a few charged particles to be accelerated 
up to high Lorentz factors to initi ate an e~ - e + pa ir cas- 
cade either above the polar caps (Ha rding! . Il98lh or in 
a part of the outer magnetosphere close to the fron- 
tier of the open magnetic fie ld lines region called the 
ou ter gap (ICheng et all\l 9S6). In the outer gap model 



of IZhang and Chengl ( 1997 ). primary e -e + pairs have 



a typical energy E p ~ 5 10 6 P 1 / 3 MeV (where P is 
the pulsar's period in seconds). Photons with energy 
E > E crit ~ 3 B^ 1 MeV (see Sec. IV.B.lb) will generate 



11 ms pulsars spin-up by accreting angular momentum from a com- 
panion star. 

12 The light cylinder is the surface inside which the closed field 
lines have a rotation velocity smaller than the speed of light; it 
separates the region of closed and open magnetic field lines. 
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TABLE VIII Properties of magn etized neutron stars (data for 
pulsars a nd ms pulsars a r e from ILorimerl . 120051 and for magne- 
tars from I Gill and Hevll 120071 although iKeane and Kramerl 
120081 suggest somewhat higher birthrates) . 



Pulsars ms Pulsars Magnetars 



Magn. field 


(B) (G) 


10 i2 


3 x 10 s 


3 x 10 14 


Period 


(P) (s) 


0.5 


3 x 10" 3 


10 


Birthrate 


R (yr- 1 ) 


1.5 x 10" 2 


1(T 5 


2 x 1(T 3 


Lifetime 


(r) (yr) 


10 7 


3 x 10 9 


2 x 10 4 


Total number 


N 


1.5 x 10 5 


3 x 10 4 


40 


e + yield a 


h e ± (s _1 ) 


4 x 10 37 


5 x 10 37 


4 x 10 40 


Total e + yield 6 


N e ± (s- 1 ) 


5 x 10 42 


1.5 x 10 42 


1.6 10 42 



their scaleheight 13 >300 pc and they should have a sub- 
stantial bulge component, albeit with a bulge/disk ratio 
B/D< 0.5. 

The main problem with compact magnetized ob- 
jects as candidate e + sources is the expected high 
energy of the produced posi trons (E~30 MeV, 
Ide Jager and Djannati-Ataj 120081 ), which violates the 
constraint from the continuum MeV emission observed 
in the inner Galaxy. 



4. X-ray binaries and micro-quasars 

X-ray binaries (XRBs) involve a compact object (neu- 
tron star or black hole, hereafter the primary) ac- 
creting matter from a normal star (the secondary) 
through an acretion disk. They are classified as high- 
mass (HMXRBs) or low-mass (LMXRBs) depending on 
whether the mass of the secondary is of spectral type 
earlier or later than B (heavier or lighter than ~4 M , 
respectively). About 300 XRBs are currently detected in 
the Galaxy, for a total X-ray luminosity of ~10 39 erg 
s _1 ; their luminosity function suggests, however, that 
there are more than a thousand of them. LMXRBs are 
~10 times brighter on average and slightly more numer- 
ous than HMXRBs. The corresponding scale heights are 
^410 pc in the form er case and ~150 pc in the latter 
([Grimm et q/1 l2Q02h , 

Some XRBs exhibit radio emission, which is usually 
attributed to synchrotron radiation emitted by leptons 
(electrons, and perhaps positrons); leptons are launched 
along diametrically opposite jets fuelled by the accretion 
energy. If the jets are confirmed by imaging, the system is 
called a microquasar (/iQ), for it has a similar structure 
to quasars (t he latte r bein g a million times larger and 
brighter, i.e. iMirabell 120081 and references therein). 

The physics of microquasars is extremely complex and 
there is no generally accepted model at present. Em- 
pirical evidence suggests that at high accretion rates, 
X-ray emission peaks at ~1 keV (high- soft state), 
whereas at low accretion rates X-ray emission appears 
at higher energies, with a power-law spectrum and an 
exponential cut-off at ^1 keV (low-hard state, e.g. 
iMcClintock and RemillardL 120061 ). It is expected that 
a persistent j e t wil l be present in the low-hard state 
(Fender et al and correlations between the ra- 

dio and X-ray luminos ity and the accretio n activity have 
been proposed (e.g. [Corbel et al. , 2003). It should be 
noticed that the issue of positron production in micro- 
quasars suffers from two major uncertainties: i) while 
electrons are certainly responsible for the observed syn- 
chrotron emission, it is not yet clear whether the posi- 
tively charged component of the jets consists mainly of 



Storv et al. (2007) adopt a scaleheight of 200 pc for ms pulsars, 
but in view of the age of those systems - and the additional effect 
of a kick velocity from the explosion - this should be >300 pc. 



a: Individual source yield from Eq. [9] b: Galactic yield from 
N e ±—h e ±R (r), assuming £=1 (see text). 



a pair cascade involving a total of N e ±=E p / E CT - lt pairs, 
most of which will be reflected by the magnetic mirror 
effect and then move towards the light cylinder. 

In terms of the surface magnetic field and the pulsar 
period the total e + production rate of the cascade is: 

n e± ~ 2.8 10 37 BfH P- 8 / 21 s" 1 (9) 

where £?d,i2 is the dipole magnetic field in 10 12 Gauss. In 
the case of ms pulsars the dipole assumption at the stel- 
lar surface is not valid anymore and the magnetic field 
should be rescaled as — > Bd(R*/£) 3 , where R* is the 
star radius and t ~1 km is the curvature radius of the 
magnetic field in the stella r surface (approxim ately equal 
to the stellar crust radius. [Wang et aD . l2QQ6h . The effec- 
tive rate of positrons injected into the pulsar wind is a 
fraction £ of h e ±. This fraction is probably lower than 
one in normal pulsars and close to one in ms pulsars: 
due to the lower magnetic field of the latter (B ~ 10 8-9 
G) the light cylinder is much closer to the neutron star 
surface and particl es are expected to escape more easily 
([Wang et all I2QQ6I ) . Notice that in the extreme magnetic 
field conditions of magnetars the prod uction of most of 
the e + -e~ pairs is probably suppressed ( Harding and Lail 
l2006h . 

The total Galactic injection rate N e ± of one partic- 
ular class (normal pulsars, ms pulsars, magnetars) is 
N e ±=h e ±N , where n e ± is the average e + production rate 
of one source and N =R (r) is the number of sources in 
the Galaxy (where R is the birthrate and (r) the typical 
lifetime of the sources, see Table [VTlT]) . 

In view of their young age, pulsars and magnetars are 
expected to have a radial distribution closely following 
the one of the star formation rate (Fig.[TQ|). a small scale- 
height (~100 pc) and an insignificant population in the 
bulge. Millisecond pulsars are expected to have a dif- 
ferent radial distribution, since they originate in binary 
systems of all ages. For that reason, their radial distri- 
bution should be closer to the one of SNIa (Fig. [TP]) . 
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ions or positrons; and ii) even if positrons are largely 
present, it is not yet known whether the e~-e + pairs 
of the jets are ejected at ultra-relativistic velocities or 
not: current wisdom is that the jets are mildly relativis- 
tic ( Gallo et all |2QQ3[) but al t ernat ive views have been 
expressed (e.g. iFoellmi et all . |2QQ8[ ). This has obvious 
implications for the inflight e + annihilation and the pro- 
duction of >1 MeV 7-ray continuum. 

Positrons can be pair-created in the vicinity of the com- 
pact object, either in the hot inner accretion disk, in the 
X-ray corona surrounding the disk, or at the base of the 
jets; the latter may channel a fraction of the e~-e + pairs 
out of the system. Alternatively, if the jets consist mainly 
of relatively cold e~-e + or e~-p plasma, they may create 
new e~-e + pairs at t he ter mination shock with the ISM. 
iHeinz and Sunvaevl ([2002) noticed that the total kinetic 
luminosity of microquasar jets in the Galaxy, evaluated 
at Lkin ~ 3 10 38 erg s _1 , can produce up to 4 10 43 e + 
s _1 , i.e. more than required from observations. This es- 
timate requires ~5% of the kinetic power to be converted 
to e~-e + pairs, i.e. a reasonable conversion efficiency. 
Soon after the fi rst data release o f the 511 keV image by 



SPI/INTEGRAL,lPran^Hi (|2Q0j) noticed that i) the ob- 
served distribution of LMXRBs in t he Galaxy is s t rongl y 
peaked towards the central regions (jGrimm et all . E002) , 
similar to that of the 511 keV emission, and ii) their total 
X-ray luminosity is ^10 39 erg s _1 , a hundred times larger 
than the corresponding mass-energy of the observed 10 43 
e + s _1 in the Galaxy. He also noticed, however, that 
most of the strongest sources (accounting for 80% of the 
total Galactic X-ray flux) are evenly distributed in the 
Galactic plane, with no preference for the bulge; he con- 
cluded that, if the positron emissivity scales with their 
X-ray flux, then LMXRBs cannot be the origin of the 
bulge Galactic positrons. The argument is invalid, how- 
ever, if the timescale for the variability of the X-ray flux 
is much smaller than the slowing down timescale of the 
positrons, and/or if positrons annihilate far away from 
their sources. 

Various features of the scenario of microquasars as 
positr on producers were studied in iGuessoum et all 
(200 6]), on the basis of i) existing theoretic al models (e.g. 
iBeloborodovl Il999fc lYamasaki etaJL Il999h and ii) global 
energetic considerations (of XRB luminosities correlated 
to jet power and to positron ejection rates). They found 
that such considerations lead to average values up to 10 41 
e + s _1 for a jet. Interestingly enough, this is roughly the 
current upper limit for e + produ ction rates in XRBs wit h 
SPI/INTEGRAL (see Table 1 in lGuessoum et qZll2006h . 
If a hundred microquasars exist in the Milky Way (not an 
unreasonable extrapolation from their currently known 
population of two dozen), then these objects may con- 
tribute substantially to the observed 511 keV emission. 
The distribution of the known microquasars shows indeed 
some clustering towards the inner Galaxy, but the data 
are insufficient for statistically significant conclusions. 

A similar inve stigation is performed in 
iBandvopadhvav et all (|2QQ9h . who consider hadronic jets 



(also containing e~-e + pairs) launched by all LMXRBs, 
down to the lowest X-ray luminosities. Extrapolating 
from the results of recent deep X-ray surveys of the 
central bulge, they estimate that a bulge population 
of 300-3000 LMXRBs would inject mass in jets at a 
rate of 10 1T — 10 18 g s _1 ; the observed e + production 
(=annihilation) rate of 2 x 10 43 s _1 requires then a yield 
of 40-400 e~-e + pairs per proton. As an example the 
authors discuss observations of the giant galaxy M87, 
and they argue that such a high ratio is justified by the 
observational finding that the jet plasma in that galaxy 
is cold, i.e. the energy spectrum of the electrons has no 
measurable low energy cut-off and it is thus dominated 
by low energy particles (electrons and positrons). They 
also argue that this effect allows one to satisfy the 
observational constraint in the MeV range ( Sec. II. B. 2). 
However, an analysis of the M87 jet by iDunn et 



(2006|) concludes that the jet is e -e + dominated only 



under the assumption of a low-energy cut-off E m i n ~0.5 
MeV, but higher E min (implying smaller pair fractions 
per proton) cannot be excluded and are even suggested 
by polarization measurements. It is clear then that 
neither the energies nor the abundance of positrons in 
the jets are known at present. 

Note that in bo t h IGuessoum et al\ (|2006f ) and 
IBandvopadhvav et all ( 20091 ) the bulge /disk ratio can 
only be reproduced if it is assumed that a fraction of 
the disk positrons do not annihilate in the disk, but 
leave it altogether. These argume nts will be discussed in 
Sec. IV . D along with the claim of Wei denspointner et all 
(|2008af ) that their recent finding of asymmetric 511 
keV emission from the inner disk (Sec. II. B. 3) favours 
LMXRBs in the hard state as the e + sources. 



5. Positron production by the Galactic black hole 

As already discussed in Sec. II. B, the supermassive 
black hole (SMBH) in the Galactic center (GC) had been 
suggested already in the 1980s as a e + source, on the 
basis of the varia bility found in the HEAO-3 data by 
iRiegler et all ( 198lh . Variability is not an issue anymore, 
but the first year data of SPI/INTEGRAL revived the 
idea of the SMBH as e + source because of the difficulty 
met by other candidate sources to explain the 511 keV 
image. 



Compared to the situation in the 1980s, the SMBH 
models have to cope with two new requirements: 

i) The emission does not originate from a point source 
in the GC, but from a region extended over the whole 
bulge. This implies that positrons from the central source 
have to travel distances comparable to the bulge radius, 
i.e. ~2 kpc. 

ii) Sgr A*, the mult i- wavelen gth emission sourc e at the 
GC, is notoriously weak (e.g. Eck art et aZ.L l2QQ8h : its X- 
ray luminosity is ~10 35 erg s _1 and its bolometric lumi- 
nosity is estimated to ~10 36 erg s _1 , i.e. ~3 10 -9 times 
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the corresponding Eddington luminosity 14 . The X-ray 
emissivity of Sgr A* is ~10 4 times weaker than the com- 
bined emissivity of the population of Galactic XRBs (see 
Sec. IV. B. 4) or even some individual XRBs; if positron 
production is correlated to X-ray emissivity, it is difficult 
to conceive Sgr A* as important e + source. 

It turns out that a viable solution may emerge, sat- 
isfying both constraints, if one drops the assumption of 
steady state (i.e. equality between e + production and an- 
nihilation rates): since positrons need time to slow down 
and to fill the bulge, one may invoke a much higher activ- 
ity of Sgr A* in the past. The question is then whether 
that high activity was due to a rare event (in which case 
the current low activity represents the normal state of 
Sgr A*) or whether the past high activity was the norm 
(in which case today's low activity is of low probability). 
Proposed models explore both possibilities (Fig. ITT]) : 

a) Tidal disruption of nearby stars and subsequent ac- 
cretion of their material may boos t the a ctivity of Sgr A* 
for timescales of 10-100 yr ([Reesl 1 19881 ). In view of the 
star density in the vicinity of Sgr A* such events may 
occur e very r caT) ~ 10 4 — 10 5 yr (e.g. ISver and Ulmerl , 
1999h. iFatuzzo et ali 



IF - 10 5 yr (e.g. M , 

] (l200lh . ICheng et ali (I2006D and 



Cheng et al suggest that such events are at the 



origin of past high activity of Sgr A*. 

b) Quasi steady-state accretion of surrounding gas was 
^10 4 times higher in the past, but it was interrupted 
by some external facto r which destr oyed the accretion 
flow ~300 years ago (iTotanil . l2QQ6h . This timescale 
results from the X-ray emission observed from Sgr B 
and Sgr C (both at distances of ~75 pc from Sgr A*, 
see Fig. [T8|) which has been interpreted as the delayed 
reflexion of (a much stron g er) past emission from Sgr 
A* (e.g. iMarkevitch et~ali, 1993. iKovama et oD, [2006, 
iMurakami et all l2000l iReynivtsev et all 120041 ) Build- 
ing on those ideas, ITotanil (|2006f ) suggested that it is the 
expansion of the SN remnant Sgr East which destroyed 
the high accretion flow. The age of that remnant (a few 
10 4 yr) is smaller than the estimated timescale between 
two SN explosions in the GC vicinity (~10 5 yr), but the 
probability that we observed Sgr A* just 300 yr after its 
crossin g by the expa nding shell of Sgr East is rather small 
f~l%. ITotanil . l2QQ6h . 

It should be noted that the timescale for variability 
of 511 keV emission is the longer of the two timescales 
Tprod (for variability of e + production in the above mod- 
els) and r ann (for e + slow down and annihilation). The 
former is ~10 5 yr, while the latter depends strongly on 
the physical conditions of the ISM (Sec. V). However, if 
the positrons survive annihilation in the GC vicinity and 



14 The Eddington luminosity is the limiting value for which gravi- 
tational attraction of a point source of mass M (accreting from 
surrounding material) is matched by repulsive radiation pres- 
sure due to accretion luminosity L. It is given numerically by 
LEdd ~l-3 x 10 38 (M/Mq) erg s _1 (assuming spherical sym- 
metry and Thompson scattering of radiation). 
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FIG. 17 Illustration of the two schemes conceived for the 
positron production activity of Sgr A*. In (1), protons are 
injected/accelerated every ^10 5 yr following the tidal dis- 
ruption of a star, and their collisions with the ISM produce 
positrons (after 7r + decay). In (2), Sgr A* produces quasi- 
continuously positrons at high rate, except when its accre- 
tion flow is interrupted by the passage of the shock front of 
a nearby SN explosion (expected to occur on a timescale of 
^10 5 yr); the last SN explosion created the SN remnant Sgr 
East a few 10 4 yr ago and its expansion interrupted the ac- 
cretion flow 300 yr ago (according to Totani 2006). Despite 
the discontinuous e + injection in both cases (1) and (2), the 
annihilation of positrons and the resulting 511 keV emission 
are in a quasi-steady steady state (3) if positrons diffuse in 
the bulge, because of the long timescale of the latter process 
(~10 7 yr). 



manage to fill the bulge (and this is a big if), the corre- 
sponding diffusion timescale can be estimated by quasi- 
linear diffusion theory as r^iff ~10 7 yr (e.g. Jean et al. 
2006). It is this long timescale that determines variabil- 
ity of 511 keV emission and makes it essentially constant 
(Fig.H7|). despite the variability of the e + production rate 
in both cases (a) and (b). We shall discuss the problem of 
e + propagation in the bulge (a crucial issue for all models 
producing bulge positrons in the SMBH) in Sec. VI. 

Production of positrons can be envisioned by either 
a) collisions of protons accelerated by the SMBH with 
the ISM and subsequent 7r + decay or, b) direct pair- 
production (Sec. IV.B.l). We briefly present below a 
few models which constitute specific realisations of the 
aforementioned ideas. 

a. Production of positrons via p-p collisions 

Models of this class usually try to reproduce simul- 
taneously observations of 511 keV emission from the 
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FIG. 18 Schematic representation of the inner Galaxy at different linear (bottom and left axis) and angular (top and right axis) 
scales. The main features discussed in Sec. IV. B. 4 (SMBH scenario) appear on the figures, as e.g. CND (Circumnuclear disk, 
left), CMZ (Central Molecular Zone, middle) and the tilted HGD (Holed Gaseous Disk, containing atomic and molecular gas, 
right). If indeed Sgr A* is the main positron source, (most of) its positrons have to diffuse through the CMZ into the bulge. 
Notice that the size of the CMZ corresponds to the size (FWHM) of the inner bulge in the first model of SPI data fitting (Table 
I). The size of the outer SPI bulge in that same model is indicated by a dotted circle in the rightmost panel. 



bulge 



of higher energy 7-rays from the central re- 



protons every ~10 7 yr) it is impossible to explain simul- 
taneously the 511 keV and high energy 7-ray emissions: 
although positrons take a long time to diffuse and annihi- 
late, protons interact rapidly and the corresponding 7-ray 
emission fades out in less than 1Q 5 yr. I n a subsequent 



version of that work, Ch eng et al\ 



gions. Such emission has been detected at >500 
MeV with EGRET/CGRO (|Maver-Hasselwander efalL 
Il998h and at TeV energies with various ins truments 
(Aharonian et al ], 120041 : lAharonian elo/l [2006, and ref- 
erences therein). 

The generally large timescale for e+ thermaliza- the idea proposed by jFatuzzo et al 
tion/annihilation (see Sec. V) creates problems to any 
scheme trying to explain simultaneously current obser- 
vations of 511 keV and higher energy 7-rays, if acceler- 
ated protons are assumed to be at t he origin o f bot h 
emissions. This was realized in iFatuzzo etaR (|200lh . 
who explored the fate of positrons produced by energetic 
protons, required to explain the GeV emission detected 
by EGRET from the source 2EG J1746-2852; the latter 
is located in the inner arcmin of the Galaxy and its 7- 
ray emiss ivity is asso c iated to proton acceleration i n Sgr 
East by iMelia et all (|l998f ). IFatuzzo et all (|2Q0lh find 
that in the physical conditions of Sgr East, the e + ther- 
malization timescale is much longer than the age of Sgr 
East (or that the current e + thermalization/annihilation 
rate is much lower than the e + production rate, which 
is inferred from the high-energy 7-ray emissivity of 2EG 
J1746-2852). They conclude that recently accelerated 
protons in Sgr A East cannot be the source of the ob- 
served annihilation radiation, unless some more efficient 
e + cooling mechanism is at work. Alternatively, they 
suggest that more positrons from previous episodes of 
activity in the Galactic center (i.e. from tidal disruption 
of stars) have been "stockpiled" in that region and are 
annihilating now. 

ICheng et al\ ( 2QQ6h explore further those ideas, but 
they consider in more detail the propagation of protons 
and positrons in the Galactic bulge. They find that in 
the case of an energetic but rare event (such as the dis- 
ruption of a 50 M star, releasing 10 54 erg in energetic 



200 71) return then to 
(|2QQlh . namely that 
the observed 511 keV emission results from the annihi- 
lation of positrons produced and stockpiled in the bulge 
by dozens of tidally disrupted low mass stars in the past 
^10 7 years: each one of those events is less energetic 
(-10 52 erg) but more frequent (corresponding timescale: 
~10 5 yr) than the rare massive star disruption. Because 
of the short timescale of e + production with respect to 
the timescale of e + annihilation, the resulting 511 keV 
emission is quasi-steady (see Fig.[T7j). Thus, although the 
invoked e + production event is of low probability (small 
frequency), the observed intensity of the 511 keV emis- 
sion is not. 

Obviously, most of the positrons are pro duced in the 
high- d ensity region where protons interact. ICheng et all 
(2007) argue that positrons do not annihilate in those 
regions, but they avoid them (because of the screening 
effect produced by magnetohydrodynamic waves excited 
near molecular clouds) and they propagate through the 
intercloud medium. In general, models where positrons 
result from the production and decay of 7r + will have 
typical energies of ~30 MeV, resulting in too much 
emission at >MeV ener gies from in-flight annihilation. 
IChernvshov et all (|2QQ9h find that this difficulty may be 
circumvented if the magnetic field in the galactic bulge is 
high enough (>0.4 mG), because in that case positrons 
lose their energy rapidly (before significant in-flight anni- 
hilation) through synchrotron emission; this cooling pro- 
cess takes place within the timescale of e + production, i.e. 
long before e + annihilation, and should be currently un- 
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detectable in radio. However, as discussed in Sec. IIII.D) 
observations do not favour, a present, such high values 
for the magnetic field of the inner Galaxy. 



b. 



-e pair production by photons 



As in the case of XRBs, pair production around the 
Galactic supermassive black hole may occur either in the 
inner hot accretion disk, in a corona above it or in the 
jets; the latter case has not been considered up to now, 
in the absence of relevant observational evidence in Sgr 
A*. 

iBeloborodovl (Il999h studied a detailed model of pair 
production, resulting from collisions of 7-rays from the 
hot inner disk with X-rays from the outer disk; the result- 



data suggests a narrow bulge component (Table I), the 
size of which (3° FWHM) corresponds to the size of the 
CMZ: ~l/3 of the bulge positrons may annihilate there 
and the remaining 2/3 may diffuse in the outer bulge 
(11°FWHM in Table I, see also Fig. [18]). This picture 
may have difficulties, however, with the results of the 
spectroscopic analysis of SPI (Sec. V.E) which finds 
that bulge positrons annihilate mostly in a warm medium 
(neutral or ionized), not in a molecular one. 



C. Dark matter and "non-standard" models 



+ pairs are blown outwards from the radiation L General P ro P erties of dark matter P artides 



mg e~-e 

pressure of the disk, in a mildly relativistic wind. The 
model is not specifically designed for the case of Sgr A* 
but rather for extragalactic black holes in Active Galac- 
tic Nuclei (AGN), but its results can be extrapolated to 
the case of Sgr A*. The maximum pair production is 
obtained for a dense, optically thick, wind. In that case 
positrons annihilate mostly inside the wind flow and pro- 
duce a very broad annihilation line, unlike the one ob- 
served with SPI/INTEGRAL. Positrons can escape and 
annihilate in the ISM in the case of a less dense, opti- 
cally thin, wind but in that case their production rate is 
substantially smaller. 

Totani (2006) considered pair production in the hot, 
inner accretion disk during past phases of higher ac- 
tivity in Sgr A*. The invoked past accretion rate, 
rh ~lO _4 M /yr, is not extravagant and could easily re- 
sult from the material released by the regular tidal dis- 
ruption of nearby low mass stars (as in Gheng et al . 2007 ) 
or from winds of nearby massive stars ( Quataertl [2004) , 
although Totani (2006) assumes that it originates from 
the ionized "halo" surrounding Sgr A* (see Fig. [T8j left). 
Totani (2006) considers pair production in the frame- 
work of the so-called Radiatively Inefficient Advection 
Flow (RIAF) model for accretion disks; this model, de- 
coupling accretion from emerging luminosity, has been 
appli ed with con s iderable success to the case of Sgr A* 
(e.g. [Yuan et all I2QQ4J : IXu et all . 120061 ). Totani (2006) 
finds that the very high temperatures of the inner disk 
(T^IO 11 K) implied by the RIAF models are essential 
for a high rate of e + production, which he evaluates up 
to ~10 43 e + s _1 , i.e. close to the observationally inferred 
annihilation rate. 

In all models of e + production from Sgr A*, e + anni- 
hilation occurs on much longer timescales than e + pro- 
duction and varies much less in time than the latter. An 
advantage of direct pair production models with respect 
to those involving energetic proton collisions is the low 
energy of the positrons produced, allowing them to sat- 
isfy the constraint of the observed MeV continuum. 

In all cases, it has to be demonstrated that positrons 
may diffuse from Sgr A* throughout the bulge without 
excessive annihilation in the dense inner regions (the Cir- 
cumnuclear Disk), which would give a strong, point-like 
emission. On the other hand, the latest analysis of SPI 



In the past 40 years or so, particle physicists searched 
for possible dark matter (DM) candidates meeting 
three basic requirements, namely: stability (at least on 
timescales comparable to the age of the Universe), charge 
neutrality (to avoid electromagnetic - ELM - interactions 
and prevent DM to shine) and with a non negligible mass 
(so that it can contribute gravitationally) 15 . The absence 
of electric charge favours DM models with weak interac- 
tions. However weakly interacting particles may suffer 
from a prohibitive "free-streaming" effect 16 , depending 
on their mass: for example, if DM is composed of mass- 
less neutrinos, the formation of Milky Way size galaxies 
is strongly suppressed. Hence there is a lower limit on the 
mass of weakly interacting DM candidates to explain the 
formation of the smallest ojects observed in the Universe, 
which is of ~a few keV. This leads to the notion of weakly 
interacting massive particles (WIMPs) and the idea of 
collisionless dark matter. The existence of WIMPs with 
a mass in the GeV-TeV range is compatible with the 
absence of signal in present DM direct detection experi- 
ments. However, the fact that none of thes e expe riments 
has observed a positive signal yet (iBaudisl l2QQ7h 17 may 
lead to different interpretations, as discussed below. 

If the DM number density today n^ m were similar 
to the relativistic particle density (n 7 ~ 400 cm -3 for 
photons), the DM mass-energy density pdm — ^dmX^dm 
would exceed the critical density p c by several orders of 
magnitude, for DM particle masses m^ m >1 keV. Given 
that particles lighter than keV are forbidden by the free- 
streaming argument, DM particles which have been in 
thermal equilibrium with radiation at some stage of the 



15 These criteria are in fact supported by recent observations (as- 
suming a Friedman- Robertson- Walker Universe). For example, 
extrapolation of the physics of ordinary baryons to DM suggests 
that DM ELM interactions should damp the DM primordial fluc- 
tuations on a cosmological scale, and prevent the formation of 
small scale structures (smaller than a Milky Way size galaxy). 

16 Free streaming refers to the motion of non-interacting particles 
endowed with some initial velocity across the Universe. It has 
the effect of erasing irregularities on scales smaller than the free- 
streaming length. 

17 Apart, perhaps, from DAM A/LIBRA, (|Bernabei et al.l \2003). 
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cosmic evolution should subsequently disappear. This 
may occur in two ways: either through an extremely 
small decay rate, which ensures a lifetime comparable to 
the age of the Universe (and thus implies quasi-stability) 
or through annihilation processes. 

The latter mechanism has received a lot of attention in 
the last three decades. The requirement pdm °° Pc implies 
that the annihilation cross-section of DM particles should 
be c omparable t o the weak interact i on cross section, see 
e.g. iHutl ([l977l ). |Prvor et all (|l980f ). I Abbott and Sikhdel 
(|l983h 18 . In some cases, this requirement can also be 
used to constrain the DM ma ss. A s was pointed out by 
iLee and Weinberg! (|l977h and <|l977h , the annihila- 
tion rate (< av >oc v 2 ) of fermionic particles with typical 
weak interactions has a square dependence on DM mass 
rn 2 drn . The observed pdm implies then that DM particles 
should be heavier than a few GeV. This constitutes an 
extra motivation for considering WIMPs as DM candi- 
dates. 

About twenty five years ago, such properties (weak 
but non negligible interactions and rridm >1 GeV) sug- 
gested that direct detection of DM would be relatively 
easy. However: 

1) The DM spin independent interactions with matter 
(as measured by direct detection experiments) is at least 
eight orders of magnitude weaker than t he weak interac- 
tions ([Angle et a/.U2QQ7l : iLemranil . 120061 ) 

2) Indirect detection experiments find no "smoking gun 
"evidence (i.e. the emission of a monochromatic line at 
an energ y E = rridm) allowing for a c l ear id entification 
of rridm ((Abdo et all 120091 : lAharonianL l2QQ9h . 

3) No signature of new physics, which would indirectly 
validate the exist ence of DM particles, was found at LEP 
or TEVATRON (lAbbiendi et all Il999b iPiperl 120091) . 

The aforementioned facts could either mean that DM 
has much weaker interactions than the Standard Model 
or point towards very heavy or very light DM particles 19 . 
Besides, the absence of ELM interactions could also im- 
ply strong interactions but at differe nt energy scale than 
previously considered ( Boehml l2QQ8h . 

The aforementionned constraints, individually taken, 
can be easily accommodated in existing DM models (e.g. 



18 The abundance of DM particles today fix is fixed in the so- 
called "freeze-out" epoch, when the expansion rate of the Uni- 
verse H (a function of fix) equals the DM annihilation rate 
r = nx < av >, where nx is the DM particle abundance and 
< av > their annihilation cross-section; this leads to a "relic" 
abundance of DM particles of fix = /(< &v >) and for the ob- 
served fix ~0.25 one obtains < av >~ 10 -26 cm 3 s _1 , which is 
close to the value of weak interaction cross sections (the actual 
value depen ds on the nature of the DM candidate and its mass, 
see e.g. iBertone et all I2004T ). 

19 If DM is lighter than a few GeV, its interaction with matter 
would be essentially "invisible" for current detectors (MeV par- 
ticles would require for example detectors with eV energy thresh- 
old, while they are currently in the keV range). If DM particles 
are too heavy, their number density in our Galaxy is too small 
to generate a significant number of events in a detector. 



WIMPs). However, when combined together, they ac- 
tually eliminate many proposed models. For example, 
to reduce the tension between direct detection experi- 
ments, which are now sensitive to elastic cross sections 
of the order of 10 _43 cm 2 (corresponding to < av >~ 
10 -33 cm 3 s _1 ), and the annihilation rate av ~ 3 10 -26 
cm 3 s _1 imposed by the "relic density" criterion (corre- 
sponding to a ~ 10 -36 cm 2 ), o ne often ha s to "d ecouple" 
the corresponding processes ([Ellis et al I l2QQ0h 2Q . The 
lack of evidence could also mean that rridm <<1 GeV or 
m^m >>1 TeV. In the case of sub-GeV DM, one may 
naturally circumvent the Hut-Lee- Weinberg limit if DM 
has been produced out of thermal equilibrium or if it is 
made of scalar particles with non chiral couplings to Stan- 
dard Model particles. Non-thermal production also helps 
very heavy DM particles to avoid conflict with the relic 
density criterion. Thus, the acceptable range for rridm 
still lies from the sub-eV range (with axions) to several 
TeV (such as "excited" DM, Kaluza-Klein particles) and, 
in fact in some case, to even higher masses. Hence, keV, 
MeV and sub TeV candidates (such as, respectively, ster- 
ile neutrinos, light dark mat ter, and neutra linos) remain 
potential solutions (see, e.g. iBertond , 120071 ). 

Progress in the field of DM may come from indi- 
rect detection. Indeed, cosmic ray spectra "anoma- 
lies" (with respect to "standar d" astrophysical pre dic- 
tions) appear puzzling enough ([Picozza et all l2QQ7h to 
open u p new particl e phys ics scenarios ( Adriani et all 
l2009l l?l: IChang et all l2QQ8h - For example, recent re- 
sults from the PAMELA satellite indicate an excess of 
CR positrons above 10 GeV over the background ex- 
pected from CR intera c tions with interstellar matter 
([Moskalenko and Strong) . Il998l h but no corresponding 
excess in antiprotons 21 . However, the all lepton spectrum 
as measured by t he Fermi- L AT cam e out flatter than pre- 
viously thought ( Abdo et a/.L 12009). The sources of these 
particles (or electrons and positrons) should be less than 
~1 kpc away, since electrons and positrons in this energy 
range suffer heavy energy losses through inverse Comp- 
ton and synchrotron processes. Conventional astrophys- 
ical sources (e.g. nearby pulsars, see Sec. IV. B. 3) could 
explain these excesses, but the possibility of the first indi- 
rect detection of DM annihilation created a lot of excite- 
ment in the particle physics community (see lEssig et all 
120091 and references therein). 

The recent excitement illustrates the new trends in the 
DM particle physics community and shows how far we 



20 This decoupling can be made by invoking either "co- 
annihilation", i.e. annihilation of DM with another particle, 
present during the dark matter transition to the non-relativistc 
regime (Griest and Seckel, 1991), or by fine tuning the DM pa- 
rameters so that the annihilation cross section is enhanced but 
the elastic scattering cross section remains very small. 

21 The balloon-borne experiment ATIC reported an excess of elec- 
trons plus positrons in the 300-8 00 GeV range, but the excess was 
not confirmed by Fermi-LAT (jAbdo et all l2009t ) and by HESS 
(|Aharonianl 120091 ) 
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are from the determination of the nature of DM. This 
also demonstrates that considering only neutralinos or 
Kaluza-Klein particles as DM may be too restrictive. In 
this review, we shall focus only on the DM candidates 
which have been explicitly invoked to explain the Galac- 
tic 511 keV emission. 



2. Specific dark matter candidates for e + production 

Positrons produced by DM annihilation or decay may 
annihilate in flight, before losing a large fraction of their 
energy. Inflight e + annihilation would provide an addi- 
tional source of continuum 7-ray emission towards the 
Galactic centre (Sec II.C.l) which can be used to con- 
strain the e + energy at injection and, consequently, the 
mass of annihilating or decaying DM particles. Assum- 
ing that the contribution of known astrophysical sources 
to the measu red continuum in this ene rgy r ange is well 
under stood, iBeacom and Yiiksell (|2006f ) and ISizun et all 
(2006) obtained a mass upper limit of m^ m ~ a few MeV 
(see Fig. [6]). Among the many proposed DM scenarii, 
those that may satisfy such a constraint can be classified 
in two main categories: 

i) Light DM particles of ^MeV mass, either annihilat- 
ing (IBoehm et all l2004bt [Gunion et al , 20061). or decay- 
ing ([Hooper and Wand, l2004t iPicciotto and Pospelovl 
l2005l ) or even both ([Pospelov et a/ll2008l ). 

ii) Heavy DM particles in the ~GeV - TeV range, 
de- exciting (or decaying into another particle) with a 
mass difference of a few MeV between initial and final 
state s ( Finkbeiner and Weinerl [20071 : [Pospelov and Ritzl 
l2QQ7h . 

Other, more intricate possibilities, involving dark mat- 
ter, cosmic strings, primordial black holes and other ex- 
otica, will be briefly presented in Sec. IIV.C.3I 

Low mass annihilating DM particles were initially pro- 
posed to illustrate a new damping effect 22 but also 
as a counter exampl e of the " Hut-Le e- Weinberg limit" : 
Boeh m et al\ (|2004ah pointed out and IBoehm and Favetl 
(2004) showed that this limit is only valid in the case 
of fermionic DM candidates interacting with Fermi (i.e. 
weak) interactions. But, if one or both of those assump- 
tions are relaxed, a very different conclusion may be ob- 
tained: for instance, if the DM particle is a scalar (spin-0) 



matter relic density in that case. Such particles are ex- 
pected to annihilate into e~-e + pairs (Fig. [19]) either via 
a heavy charged parti c le exch a nge or a new neutral gaug e 
boson (IBoehm et all l2004al : IBoe hm and Favetl, I20C 



which annihilates into a e~-e + pair via the exchange of a 
fermionic particle F, then the relic density criterion con- 
strains the characteristics (mass and couplings) of the F 
particle instead of the mass of the DM particle. Parti- 
cles substantially lighter than a few GeV (and down to 
the MeV range) may then account for the observed dark 



IBoehm and Ascasibarl (|2004l ) and lBoehm and Silkl (j2C 
found that the properties of such particles, if they are at 
the origin of the 511 keV line, should affect the value of 
the fi ne structure constant. In fact , usin g this very argu- 
ment. IBoehm and Ascasibarl (2004) and lBoehm and Silkl 
(|2QQ8h could exclude DM particles heavier than 7 MeV 
(assuming a NFW dark matter halo and the correspond- 
ing best fit cross section). 

The light DM particle idea, at least in it simplest form, 
has been challenged by an analysis of t he explosion of 
supernova SN1987A (iFavet et all l2006h . which puts a 
lower limit of ^10 MeV to the mass of the particle. Such 
a limit is very close to the upper limit allowed by the 
observed ~MeV continuum. 

Low mass (<100 MeV) decaying DM candidates were 
proposed by several groups, after the release of the 
511 keV map by the SPI collaboration, in order to ex- 
plain the large amount of low energy positrons in the 
Galactic bulge. Decay into e + -e~ pairs would be one 
of the dominant decay modes of such particles since, 
apart from the neutrino and photon channels, electrons 
would be the only other kinematically accessible chan- 
nel. Depe nding on the mode l , thes e particles may (the 
axinos in iHooper and Wang . 2004) or may not (ster- 
ile neutrinos in Picciotto and Pospelov . 20051 moduli in 
iKasuva and Kawas aki. 2006| or in ICraig and R^bvll2009h 
be the major contributors to cosmic DM density. 

The second category of DM particle candidates in- 
vokes electroweak scale WIMPs, with masses in the 
100 GeV - 1 TeV range and possessing almost mass 
degenerate excited states, i.e. the difference between 
the excite d and ground states sh o uld b e of the order 
of ~ MeV. iFinkbeiner and Weinerl ([20071 ) noticed that 
the velocity dispersion of DM particles in the gravita- 
tional potential of the inner Galaxy is of the order of 
a few 100 km s" 1 , endowing a 500 GeV WIMP with 
a kinetic energy >511 keV. Inelastic scattering between 
WIMPs could raise one or both of them in their ex- 
cited state (s) and de-excitation to the ground state could 
proceed via emission of a e~-e + pair. This scenario 
has also been invoked to explain other observables such 
as the DAM A /LIBRA signal or the WMAP "haze" 23 
(jHooper et all |2007|, |Fmkbeiner et all l2008h and th e 
PAMELA e + excess (jArkani-Hamed and Weinerl . [2QQ8h . 
A recent investigation of this idea suggests, however, that 
the model parameters have to be pushed to their extreme 
values for the e + production rate to agree with observa- 



The mixed damping effect is in fact analogous to the "Silk 
damping", with dark matter playing the role of baryons 
and neutrinos replacing the photons (|Boehm et al .], 120011 : 
IBoehm and Schaeffer, 2005|) . 



Observations with the Wilkinson Microwave Anistropy Probe 
(WMAP) revealed an excess of microwave emission in the inner 
20° of the Galaxy, which cannot be accounted for by conventional 
astrophysical explanations, such as thermal Bremsstrahlung from 
hot gas, synchrotron emission, etc. 
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FIG. 19 Possible Feynman diagrams for light annihilating 
dark matter particles. 



tions riChen et all . l2QQ9h . 

Independently of their "naturalness" (or lack of) as 
extensions of the standard model, the various proposed 
scenarios for DM particles at the origin of Galactic 511 
keV line also differ as to the predicted spatial profile of 
the resulting e + population. Assuming that positrons an- 
nihilate close to their production sites, this may constrain 
and discriminate between the various models. 

The rate of positrons produced locally by 
annihilation/decay /de-excitation of DM particles is 

given by: h e +(r) oc T x > e ~e+ n x(r), where nj(r) is 

the number density of DM particles at distance r from 

the Galactic center 24 . The interaction rate T x > e ~e+ m 

the case of decaying DM particles has a maximum value 
corresponding to the inverse of the age of the Universe. 
In the case of annihilating or collisionally excited DM 

particles, T x > e ~e+ = < crv > nj, where < crv > is 

the annihilation or excitation cross-section folded with 
the velocity distribution of the DM particles. 

In the case of decaying DM particles, h(r) e + oc n(r)x-, 
i.e. the positron production profile follows closely the 
DM density profile; notice that this also the case for all 
"conventional" astrophysical sources studied in the pre- 
vious sections. On the other hand, in the case of anni- 
hilating or de-excited DM, n e +(r) oc< crv > n 2 v (r), i.e. 
the positron profile is generically more centrally concen- 
trated (because of the n x (r) term), but this shape can 
be also modulated by the possible velocity dependence of 
the interaction cross-section cr(v), since the typical parti- 
cle velocities depend on the gravitational potential 3>(r) 
of the DM halo. In the latter case, comparison of the 
model to the data requires a more elaborate analysis; an 
example of such an analysis will be given in Sec. IV. D. 3. 



3. Other "exotica" 

Several rather "exotic" objects have been invoked as 
sources of the Galactic positrons producing the 511 keV 
emisssion. Here we provide a non-exhaustive liste of 
them. 

For example, iHuh et cd] ([2007) suggested MeV milli- 
charged (fermionic) particles. As in the first version of 



MeV annihilating dark matter, this scenario assumes a 
new light (U(l)) boson exchange, called exphoton. Due 
to the presence of kinetic terms, the dark matter would be 
millicharged, while the rest of the scen ario ressembles the 
fermi onic candidates introduced in ( Boehm and Favetl 
120041 )). The required parameter range to explain the 511 
keV line appears to be compatible with the constraints 
from the relic density requirement and collider experi- 
ments. 

The idea of Q-balls in gauge mediated super- 
symmetry breaking scenar ios was proposed in 
iKasuva and Takahashil ([2005). Q balls can be de- 
picted as stable localized field configurations, their 
stability being guaranteed by a conserved charge Q 
associated to a U(l) sy mmetry. For example, Q coul d 
be the electric charge; in IKasuva and Takahashil (2005), 
Q is in fact the lepton number. These objects may have 
a long enough lifetime and yet a small energy density, 
possibly enabling them to be present in our Galaxy and 
to explain the 511 keV emission. 

Macroscopic objects, like supe rconducting dark mat- 
ter (lOaknin and Z hitnitskvl 120051 ) or compact composite 
objects (Zhitnitsky, 2006) were also proposed. These ob- 
jects are hypohesized to form during the QCD phase tran- 
sition and could be schematically depicted as " quark" - 
balls. They would introduce a link between the dark mat- 
ter and baryonic energy densities nowadays and eventu- 
ally explain why these two quantities are of the same or- 
der of magnitude. Positronium formation would happen 
through electrons or baryonic matter interactions (anni- 
hilations) inside the compact composite objects. How- 
ever, the ability of such a scenario to explain the 511 
keV e mission has been criticized in ICumberbatch et~cd\ 
(2008), who find that positronium formation is hardly 



24 Generally, nx (r) is assumed to be spherically symmetric; devia- 
tions from spherical symmetry, due to triaxiality of the DM halo, 
are negligible with respect to other uncertainties of the problem. 



possible at all in such objects. 

If a tangle of light superconducting strings exist in the 
Milky Way, it may also act as a low en ergy positron 
source (Ferrer and Vachaspatil 120051 [2008) . If the string 
curvature radius is smaller than a characteristic scale, the 
string can move with respect to the magnetized plasma 
in the galaxy. The possible interaction of the string with 
the magnetic field can then generate a current composed 
of zero modes of charged particles. Owing to this mech- 
anism, positrons could propagate along the string. They 
eventually leave the string if they become non relativis- 
ts and scatter with counter-propagating particles in the 
string. Depending on the string curvature radius, the 
positrons would then have an energy of E < MeV, al- 
though this energy could reach the GeV scale if there are 
superconducting strings at the TeV scale. 

Finally, Titarchuk & Chardonnet (2006) proposed that 
X-rays from the SMBH collide with 10-MeV 7-rays from 
small- mass black holes (10 17 g) to give pairs; this can 
produce about 10 42 e + /s, about an order of magnitude 
less than needed. This rate is obtained by taking the to- 
tal X-ray and 7-ray luminosity of the inne r GC regions : 
L x ~ 2 10 39 ergs/s and L 1 ~ 4 10 38 ergs/s ([Strong et all . 
l2000f h then assuming the gamma radiation comes from 
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an optically thick medium and that its spectrum is there- 
fore a blackbody one with a temperature of T 7 = 10 MeV; 
the two energy distributions (X and 7) are convolved to 
compute the pair production rate. A simpler scenario, in- 
vo lving evapora t ing; pr imordial black holes was proposed 
bv lBambi etal\ (|2QQ8h . 



D. Assessment of sources 

In this subsection we summarize the pros and cons of 
each one of the candidate positron sources presented so 
far, in the light of the observational constraints of Sec. 
II. B, namely: i) the total e + annihilation rate ( > 2 10 43 
s _1 ), ii) the typical energy of the injected positrons, or 
the equivalent mass of annihilation DM particles (<3- 
7 MeV) and (perhaps, most significantly) iii) the mor- 
phology of the 511 keV emission (with a bulge/disk ratio 
B/D>1 in the case of a thin disk emission). A fourth 
constraint, namely the longitudinally asymmetric disk 
emission, should be added to this list, once robustly es- 
tablished by further data and analysis. 



1. Positron annihilation rate 

Assuming a steady state regime, the e + annihilation 
rate has to be equal to the average e + production rate 
during the lifetime of e + in the ISM. 

The only source definitely known to provide substan- 
tial amounts of e + at a well constrained rate is the ra- 

+ g-1 xhe decay of 



dioactive decay of 26 Al: 0.4 1 43 e T s 
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Ti probably provides another 0.3 10 e + s . GCRs 
probably provide 0.1 10 43 e + s _1 . Nova models (as con- 
strained against several observables such as eject a abun- 
dances, velocities etc.) may provide a e + yield from 
22 Na decay not be much below the reported value of 
10 41 e + s _1 . The e + of all other candidate sources is 
entirely speculative at present. Values discussed in pre- 
vious sections should be considered as optimistic rather 
than typical values. O bserved upper limits of individual 
sou rces (see Table 4 inlKnodlseder et al. 1. 120051 and Table 
1 in lGuessoum et a/"ll2006h are of little help to constrain 
positron sources. No useful observational constraints ex- 
ist up to now on the e + yields of hypernovae/GRBs, pul- 
sars, ms pulsars, magnetars, LMXRBs, microquasars, the 
SMBH at the Galactic center or dark matter annihilation. 
SNIa remain an intriguing, but serious candidate, with a 
potential Galactic yield of 2 10 43 e + s _1 . 



2. Positron energy 
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FIG. 20 Maps of the Galactic 511 keV emission (flux in cm -2 
s _1 sterad -1 ), as observed from SPI (in all panels, thin iso- 
contours from Wei denspointner et all l2008ah and from obser- 
vationally based or theoretical es timates. A) Observed 26 Al 
(and, presumaby, 44 Ti) map ffrom lPluschke et al. , 2001 ) ; B) 
Accreting binary systems (SNIa and, presumably . LMXRBs , 
see t ext); C) Observed Hard LMXRBs (from iBird et all 
2007). The robustly expected e + annihilation from radioac- 
tivity in the disk (upper panel) is not yet fully seen by SPI. 



croquasars and the SMBH at the Galactic center. Con- 
versely, pair creation involving very high energy photons, 
as in e.g. pulsars or magnetars, will produce positrons 
of too high energy. The same holds for energetic p-p col- 
lisions in Galactic cosmic rays or in the baryonic jets of 
LMXRBs, microquasars and the Galactic SMBH. Those 
processes produce e + of >30 MeV, thus may be discarded 
as major e + sources in the Milky Way. Also, that same 
constraint limits the mass of putative decaying or an- 
nihilating DM particles to <10 MeV, while it does not 
constrain the mass of de-exciting DM particles. 



Radioactive decay produces positrons of E<1 MeV, 
naturally fulfilling the observational constraint on con- 
tinuum 7 rays from in flight annihilation. The same ap- 
plies to pair creation through 7 — 7 collisions in the inner 
accretion disk or at the base of the jets of LMXRBs, mi- 



3. Morphology 

None of the e + sources studied in this section repro- 
duces the large B/D>1 ratio inferred from SPI data, as 
can be seen in Figs. [20] and [2TJ where we present flux 
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FIG. 21 Intensity of 511 keV emission as a function of Galac- 
tic longitude. All fluxes are integrated for latitudes \b\ <15°. 
In all panels, the thick solid curve corresponds to SPI obser- 
vations, i.e. the map of Fig. [20] (Note: We emph asize that 
SPI maps and fluxes are provided here for illustration pur- 
poses only; quantitative comparison of model predictions to 
data should only be made through convolution with SPI re- 
sponse matrix.). The thick dotted histogram (top and mid- 
dle) is the observed longitude distribution of LMXRBs (from 
Grimm et al. 2002); the latter resembles closely the theo- 
retically estimated longitude distribution of SNIa (thin solid 
curve in the upper panel), which has been normalised to a 
total emissivity of 1.6 10 43 e + s~\ with Bulge/Disk=0.45 
(maximum Bulge/Disk ratio for SNIa from Table IV|) . Also, 
in the upper panel, the lower dashed curve corresponds to the 
expected contribution of the 26 Al and 44 Ti /3 + -decay from 
massive stars. The thin solid histogram in the middle panel is 
the observ e d lon gitude distribution of Hard LMXRBs (from 
Bird et al 1 1200/t ) and it has the same normalization as the 
thick histogram. In the bottom panel, the SPI 511 keV profile 
is compared to profiles expected from dark matter annihila- 
tion (after Table IVTl) . 



skymaps and l ongitude profiles, respe c tively, comparing 
the SPI data of lWeidenspointner et all ( 2008al ) to various 
expected source profiles, either theoretical or observed 
ones. The comparison is made under the explicit as- 
sumption that positrons produced from the correspond- 
ing sources annihilate close to them. 



The best-established e + sources, /3 + -decay from 26 Al 
and 44 Ti produced in massive stars, yield B/D<0.2, as 
derived from the observed distribution of the 1.8 MeV 
line (normalised here to a total disk emissivity of 0.7 
10 43 e+ s- 1 , Sec. IV.D.l) (Figs. [20]and[2lJ top). Notice 
that similar profiles are expected for pulsars, magnet ars, 
hypernovae and gamma-ray bursts (albeit with different 
normalizations). 

Binaries involving low mass stars, such as SNIa, novae 
and LMXRBs, are expected to have a steeper longitude 
profile, with a maximal B/D<0.5 (assuming the bulge 
and disk masses of Sec. III). Using data from Fig. [10] 
and Table [V] and adopting an exponential density profile 
for the bulge (with scalelength of 400 pc and normalised 
to 1.4 10 10 Mq) one may estimate an expected sky distri- 
bution and corresponding longitude profile of SNIa, also 
displayed in Figs. l20l (middle) and I2T1 (top) . where it is 
assumed that the e + escape fraction from SNIa is 3%. 
Fig. [21] (top) indicates that the theoretically expected 
SNIa longitude profile ressembles th e observed profile of 
LMXRBs (from lGrimm et all |2Q()3 ) . This similarity re- 
flects the fact that both classes of sources involve an old 
stellar population, which is proportionally more abun- 
dant in the inner Galaxy and the bulge than in the rest 
of the disk. Novae are also expected, on those grounds, 
to have a similar distribution (albeit with a much lower 
normalization constant). 

The upper panel of Fig. [21] clearly shows that 

a) The expected outer disk (I >20°) contribution of 
massive star radioactivity ( 26 A1 and 44 Ti) is not yet de- 
tected, due to insufficient sensitivity; SPI/INTEGRAL 
or a future instrument should reveal that component, or 
else it should be concluded that positrons diffuse far away 
from their sources. 

b) SNIa or LMXRBs or microquasars can explain only 
about half of the strong 511 keV emission from the inner 
Galaxy, assuming they produce as many positrons as in 
Table El) 

c) Any remaining annihilation 7-ray emission requires 
a supplementary source (dark matter or the central 
SMBH, provided its positrons can diffuse to kpc dis- 
tances and fill the bulge); alternatively, it may be as- 
sumed that SNIa or LMXRBs or microquasars produce 
twice as many positrons as assumed in Fig. [211 but half 
of the disk positrons are transported to annihilate outside 
the Galactic disk (see Sec. VI. C). 

In the middle panel of Fig. [21] we compare the lon- 
gitude profile of the observed 511 keV emission to the 
one of the hard LMXRBs (emitting in the 20-100 keV 
range], as seen with IBIS /INTEGRA L ffrom lBird et all 
l2QQ7h IWeidenspointner et all (|2QQ8ah notice that the lat- 
ter distribution exhibits a pronounced asymmetry, with 
source number ratio N(l < 20°) /N (I > 20°)=1.7, which 
matches well the asymmetry in the 511 keV flux re- 
ported in the same paper. They suggest then that hard 
LMXRBs may be at the origin of the disk emission. 

We n ote that the study of the sa me SPI/INTEGRAL 
data bv lBouchet et all (|2008l . l2010h with different meth- 
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TABLE IX Properties of candidate positron sources in the Milky Way 



Source 


Process 




e + rate 6 


Bulge/Disk c 


Comments 






fMeV") AT 


+ fl0 43 s" 1 ) 


B/D 




l\/TQaaiT7-<=> dnrc ^ A 1 
IVldbblvt; bUdlb. r\L 






n 4 




IV, JD I U . \J Ubcl VcLUlUIlcLliy llllcl I cLl 


Supernovae: 2 ^Ti 


/3^-decay 


rvj 1 


n q 
U.o 


<.u.z 


AT: Robust estimate 


QlVTo- 56 isj; 
oiNia. iM 


/3^-decay 


_ 1 
~1 




<u.o 


Assuming / e + esc =0.04 


Novae 


/3 + -decay 


-1 

~1 


U.U2 


<(J.O 


Insufficent e + production 


Mypernovae/ljrxl-D: 1N1 


/3 + -decay 




7 


<U.z 


Improbable in inner MW 


Cosmic rays 


p-p 


~30 


0.1 


<0.2 


loo high e^ energy 


LMaKds 


7-7 


~1 




<(J.O 


Assuming L e + ^0.01 L b s ,x 


Microquasars (/iQs) 


7-7 




1 


<U.O 


e + load of jets uncertain 


Pulsars 


7 - 7 / 7 - 7s 


>30 


0.5 


<0.2 


Too high e + energy 


ms pulsars 


7 - 7 / 7 - 7s 


>30 


0.15 


<0.5 


Too high e + energy 


Magnetars 


7 - 7 / 7 - 7s 


>30 


0.16 


<0.2 


Too high e + energy 


Central black hole 


P-P 


High 


? 




Too high e + energy, unless B >0.4 mG 




7-7 


1 


? 




Requires e + diffusion to ^1 kpc 


Dark matter 


Annihilation 


1 (?) 


? 




Requires light scalar particle, cuspy DM profile 




Deexcitation 


1 


? 




Only cuspy DM profiles allowed 




Decay 


1 


? 




Ruled out for all DM profiles 


Observational constraints 




<7 


2 


>1.4 





a: typical values are given, b: e + rates: in roman: observationally deduced or reasonable estimates; in italic: speculative (and 
rather closer to upper limits), c: sources are simply classified as belonging to either young (B/D <0.2) or old(<0.5) stellar 

populations. 



ods does not find significant disk asymmetry. Obvi- 
o usly, the important (and intriguing) observational result 
of IWeidenspointner et all (|2008af ) needs confirmation by 
further observations/ analysis. Assuming the asymmetry 
is real, what might the implications b e? The interpreta- 
tion of lWeidenspointner et all (|2QQ8al ) implicitly assumes 
that: 

i) Positrons annihilate relatively close to their sources, 
such that the annihilation morphology reflects the source 
morphology. 

ii) Among all LMXRBs only the bright and hard 
LMXRBs of the IBIS/INTEGRAL catalogue ([Bird et all 
2007) are important e + contributors. 

iii) Those X-ray bright and hard-spectrum LMXRBs 
have the same average positron yields, which are there- 
fore not correlated to their currently observed (but 
known to be widely varying) X-ray luminosities; in 
that way, the collective e + production of hard-spectrum 
LMXRBs is just proportional to their total number, not 
to their total X-ray brightness. 

Assumption (i) underlies all efforts to match the ob- 
served 511 keV morphology with some particular class 
of sources. However, if it is adopted, and if it is as- 
sumed that the observed disk emission is due to 0.7 10 43 
e + s -i re i easec [ Dv hard LMXRBs, one has to explain 
why the robustly established e + production of 26 Al and 
44 Ti is not detected by SPI. Indeed, the corresponding 
e + production rate is quite high (0.7 10 43 s _1 , with small 
uncertainty) and positrons are released in the dense en- 
vironment of massive stars and CCSN. In comparison, 
positrons from LMXRBs are released away from the disk 



(in view of the ~400 pc scaleheight), i.e. in less dense en- 
vironments, and could travel and annihilate further away 
from their sources than those of massive star radioactiv- 
ity. If both radioactivity and LMXRBs release 0.7 10 43 
e + s -1 , the former should dominate the observed 511 keV 
emission (the latter having a lower surface brightness), 
and n o significant asymmetry should be seen ([Prantzosl 
l2QQ8h , 

Assumption (ii) has been criticized in 
Bandvopadh vav et al. I (|2009tl . who note that a lower 
sensitivity cut-off than the one of IBIS would lead to a 
different spatial distribution of the hard LMXRBs, in 
view of the steeply rising luminosity function of those 
sources. Besides, in view of the time variability of 
LMXRBs, the present day asymmetric profile (merely 
a snapshot in time) does not guarantee that the same 
morphology characterizes the total number of hard 
LMXRBs that may contribute to e + production during 
the e + lifetime (~10 6 yr) . Notice a lso th at in the 4th 
IBIS source catalogue (jBird et al I l2010h there is no 
strong evidence for a LMXB distribution asymmetry in 
the Galactic plane. 

Finally, assumption (iii) is far from obvious. This 
assumption certainly applies to e.g. SNIa, which are 
assessed to have an average 56 Ni yield of 0.7 M and 
to constitute a relatively homogeneous class of objects. 
One may certainly imagine that LMXRBs also pro- 
duce, on average, the same yield of positrons, at least 
on timescales comparable to the positron annihilation 
timescale. However, if LMXRB positrons are produced 
in the inner accretion disks by processes depending on 
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parameters of the binary system (e.g. temperature, de- 
pending on black hole mass) then only a few of those 
systems may be important e + producers; their spatial 
distribution may not be represented at all by the one of 
all hard LMXRBs. 

The morphology of the observed 511 keV emission pro- 
vides also some interesting constraints in the case of dark 
matter particles as positron sources (under the assump- 
tion of negligible e + propagation). An ill ustration of such 
an an alysis is provided in the work of lAscasibar et all 
( 20061 ). who convolved the positron maps predicted for 
various light DM particle scenarios and types of DM halo 
profiles with the response function of SPI. Comparison to 
the data showed that: i) Particle candidates with velocity 
dependent cross section are excluded as the main source 
of 511 keV emission, ii) Fermionic DM candidates are also 
excluded, since they would need to exchange too light 
charged particles, and iii) Decaying dark matter cannot 
be the main source of low energy positrons, because the 
resulting flux profile is too flat, compared to SPI data. 
Notice that this latter feature is a generic property of all 
models involving decaying particles, where the positron 
production (and annihilation) rate is proportional to the 
DM density profile: even "cuspy" profiles, such as the 
NFW (see FigJT2j) do not provide a 7-ray flux profile suf- 
ficiently peaked towards the inner Galaxy. Annihilat- 
ing or de-exciting DM produces positrons at a rate pro- 
portional to the square of the DM density profile (Sec. 
IV.C.2) and leads to a much more peaked 7-ray profile. 
lAscasibar et all (|2QQ6h found that light scalar annihilat- 
ing particles remain a possible candidate, provided the 
DM halo is at least as cuspy as the NFW profile with 
7 ~ 1 (see bottom panel of Fig. [21]); however, as stressed 
in Sec. IIII.E| astrophysical evidence favors flatter DM 
halo profiles. 

The proximity of the Galactic center and the ex- 
pected high density of DM particles there, make it the 
prime target for the detection of all kinds of radia- 
tion emitted indirectly by DM (either decaying, anni- 
hilating or de-exciting). However, because of the un- 
certainties presently affecting the density profile of DM 
haloes (see Sec. III.E) and the possible contamination 
of the signal by more conventional astrophysical sources, 
other potential targets have been seeked for. The dwarf 
spheroidal (dSph) satellites of the Milky Way, with their 
high mass /light rati o and relative proxi mity, may con- 
stitute such targets. iHooper et all (|2004l ) suggested that 
the light DM hypothesis could be tested on the nearby 
(25 kpc) dSph galaxy Sagittarius, which appears to be 
dominated by dark matter. A search for the expe cted 
annihilation signal at 511 keV ( Cordier et all . [2004) was 
unsuccessful. 



already emphasized, e + production rates of all those 
sources are extremely uncertain (except those of 26 Al, 
44 Ti and GCRs) and the values listed above should be 
considered as optimistic rather than typical ones. Only 
in the case of novae may the estimated production value 
be used to eliminate those sources as important e + pro- 
ducers. Source morphology and high energy of pro- 
duced positrons appear to exclude pulsars, magnetars 
and GCRs as major contributors to the observed 511 
keV emission from the bulge. Source morphology alone 
would exclude hypernovae and GRBs. The high energy 
of produced positrons disfavors ms pulsars, as well as p-p 
collisions from any source (micro-quasars, LMXRB jets, 
the central SMBH). This still leaves several potentially 
important e + contributors, but none of them has the ob- 
served morphology of 511 keV emission. 

Thus, assuming that positrons annihilate near their 
sources, one has to conclude that 

i) Either an unknown class of sources dominates e + 
production, or 

ii) Positrons are produced by a combination of 
the sources of Table El e.g. (a) 26 Al+ 44 Ti for 
the disk and dark matter for the bulge, or (b) 
26 Al+ 44 Ti+LMXRBs(or microquasars) for the disk and 
the bulge plus a contribution from the central SMBH for 
the inner bulge, or (c) some other combination. 

Alternative (ii) bears an interesting "philosophical" is- 
sue: how is it possible that two widely different classes of 
sources have so similar e + yields (to within a factor of a 
few), such as required to fit the observations? However, 
such " coincidences" are not unusual in astronomy 25 these 
days. 

A more important issue arises in solutions involving 
the central SMBH as the main e + producer within the 
Galaxy's bulge: If its positrons can diffuse to kpc scales in 
the dense environment of the inner bulge, then positrons 
should diffuse to even larger scales in the less dense envi- 
ronment of the disk; this would be even more true outside 
the spiral arms, where most of the SNIa/LMXRB//iQ 
positrons are expected to be released. Positron escape 
from the disk to the halo would alleviate the morphology 
problem, by reducing the disk 511 keV emissivity and 
thus increasing the B/D ratio of those classes of sources. 
However, although some of the basic physical processes 
underlying e + propagation are well understood, there is 
no clear global picture of how far positrons can propa- 
gate in the magnetized, turbulent ISM of the Galaxy. We 
turn to those two issues in the next two sections. 



4. Summary of candidate sources 

The main features of the candidate e + sources dis- 
cussed in this section are summarized in Table IIXI As 



25 Among the most famous "coincidences" are: (i) the contributions 
of baryonic and non-baryonic matter, as well as those of dark 
matter and dark energy, to the cosmic density; (ii) the solar 
abundances of s- and r- nuclei (both of ~10 -6 by mass fraction); 
(iii) the approximately equal contributions of CCSN and SNIa 
to the solar Fe content. 
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FIG. 22 The processes leading to gamma-ray production 
from positron annihilation (taken from lGuessoum et aZ.Lll99lL 
slightly modified to account for the annihilation in flight of 
relativistic positrons) . 



V. POSITRON INTERACTIONS WITH MATTER AND 
ANNIHILATION 

Positrons are initially produced with kinetic energies 
higher than those of the ISM (see previous section). Most 
of them slow down to low energy before they annihilate 
with bounded or free electrons of the ISM. However, when 
their initial kinetic energy is above a tenth of MeV, a sig- 
nificant fraction of them may annihilate in-flight. Figure 
Fig. [22] summarizes the processes that lead to gamma- 
ray production from positron annihilation. The follow- 
ing sections list the interactions that are responsible for 
the energy losses of positrons in the ISM and present the 
different ways in which they annihilate with electrons. 



A. Energy losses 

As charged leptons, positrons interact via the electro- 
magnetic force with all basic constituents of the ISM, 
namely: electrons, ions, atoms, molecules, solid dust 
grains, photons and magnetic fields. Since their initial 
kinetic energy is generally larger than the kinetic energy 
of the targets in the ISM, positrons lose energy in these 
interactions. The energy loss rate and the kind of inter- 
action depend on the energy of positrons and the density 
of target particles. Fig. [23] presents the energy loss rate 
as a function of positron energy; the contributions of each 
type of interaction are shown separately. 

Ultra-relativistic positrons (E>10 GeV) lose their en- 
ergy mainly by inverse Compton scattering with cosmic 
microwave background (CMB) photons and interstellar 
radiation fields. When the interaction occurs with an 
isotropic photon gas in the Thomson scattering regime 
(i.e. hv <C m e c 2 ), the energy loss rate (in eV/s) can be 
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FIG. 23 Energy loss rate for positrons in ISM conditions. For 
synchrotron losses, the pitch angle is here taken as n/2. 



calculated from ( Blumenthal and Gouldl . IT970I ): 



dE 
~dt 



2.6 x 10" 14 w rad f ^ 



(10) 



IC 



where ?/ ra d is the radiation energy density (eV/cm 3 ), 
and /3 = v/c its velocity relative to the speed of light. 
The radiation energy density depends on the position of 
the positron in the Galaxy; it ranges from 0.26 eV/cm 3 

CMB) to 11.4 eV/cm 3 in the Galactic Centre region 

iMoskalenko et all [20061 ). 

Ultrarelativistic positrons also lose their energy by 
emitting synchrotron radiation when they spiral along 
a magnetic field line. The energy loss rate depends on 
the magnetic field intensity, the positron's kinetic en- 
ergy and the pitch angle 26 . I ts expression (in eV/s) is 
(jBlumenthal and Gouldl . 1 197Qh : 



^f) = -9.9 x IO" 16 B 2 7 2 p 2 sin 2 (a) , (11) 
at J SY 

where B is the magnetic field (in fiG) and a is the 
pitch angle. For a positron moving through a ran- 
domly oriented magnetic field, the mean energy loss rate 
is obtained by replacing sin 2 (a) by its average value 
2/3. The synchrotron energy loss rate is proportional 
to the square of the magnetic field while the inverse- 
Compton one is proportional to the radiation energy den- 
sity. Therefore, synchrotron losses dominate for ultrarel- 
ativistic positrons in an environment where B > 6.3/iG 
y^rad, e.g. near Pulsars, the Galactic Centre region, etc. 

In the 1-10 GeV energy range, positrons lose their en- 
ergy mainly by emitting bremsstrahlung radiation in in- 
teractions with ions, electrons and atoms. The energy 



26 The pitch angle of a charged particle is the angle between the 
vectors of the particle's velocity and of the local magnetic field 
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loss rate depends on the target mass, charge, and den- 
sity. For relativist ic positrons, it is equivalent to the 
electron one. Methods for calculating the average en- 
ergy loss rate an d the brem sstrahlung energy spec trum 
are described in Havakawal (|l969h . lEvansI (|l975h and 
iBlumenthal and Gouldl (|l970h . A recent update, with a 
more accurate treatment of the differential cross section 
for electrons with midly relativis tic en e rgies in a fully 
ionized plasma, was presented in lHaugl (|2QQ4l ). A good 
approximation of the bremsstrahlung energy loss rate (in 
eV/s) of relativi stic positrons in a fully ionized gas is 
(iGinzburd . Il979l p 408): 



dE_ 
~dt 



3.6 x 10 _ii Z(Z+1) n<y 



BR 



Zn(2 7 ) - - 



(12) 

where Z and n are the nuclear charge and the number 
density (in cm -3 ) of the ion, respectively. In a neutral 
hydrogen gas, th e energy loss rate can be estimated via 
(|Ginzburj . 119791 p. 386): 



~dt J 



A n 7, 



(13) 



BR 



where A= -4.1x10 10 for hydrogen and A=-l.lxl0 9 
for helium, respectively. 

Below 1 GeV, positrons lose their energy mainly via 
Coulomb scatterings with free electrons and/or inelas- 
tic interactions with atoms and molecules. The for- 
mer process is a continous energy loss, whatever the en- 
ergy of positrons. At high energy, the target electrons 
can be considered at rest an d the energy lo ss rate de- 
pends mostly on their density. iDermerl ( 1985h calculated 
the rate of e-e Coulomb collisions in relativistic ther- 
mal plasmas. He also treated the case of collisions in 
cold plasmas. T his energy loss can be approximated by 
(jGinzburd . fl979l p 361): 



TABLE X Energy thresholds of inelastic reactions produced 
by positrons. 



Process 
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where n is the neutral atom density, Z is the number 
of electron of the atom and I is its ionization potential 
(e.g. 13.6 eV for H and 24.6 eV for He). Below 1 keV 
the interaction between positrons and atoms/molecules 
should be estimated via Monte Carlo simulations since 
positrons release a large fraction of their energy in one 
interaction and, below ~ 100 eV, they can pick-up an 
electron from an atom or a molecule to form a positron- 
ium "in flight" (see Table IXl and the next section). In this 
case, the cross sections for such collisio ns should be used 
in ev aluating the energy loss rate (see iGuessoum et all 
l2QQ5h . The Monte-Carlo method to calculate the inter- 
action probabilit y as a function of the positron's energy 
was presented by iBussard et al\ (|l979f ) . 

Equations [T01IT51 allow one to estimate to a good 
aproximation energy losses of positrons. More accu- 
ra te expressions, valid to energies > 1 GeV, are given 
in lStrong and Moskalenkol (| 19981 ). 



f dE\ 
K dt J cou 



-7.7 x 10- 



ne 



In 



73.6 



(14) 

where n e is the electron density. 

At low positron energy (E<10kT, where T is the am- 
bient temperature), the electrons of the ISM cannot be 
considered at rest anymore. The energy loss rate de- 
pends on their temperature a n d den sity in the plasma 
and is gi ven bv iBook and Alii (|l975l ) and Huba (2006) 
(see also lMurphv et all l2005al) . 

Inelastic collisions of positrons with atoms and 
molecules can be considered as a continous process for 
a positron energy >1 keV, and the energy loss rate can 
be evaluated using the Bethe-Bloch formula. The ioniza- 
tion loss is larger than the e xcitation loss. T his energy 
loss can be approximated by ( Ginzburgl 1 19791 p 360): 



B. Annihilation in flight 

In this section, we present two kinds of annihilation 
in flight: (1) direct annihilation of relativistic positrons 
with electrons and (2) annihilation via positronium pro- 
duced in interactions of non relativistic positrons with 
atoms / molecules. 



1. Direct annihilation in flight 

When high-energy positrons (j> 10 keV) slow down, 
they may annihilate in flight with free or bound elec- 
trons. The energies of the two photons emitted in this 
process are strongly shifted as per the Doppler effect. 
This produces a continuous spectrum in the energy in- 
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FIG. 24 Gamma-ray spectra from the annihilation in flight 
in the ISM for various initial kinetic energies of positrons. 



terval mc 2 /2 < E 1 < E + mc 2 /2 with E the total energy 
of the positron. 

The probability that a positron with an initial kinetic 
energy Eq annihilates in flight before reaching an energy 
E is: 



P(E ,E) = 1 



exp 



v(E')a a (E')dE' 
\dE'/dt\ 



(16) 

where v and dE f /dt are the positron velocity and the 
energy loss rate, respectively; n e is the density of target 
electrons, and a a the annihilation cross section, which 
can be est i mated for kinetic energies larger than 75 keV 
via ([Diracl . fl930h : 



7 2 + 4 7 +l 7 ( . r~ 2 -v 

3 /n(7 + VT - 1) 



7 + 1 [ 7 2 -l ' v ' ~' Vj 2 ^ 

(17) 

where r e is the classical electron radius (r e = e 2 /m e c 2 ); 
for an evaluation of this cross section below 75 keV, 
see Gould (1989). The probability P(E ,E) depends 
on the energy loss rate and consequently on the phys- 
ical conditions of the interstellar medium in which the 
positron propagates. Fig. [24] presents the spectra of 
gamma-rays emitted by in-flight annihilation of relativis- 
tic positrons slowing down in the interstellar medium, for 
several initial kinetic energies. Fig. [251 shows the fraction 
of positrons annihilating in flight as a function of the 
inital kinetic energy of positrons, in both a neutral and 
a fully ionized medium. It is negligible (<j4%) for en- 
ergies lower than 1 MeV. Above ~ 1 GeV, the fraction 
does not change because the energy loss rate is so large 
(see Fig. [23]) that positrons do not have enough time to 
annihilate in flight at these energies. 

The implications of e + annihilation in flight for the 
observed Galactic MeV emission are further analyzed in 
Sec. IV.E. 
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FIG. 25 Probability for a positron to annihilate in flight as a 
function of its initial kinetic energy, in a neutral medium (solid 
line) and an ionized medium (dashed line). Y e represents the 
ionization fraction. 



2. Positronium formation in flight 

When positrons have kinetic energies lower than ~ 100 
eV, they can pick up an electron from an atom or a 
molecule to form a positronium in flight while they slow 
down. This reaction, also called charge exchange, is en- 
doenergetic. It can happen as long as the kinetic en- 
ergy of the positron is larger than the charge exchange 
threshold energy of Ps formation with the given atom 
or molecule (see Table. [X]). This threshold energy is 
equal to the ionization potential of the atom reduced by 
the binding energy of the Ps (6.8 eV). The cross sec- 
tions for positronium formation by charge exchange were 
measured mostly in rare gases, particularly in helium, 
and in molecular and a tomic hydrogen (see references in 
iGuessoum et aZ.L l2QQ5h . 

The fraction of positrons that form a positronium in 
flight is obtained by Monte Carlo methods which consist 
in simulating the interactions of positrons with atoms 
and molecules on the basis of the cross sections and en- 
ergy loss mechanisms presented in the previous section. 
This fraction depends on the density, on the temperature, 
and strongly on the ionization fraction of the medium in 
which positrons slow down. Fig. [26] shows the fraction of 
positrons forming positronium in flight a s a function of 
the ion ized fraction in two types of media ([Murphy et all 
l2005af ). The fraction of positrons that form Ps decreases 
quickly with increasing values of the ionized fraction not 
only due to the reduction of the density of neutral H but 
also because the energy loss rate increases quickly with 
the ionized fraction (see Fig. [23]). It makes the positrons 
slow down so rapidly that they d o not have time to ex- 
change charge with H. Me asured ([Brown and Leventhall 
1 19861 ; iBrown et all [l984) and calculated values of the 
fraction of positrons forming positronium in flight in H, 
He and H2 are summarized in Table [XTl 

When the Ps is produced in flight, its kinetic energy 
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FIG. 26 Fraction of positrons forming positronium in flight 
by charge exchange with atomic hydrogen as a function of 
the ionization fraction (Y e ) in a warm component of the in- 



terstellar medium (electron density: n e = 0.1 cm 3 . 
temperature: T e = 8000 K) and in solar flare (n e = 
cm" 3 , T e = 1.16 x 10 4 K). 
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is equal to the energy of the positron minus the thresh- 
old of the charge exchange reaction. Consequently, the 
energy of the photons emitted in the annihilation are 
shifted as per the Doppler effect. The spectral shape of 
the 511 keV line emitted in the annihilation of para-Ps 
was also derived from Monte Carlo simulations by ex- 
tracting the kinetic energy distri bution of the produced 
Ps. The most recent calculations ( Guessoum et a/ll2005l ) 
yield widths of 5.8 keV, 6.4 keV and 7.4 keV for H, H 2 and 
He, respectively. The calculated line widths for He and 
H2 are in good agr eement with previous measurements 
([Brown et all\l984 . 



For large gas density, the time scale for a Ps to col- 
lide with an ambient particle X (atom, electron, pho- 
ton) may be comparable to or lower than the lifetime 
of the ortho-Ps. In that case, the ortho-Ps can be de- 
stroyed ( 3 Ps + X — > X + e + + e~) or converted into 
a para-Ps (spin flip - 3 Ps + X X + r Ps). These 
processes tend to reduce the contribution of the an- 
nihilation via the ortho-Ps state. They are expected 
to occur in the solar atmosphere during flares (see 
iGuessoum er^.ll997HMurphv et a/.U2005ah and in nova 
envelopes ([Leising and Clavtonl Il987 ). In such media, 
the interactions of Ps with components (atoms, electrons, 



TABLE XI Fraction (in %) of positrons forming positronium 
in flight, in totally neutral media. 



References 


H 


H 2 


He 


Bussard et al. (1979) 


95 


93 




Brown & Leventhal (1986) 




89.7 ± 0.3 


80.7 ± 0.5 


Wallvn et al (19941 


98 


90 




Chaouis et al (1994) 






78 


Guessoum et al. (2005) 


95.5 


89.6 


81.7 



photons) of the plasma must be taken into account in 
evaluating the proportion of annihilations through the 
ortho-Ps and the para-Ps states and is known as "ortho- 
Ps quenching" . 



C. Thermalisation 

Once the positrons have come down to energies similar 
to those of the ambient medium, they start to "ther- 
malize", i.e. their energy distribution relaxes to the 
Maxwellian function which characterizes the interstellar 
gas (or plasma) 27 . The ISM is usually considered to con- 
sist of a few phases, each with rather well defined physical 
characteristics (temperature, density, ionization fraction) 
- see Table [IVJ The timescale needed for the energetic 
positrons to relax to the ISM Maxwellian distribution is 
to be compared to the timescale for subsequent annihi- 
lation processes; if the former timescale is longer than 
the latter, it would be incorrect to assume a Maxwellian 
distribution for both the positrons and the ISM when 
calculating the e + annihilation rates. 

To tackle this question one may simply estimate the 
relaxation timescale or perform a full statistical-physics 
treatment. The former can be simply done by using the 
energy loss rate: r = — (^^f ) _1 or ( a simpler but cruder 
evaluation) r _1 = R =< nav >, taking for the cross 
section some typical inelastic scattering value. The first 
formula (using values in Fig. [23] for the energy loss rate) 
gives r ~ 6 x 10 7 seconds for 1 keV positrons, while us- 
ing the cruder approach gives ~ 2 x 10 T seconds (taking 



a ~ 10 -16 cm 2 , E + ~ 1 keV, and n ~ 1 cm -3 ); estimates 
for the relaxation of positrons from their initial energies 
of ~ MeV , is more complicated, both because the par- 
ticles undergo many different processes and the overall 
cross section is difficult to estimates, but roughly the two 
simplistic approaches give timescales < 10 12 seconds. It 
must then be noted that all of these values are much 
less than typical annih ilation timescales in the ISM (see 
IGuessoum et all 120051 ). which range around 10 12 — 10 14 
seconds. 

The sophisticated statistical-physics approach is much 
more rigorous and conclusive, though complicated. A 
few authors have attempted such tre atments, ranging 
from very broad, highly theoretical ( Wolfe and Melial 
2006) t o others cove r ing a specific area of appli- 
cation ( Crannell et all [l976), the latter focusing on 



low-energy (50 eV) p ositrons produced a nd a nnihi- 
lated in solar flares . iDermer and Liand (1989) and 
iNavakshin and M elia (jl998h have performed thorough 
statistical-physics treatments for the thermalization of 



27 It is understood that unless the background medium is in ther- 
modynamical non-equilibrium (occurrence of irregular heating, 
sudden energy losses, etc.), it will have adopted a Maxwellian 
energy distribution characterized by a temperature T. 
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high-energy electrons (and secondarily of positrons), as- 
suming the interaction of the injected particle s takes 
place with a relativistic electron-proton plasma. iBarind 
(|l987h . modifying a simple treatment bv ISpitzerl ((1956). 
showed how to calculate the relaxation time r(E e ) of 
an electron for any temperature (i.e. relativistic or 
not). The statistical-physics method uses a Fokker- 
Planck equation to follow the energetic particles from 
their injection (at a given time and energy) to their 
thermalization either with the background plasma/gas or 
with each other, taking into account the various inelastic 
scattering and loss processes; in this approach, one eval- 
uates the energy and dispersion coefficients (^-^ and 

n~ d ^(U " ? respectively). ICheng et all (|2006f ) performed 
such a treatment for high-energy (30 MeV) positrons, 
presumably produced by proton-proton collisions at the 
Galactic center and approaching thermalization with the 
ISM gas and plasma conditions (T ~ 10 4 K and 10 6 K, 
respectively). 

These various approaches, heuristic or elaborate, have 
shown that the positrons do indeed thermalize before the 
annihilation processes become important, which makes 
valid and legitimate the usage of Maxwellian distribution 
functions in the calculations of the rates of annihilation 
and other processes that the positrons undergo. 



D. Annihilation 

Positrons annihilate by various processes during the 
two stages of their "lives": a) during the slowing-down 
time (the process that is referred to as "in-flight" anni- 
hilation, described in section V.B); b) after thermaliza- 
tion with the ISM. In this section we consider the latter 
stage and its processes. As explained in sections V.A 
and V.B, while positrons lose the bulk of their energies 
from 1 MeV to ~ 100 eV, their probability of un- 
dergoing "charge exchange" ("picking up" an electron 
from an atom or molecule and forming a positronium) 
increases steadily. Thus, by the time they have thermal- 
ized with the medium, the probability f\ that they will 
have formed a Ps can be as high as 95 %, depending on 
the ionization state of the medium (see Fig. [26]). 

Depending on the physical conditions of the medium 
(ionization state, temperature, composition), positrons 
may undergo a variety of annihilation processes, which 
are listed here by order of decreasing strength of 
their cross sections (when the physical conditions allow 
them to occur): 1) charge exchange (with atoms and 
molecules); 2) radiative (re) combination with (free) elec- 
trons; 3) direct annihilation with free electrons; 4) di- 
rect annihilation with bound electrons (of atoms and 
molecules). For dust grains, the cross section for colli- 
sion and annihilation (see below), may be larger than 
those of some of the above processes; however, when the 
abundance of dust grains is taken into account, the re- 
action rate for this process turns out to be smaller than 
the others', except in rare cases. 



u 



10° 

io- 1 
io- 2 
io- 3 

IO" 4 

io- 5 
io- 6 

lOi 



"e + + H - 
"e + + H - 
"e + + H- 



• e + + H + + e- 

• e + + H* 

• e + + H* 



"e + + H^Ps + H + 
"e + + e" -» Ps + y 
~e + + e - -» 2 y 
-e + + H-2Y + H + 



■ Maxwellian distribution 
T = 8000 K 



•r? 
i 



10- 



io- ] 



10° 10 
E(eV) 



10 2 



FIG. 27 Cross-sections for ionization, excitation, charge ex- 
change, radiative recombination and direct annihilation inter- 
actions of positrons with atomic hydrogen and free electrons. 
Also shown is the Maxwellian distribution for a temperature 
of 8000 K (in arbitrary units). 



The key features of these annihilation processes (see 
Fig. [27]) are the following: 

a) Charge exchange (with H, He, and H2): although 
difficult to measure, the cross section for this(ese) 
process (es) has been obtain e d by several experimental 
groups (see iGuessoum etHUL 120051 for references); the 
important feature of that process is its threshold energy 
(6.8, 17.8, and 8.6 eV, respectively), implying that this 
reaction (which has by far the highest cross section of all 
positron annihilation processes) can only occur at tem- 
peratures larger than several thousand K; hence this can- 
not occur in the cold media of the ISM, and of course it 
cannot take place in the hot ISM phase either (because 
the medium is completely ionized). 

The width of the line resulting from the annihilation 
of thermal positrons by charge exchange (decay of Ps) 
with H or He can be derived from the kinetic energy 
distribution of the produced Ps. This distribution is ob- 
tained simply by computing the charge exchange rate for 
a population of positrons that follows a Maxwellian dis- 
tribution. The calculated widths in a warm medium (T 
= 8000 K) a re 1.16 keV and 1.22 keV, for H and He, 
respectively ([Guessoum et aZl l2QQ5h . 

b) Radiative (re) combination with free electrons: the 
cross section for this process is too small (~ 10 -20 cm 2 at 
1 eV) to be measured experimentally; one then has to rely 
on theoretical calculations, such as the determinations 
(by dif ferent approac hes) made bv lCrannell et al. I (119761 ) 
and bv lGouldl (I1989D . 

c) Direct annihilation with free electrons: this process 
has an even smaller cross section (about an order of mag- 
nitude less than the previous one at temperatures less 
than about 10 5 K), hence it is only important in the 
hot phase of the ISM; the c r oss se ction has also been 
estima ted by ICrannell et all (|l976f ) and iBussard et all 
(|l979h from the early theoretical work of lHeitlerl (|l954). 
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The width of the line resulting from both radiative 
combination and direct a nnihilation of positr ons and 
electrons was calculated bv liCrannell et al 1 (|l976h using a 
simple argument of thermal broadening due to the pair's 
center-of-mass motion; they obtained the simple expres- 
sion Y rC)dae = 1.1 x (T/10 4 ) 1 / 2 keV, which applies to both 
processes. 

d) Direct annihilation with bound electrons: the cross 
sections (for H, He, and H2) are the weakest of all 
(see Fig. |27j), but they become important by default at 
very low temperatures (in the cold phases of the ISM) 
where free electrons do not exist and the charge exchange 
(with atoms and molecules) cannot take place due to the 
threshold energies, which are (at about 10 eV) much too 
large compared to the particles' average thermal energies 
of ~ 0.01 eV. The first work which h as performed a calcu - 
lation of this cross section is that of lBhatia et~al\ (|l977[ ): 
more detained calculations, taking into account short 
range interactions between the positron and the target 
electron (e.g. virtual forma ti on of posit ronium), were 
perfor med by llgarashi et all ( 20021 ) and lArmour et all 
(|l99Qh . for positrons colliding with H and H2, respec- 
tively. 

The widths of the lines resulting from these processes 
have bee n measured b y Brown and Leventhal (1986) for 
H and bv llwatal (|l997h for He and H2; the values obtained 
were: 1.56 keV, and 2.50 keV, and 1.71 keV for H, He, 
and H2, respectively. These values have a very weak 
dependence on the temperature. 

e) Annihilation on dust grains: the importan ce of 
this process was first pointed out bv IZurekl (|l985l ). who 
stressed the important effect this process would have 
on the Ps formation fraction, which is a quantity that 
can be inferred from observational data (see section 
II. A) and thus represent an important constraint on 
models fGuessoum et all (|l99lh then refined the calcula- 
tion of the rate, adding electric-charge and positron-grain 
reflection effects, on the one hand, and spectral consid- 
erations (li ne width and effe c t on the overall calculated 
spectrum). IGuessoum et all (2005) have done the most 
extensive astrophysical treatment of this process to date, 
despite the dearth of some crucial information on the 
processes, considering the materials that constitute the 
dust grains; in particular, the widths of the lines result- 
ing from the annihilation of positrons inside the grain 
(after capture) and from the decay of the positronium 
which is formed in/on the grain and ejected out have 
been evaluated to Tp S) i n « 2.0 keV and Tp SyOUt ~ 1.4 
keV, respectively; this is highly important in that it af- 
fects the amount of dust one will infer from the galactic 
positron annihilation l ine spectra (see the discus sion in 
IGuessoum et all I2QQ4J and IGuessoum et a/.Ll2005l ): note, 
incidentally, that due to the fact that Ps inside the grain 
undergoes "pick-off" annihilation, it always gives two 
photons, never three. 

The corresponding reaction rates, taking into account 
the (Maxwellian) energy distribution of the particles 
(positrons, electrons, atoms, or molecules) and the abun- 



dance and density of each species, are given by: 



r v = n < av > 



-F 



3/2 



- E ' kT a(E)vdE . 



(18) 

The reaction rates then allow one to determine the frac- 
tion of positrons which annihilate through each process, 
fp = r p /Y,r p (the index p generically referring to a pro- 
cess), and these fractions are then used to determine the 
spectrum of emission in a given physical medium. 

The spectrum of gamma-ray emission includes contri- 
butions from various processes: each one of them consists 
of either a Gaussian function describing the line emis- 
sion at 511 keV with a given line width (FWHM denoted 
by T) or a Gaussian (denoted below by G(E,E' ,T), of 
variable E' and center E) and an ortho-positronium con- 
tinuum (at < E < 511 keV), the latter given by the 
lOre and Powelll (| 19491 ) function Pt(E). The spectrum is 
then given by: 



S(E) 



= J dE' 3x ^P t (E')+2x ^5(E'-E ) 

{ X x f\,H/H 2 G(E, E' , Tif jH / H2 ) 
+ Yxf hHe G(E,E\T^ He ) 
+ (1 - X f hH/H2 - Y f liHe ) 

x [fce,H/H 2 G(E,E' ,T ce ^ H / H2 ) 

+ fee, He G(E,E' \T ce ^ He ) 

+ free G(E,E' \T rce ) 

+ fgr,out G(E, E' ,T gr ^ out )] } 

+ 2 (l - X fi i H/H 2 ~ Y fl,He) [fdae G(E, E , T dae ) 
J ? da,H I H 2 

G(E,E , 

1 da,H/H 2 ) 

+ fda,He G(E, Eq, 

J- da, He ) 

+ fgr,inG(E 1 EQ 1 T gr ^ n )) 1 (19) 

where X and Y are the relative abundances of H(H2) and 
He (90 % and 10% respectively, by number) and fi is 
the fraction of positrons forming positronium in flight. 
The spectra are presented in Fig. [28] for each phase of 
the ISM. In neutral media the line is broad due to the 
annihilation of Ps formed in flight. The width of the line 
is ~1 keV in the warm ionized phase where positrons 
annihilate mainly by radiative recombination with elec- 
trons. The contribution of the annihilation of positrons 
in grains is negligible in all the media except in the hot 
phase where it produces the ~ 2 keV width line super- 
imposed on the broad line (~ 11 keV). The latter results 
from positrons that annihilate directly or via the radia- 
tive recombination process with electrons at a T ~10 6 
K. 

The "global" spectrum of annihilation of positrons in 
the ISM is then constructed by combining the spectra for 
each phase, considering the relative contributions (densi- 
ties and filling factors) of each. For this, a mo del of the 
ISM, like those briefly sketched in Sec. Ill (e.g. iFerrierel 
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FIG. 28 Annihilation spectra for the five ISM phases. 
Adopted temperatures (see Table HV|) are 10 K (Molecular), 
80 K (Cold), 8000 K (WNM and WIM) and 10 6 K (Hot) and 
ionization fractions are for the neutral phases (Molecular, 
Cold and WNM) and 1 for the ionized phases (WIM and Hot). 



119981 : iFerrierel Il999l l200lt iMcKee and Ostrikert Il977h is 
needed. The resulting "global" spectra can then be com- 
pared with the observational data. Of course, it would 
be at least as interesting and useful to compare the in- 
dividual phase spectra with observational data, but this 
has been impossible so far, because of insufficient spatial 
resolution and sensitivity of the detctors. 



E. Spectral analysis of observed emission 

The previous sections described the annihilation pro- 
cesses and the resulting characteristics of the possible 
annihilation emissions. One must distinguish the emis- 
sion produced by the in-flight annihilation of relativis- 
ts positrons, which is characterized by a continuous 
spectrum in the MeV domain (see section V.B.I) from 
the 511 keV line and ortho-Ps continuum emissions pro- 
duced in the annihilation of low-energy positrons (sec- 
tions V.B.2 and V.D). 

As already discussed in Sec. II. B. 2, the observed MeV 
continuum in the direction of the inner Galaxy signifi- 
cantly constrains the energies of injected positrons: they 
have to be lower than a few MeV, otherwise the con- 
tinuum emission would be much higher than observed 
by COMPTEL (see Fig. El and EJ). This allows one 
to eliminate several classes of sources in the steady-state 
regime, such as pulsars, ms pulsars, magnet ars, or ener- 
getic proton collisions (either from cosmic rays or from 
the central black hole), as major positron source s . The 
sam e argument was us ed by iBeacom" and Yiiksell (2006) 
and ISizun et~al\ ([20061 ) to constraint the mass of the de- 
caying or annihilating dark matter particles which could 
be the sources of positrons in the spheroid (see section 
IV.D) to lower than a few MeV. Note that the MeV 
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FIG. 29 Confidence regions of the fit of the temperature and 
ionization fraction to the SPI data obtained after one year of 
mission. The best fit values are T = 7800 K and y e =0.1. 



continuum does not constrain the mass of de-exciting 
dark matter particles, as disc ussed in Sec. IV.D. More- 
over, [Chern^shov^t^ (|2009h showed that the injection 
of positrons with initial kinetic energy higher than sev- 
eral GeV is allowed if the injection is non-stationary 
(e.g. through intermittent emission from the central 
black hole) and if the magnetic field is higher than 0.4 
mG in this region (see also Sec. IV. 5. a). 

Concerning the annihilation of low-energy positrons, 
the shape of the annihilation line and the relative inten- 
sity of the ortho-Ps continuum are closely related to the 
abundances and thermodynamical conditions (density, 
temperature, ionization fraction) of the plasma in which 
positrons annihilate. The broadening of the 511 keV line 
induced by bulk motions of the gas in which positrons 
annihilate is not taken into account since we do not ex- 
pect a significant Doppler shift due to Galactic rotation 
and/or turbulence in the Galactic Centre region (Av <; 
100 km/s, AE <^ 0.17 keV). Consequently, the spectral 
characteristics of the annihilation emission offer valuable 
information on the physical conditions of the ISM where 
positrons annihilate. 

The observed spectrum can be simply characterized by 
the sum of its independent components: the Gaussians 
(to describe the 511 keV lines) and the ortho-positronium 
continuum. The fraction of positrons annihilating via Ps 
(fp s ) is derived from the ratio of the ortho-positronium 
flux to the total line flux (I27) using Eq. 03 The 

measured characteristics (widths, fp s ) can then be com- 
pared with what we expect from the physics of the anni- 
hilation of positrons in the ISM (section V.D). 

From the OSSE data, iKinzer et all (|l996h inferred 
a positronium fraction fp s ^0.97±0.03 in the Galac- 
tic Centre region. Measurements with the Ge detec- 



tor TGRS onboard the WIND miss ion (1995-1997) gav e 
a compatible value of 0.94±0.04 ([Harris et all [l998). 
From the line width - (1.8±0.5) keV - and the positro- 
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FIG. 30 Best fit of the spectrum measured by SPI using the 
warm components of the ISM and the Galactic continuum. 
Contributions from the molecular, cold and ho t components 
are not needed to explain the data (|Jean et all 120061 ) . 



nium fraction measurements, Harris et all Jl998) con- 
cluded that a scenario in which annihilation does not 
occur either in cold molecular clouds or in the hot phase 
of the ISM is favored. Using preli minary SPI data o f 



Jean et all (l2QQ3h and TGRS data of lHarris et a/1 (l 
it all (|2QQ4l 



Guessoum et all (|2004l ) showed that the bulk of the anni- 
hilation occurs in warm gas. However, the two groups did 
not exclude that a significant fraction of the annihilation 
may occur in hot gas and in interstellar dust. 

Another approach consists in fitting annihilation mod- 
els to the observed spectrum. Churazov et al. (2005) 
fitted the temperature and the ionized fraction of the gas 
where the annihilation occurs with a measured spectrum 
based on SPI observations of the Galactic centre region. 
They inferred from their analysis that the spectral pa- 
rameters of the emission can be explained by positrons 
annihilating in a warm gas with a temperature ranging 
from 7000 to 40000 K and an ionized fraction > 1 %. 
However, they did not exclude a combination of warm 
and cold gases. When a similar analysis is performed 
using a spectrum measured with three times more expo- 
sure, a temperature T = (7.81q 5) xlO 3 K and an ionized 
fraction of (10±2)% (see Fig. [29]) are obtained (Jean, un- 
published). 

Instead of fitting t he temperature an d the ionized frac- 
tion to the SPI data, Uean et all (|2006[ ) adopted the spec- 
tral models for the different ISM phases (see Fig. 130]) and 
adjusted the phase fractions fi (with z={molecular, cold, 
warm neutral, warm ionized, hot}) so as to obtain the 
best fit to the spectrum measured by SPI. This model is 
described by: 



where Si(E,x gr ) is the normalized spectral distribution 
(in keV -1 ) of the annihilation photons in phase i, I e + e - is 
the annihilation flux (photons s _1 cm -2 ), and x gr repre- 
sents the fraction of dust grains and allows for uncertain- 
ties in dust abundance and positron-gr ain reaction rates 
(x nr = 1 in the standard gr ain model o f Guessoum et all . 
2005). With this method, Llean et all (|2006l ) found that 



SlSM(E) = I e+e - X^/ t X Si(E,Xgr) 



(20) 



49^23 of the annihilation emission come from the warm 
neutral phase and 51^2% from the warm ionized phase. 
While they may not exclude that up to 23% of the emis- 
sion might come from cold gas, they have constrained the 
fraction of annihilation emission from molecular clouds 
and hot gas to be less than 8% and 0.5%, respectively, 
and the c ontribution of gra i ns to be less than 1.2%. 

Finally, iGuessoum et all (|2Q10h examined the annihi- 
lation of positrons on polycyclic aromatic hydrocarbon 
(PAH) molecules in ISM conditions. They showed that 
PAHs would play a significant role only if their abun- 
dances are higher than about 10 -6 (by number), the low 
abundances being compensated by the large enhance- 
ment in the PAHs' (resonant) annihilation cross sec- 
tions. They used the 511 keV spectrum measured by SPI 
(Fig. [30]) to constrain the PAH abundance in the bulge 
and find an upper limit of 4.6 x 10 -T , consistent with 
results obtained from IR observations on other Galactic 
regions. 



VI. POSITRON PROPAGATION IN THE ISM 

While most astrophysical sources are producing rela- 
tivists positrons (Sec. HV]) . the results of the spectral 
analysis of the Galactic 511 keV emission (line width 
AE/E <0.01 and positronium fraction fp s ~97%) imply 
that positrons annihilate at very low energies (Sec. Inland 
V.E). This implies significant deceleration of positrons on 
their way to the annihilation sites, through the various 
energy loss processes summarized in Sec. [Vj Further- 
more, various plasma processes - such as advection due 
to the Galactic wind, adiabatic deceleration, energy gain 
due to stochastic acceleration and diffusion - may affect 
significantly the propagation of positrons and the extent 
of the region of e + annihilation emission. 

There are two distinctly different regimes of positron 
transport in the ISM: collisional and collisionless. The 
former implies Coulomb interactions of positrons with 
particles in various gas phases of the ISM (in the pres- 
ence of radiation and magnetic fields), while the latter is 
essentially due to the scattering off magnetic turbulence 
in the interstellar plasma. Observations of the solar ener- 
getic particles in the heliosphere reveal that their trans- 
port is indeed dominated by the scattering off fluctuating 
magnetic fields (see Sec. VI.B.l). 

In the collisional regime, positrons change their mo- 
mentum and lose their energy through collisions with 
gas particles while propagating along magnetic field lines 
in the ISM. Fig. [31] shows the stopping distances D = 
v(E')dE' / E' (where v(E) is positron's velocity and 
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FIG. 31 Positron stopping distance D in various phases of the ISM (characterized by typical densities uh and ionization 
fractions x) as a function of positron energy E. Horizontal dashed lines display typical sizes of the corresponding ISM phases. 
In left panel, the lower dashed curve displays positron gyroradius r g in a magnetic fi eld of B — 5 jjlQ, a nd shaded aerea displays 
minimum scalelength Xmin of MHD Alfven waves (the latter being calculated in I Jean et all ((2009) ; both r g and Xmin are 
enhanced by a factor of 10 9 . 



E' represents the sum of energy losses of Sec. V.A) in 
the various phases of the ISM. For MeV positrons, D 
is much larger than the typical size of HIM, WIM and 
WNM phases, and comparable to the size of CMM phase 
(horizontal lines in Fig. I3T]) 28 . Therefore, in the colli- 
sional regime, only the CMM phase may be efficient in 
stopping MeV positrons. 

In the case of a magnetized plasma, positrons are spi- 
raling along the magnetic field lines. The gyroradius 
of a positron with Lorentz factor 7 is r g ~ 1.7 10 9 
B~q{^ 2 — l) 1 / 2 cm, where the local mean magnetic field 
B^g is expressed in fiG. For typical values of the Galac- 
tic 5-field (1-10 /iG, see Sec. III.D), r g is many orders 
of magnitude smaller than Ac in all ISM phases. In a 
magnetized, turbulent, plasma, the most efficient of col- 
lisionless processes is scattering off magnetic fluctuations 
of size vb — r g , which in duce resonant p itch angle scat- 
tering of positrons (e.g. iKulsrudl 120051 and references 
therein), or non-resonant intera ctions with fluctuations 
on sc ales just above r g (see e.g. iRagotl . 120061 : iToptvginl 
1985). The transport of energetic (>GeV) GCR is driven 
by such collisionless processes (Sec. IV. B. 2), but in the 
case of MeV positrons the situation is not clear, because 
there is no observational evidence on the level of intestel- 
lar turbulence at such small scales (although one may 
reasonably expect that it can be quite high in the vicin- 
ity of some positron sources, e.g. supernovae). 

Because of the lack of observational evidence on the 
dominant regime of propagation of low energy positrons, 
in this Section we consider both regimes (Sec. IVI.AI and 



For a complete description of the ISM gas phases see Section 

MM 



Sec. [VLBj respectively). Some implications of this anal- 
ysis for the spatial morphology of the Galactic 511 keV 
emission are presented in Sec. IVI.CI 



A. Transport by collisions 

The propagation of positrons in different gas phases 
in the collisio nal regime is tho roughly studied in the re- 
cent work of Jean et al. (2009), with a Monte Carlo code 
describing e + "cooling" and annihilation 29 . The Monte 
Carlo simulation consists in calculating the trajectory of 
positrons (injected through a point source) along mag- 
netic field lines, taking into account collisions that change 
their pitch angle and energy. Their initial velocity distri- 
bution is assumed to be isotropic, so the cosine of their 
initial pitch- angle is distributed uniformly between -1 and 
1. The initial kinetic energy of positrons is in the range 
10 keV < E < 10 MeV (appropriate for radioactivity 
emitted positrons). The simulation provide the spatial 
position of positrons at the end of their life. Analysis 
of a large number of test particles leads to the spatial 
distribution of the annihilation sites. 



The process of e+ cooling encompasses two distinctly different 
phases: slowing down and thermalization. During the slowing- 
down period, the e+ kinetic energy is above the charge-exchange 
threshold and the particles cool through binary or Coulomb in- 
teractions. Once the e + kinetic energy is between the charge- 
exchange energy and the thermal energy of the gas, the parti- 
cles enter into a thermalization phase where the elastic and/or 
Coulomb interactions dominate. The thermalization ends once 
the e + kinetic energy becom es close to the g as thermal en- 
ergy and positrons annihilate. iJean et al. (2009) found that the 
slowing-down process is always the longer of the two. 
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The simulations show that the spatial field-aligned dis- 
tributions of positrons at the end of their life are nearly 
uniform out to the maximum distance traveled along the 
field lines (corresponding to initial zero pitch angle). This 
implies that the pitch angle does not change significantly 
during the slowing-down period. For a given density, the 
most significant effect on positron propagation is found 
to be due to the variation of the ionization fraction in 
the neutral gas (warm neutral, cold neutral and molecu- 
lar): as the ionization fraction increases, Coulomb scat- 
tering quickly becomes the dominant process of e + energy 
losses. 

Fig.l32lshows the spatial extents of the final (before an- 
nihilation) distribution of positrons along and perpendic- 
ular to the regular magnetic field, respectively, as a func- 
tion of their initial kinetic energy, for each ISM phase. 
The extent along the regular magnetic field corresponds 
to approximately twice the maximum distance travelled 
by positrons and gives the size of the annihilation site 
(even though positrons originate from a point source). 
A "realistic" Galactic magnetic field was assumed, com- 
posed of an average magnetic uniform field B plus a ran- 
dom field of intensity SB c± B, due to the turbulent mo- 
tions of the gas (Sec. IIII.D|) . The effective distances 
traveled by positrons (i.e. in a straight line from their 
initial to final position) following the chaotic field lines 
are smaller than in the case of a simple, uniform magnetic 
field, but only by 25%. The reason for such a surpris- 
ingly small effect is that, because of the assumed injection 
scale of turbulence from supernovae (10-100 pc) and the 
adopted turbulence power-spectrum (Fig. [TT]h the num- 
ber of small scale fluctuations interacting with positrons 
is quite small. For that reason, the distances calculated 
by Monte- Carlo for a realistic Galactic magnetic field in 
Fig. [32] are quite similar to the stopping distances D 
calculated semi-analytically for a uniform average mag- 
netic field in Fig. [3TJ In other words, the magnetic field 
rapidly orients the motion of randomly injected positrons 
along its direction and MeV positrons travel an effective 
distance ~ S ~10 n _1 kpc (where n (cm -3 ) is the ISM 
density) from their origin before stopping by Coulomb 
collisions. As already stressed, these distances are larger 
than typical sizes of the warm ISM phases where Galactic 
positrons annihilate, according to the spectral analysis of 
511 keV emission (Sec. IV. E|) . This may imply either that 
the adopted model for turbulence is inadequate (for in- 
stance, turbulence can show anisotropic or intermittent 
behaviou r, instead o f the is otropic turbulence cascade as- 
sumed in ljean et aZ.L 12009) or that positron transport (at 
least, in the warm phases) is dominated by collisionless 
processes. The latter issue is discussed in the next sec- 
tion. 



B. Wave-particle interactions and collisionless transport 

In the case of weak homogeneous turbulence - electro- 
magnetic fluctuations over several orders of magnitude 
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FIG. 32 Minimum and maximum extents of the spatial dis- 
tributions of 1 MeV positrons slowing down to 100 eV, along 
(top) and perpendicular (bottom) to the regular Galactic 
magnetic field, taking into account the turbulent behavior 
of the field lines as w ell as realist ic values for the density in 
each ISM phase (from iJean et all 120091 . 



in wavelength, homogenously distributed in space - the 
wave-particle interactions and collisionless charged parti- 
cle transport are rather well understood. Since the gyro- 
radii of particles in interstellar magnetic fields are usually 
much smaller than the relevant spatial scales, one may 
use the so-called "gyro-phase averaged distribution" of 
particles, which depends on four variables: time t, spa- 
tial coordinate s along the field lines, momentum p, and 
pitch angle \i. The evolution of the particle distribu- 
tion, f(t, s. p, /i), can be described by a Fokker-Planck 
equation (M elrosel Il980l ), while the injected spectrum of 
positrons and its spatial distribution depend on the na- 
ture of the e + source (see Sec. IIV|) . 

Positrons undergo pitch-angle scattering described by 
the angular diffusion coefficient and stochastic accel- 
eration by interaction with plasma turbulence, described 
by the momentum diffusion coefficient D pp . In collision- 
less turbulent interstellar plasmas, the kinetic coefficients 
in the Fokker-Planck equation are dominated by resonant 
wave-positron interactions involving both cyclotron and 
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Cherenkov resonances 30 ([Melrosel Il980l ) and they can be 
evaluated in the "quasi-linear" theory, valid in the weak 
turbulence regime 31 . 

Magnetohydrodynamical (MHD) waves 32 are of prime 
interest for the transport of low energy positrons . They 
can exist only at frequencies lower than the proton cy- 
clotron frequency Q cp = (qB/m p c), where q is the ele- 
mentary charge, and m p is the proton mass. They are 
damped either by collisional effects (mainly viscous fric- 
tion and ion-neutral collisions) at low frequencies or by 
the Landau damping 33 (due to thermal protons) at fre- 
quencies approaching Q cp . Higher- frequency waves are 
potentially important as well. However, whistler waves 
produced by electron-proton plasma instabilities, which 
are the most interesting waves in this frequency domain, 
are right-handed polarized; therefore, they cannot be 
in resonance with positrons, unless positron or proton 
flows generate their own waves. We shall restrict our 
analysis to MHD waves with frequency u <C ^ se , where 
Q se = Vtce/l and ^ ce = qB/m e c are the synchrotron and 
the cyclotron frequencies of the positron, respectively. 

Collisionless processes may result in efficient decel- 
eration - or "cooling" - of fast positrons (see, e.g. 
iPetrosian and Bvkovl 120081 and references therein). The 
collisionless positron scattering due to particle-wave in- 
teraction is efficient in warm ionized phases of the 
ISM where the wa ve damping is not too strong (see 
eg. iKulsrudL l2005h . Moreover, adiabatic deceleration 
of positrons in jets or expanding shells (for example in 
SNRs) results in positron cooling, even without Coulomb 
collisions; this occurs if the positron mean free path, 
which is dominated by e + scattering by waves, is shorter 
than the typical scale of bulk plasma motion 34 . Hence, 
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FIG. 33 Mean free path of energetic particles parallel to mag- 
netic field direction, as a function of particle rigidity in the 
interplanetary medium; circles and upward pointing trian- 
gles d enote the measured values and upper limits, respectively 
(from iBieber et all 1 1994 ) . 



wave-particle interactions, both resonant and adiabatic 
non-resonant, could result in particle deceleration or re- 
acceleration, depending strongly on the local conditions. 
The poor knowledge of those local conditions (concern- 
ing essentially the small-scale magnetic turbulence in 
the ISM) imposes a case by case study and precludes 
any generic conclusions to be drawn. Studies of wave- 
positron interactions have then to rely either on in situ 
measurements in the solar wind (Sec. IVI.B.1[) or on theo- 
retical modeling of the properties of the MHD turbulence 
(Sec. lVI.B.2!l . 



The electron cyclotron resonance involves the gyro-motion of 
electrons (or positrons) perpendicularly to the magnetic field: 
the transverse electric field associated with the wave rotates at 
the same velocity and in the same direction with the particles, 
which absorb its energy and accelerate. The Cherenkov reso- 
nance involves particle motion along the magnetic field lines. 
To calculate the diffusion coefficients beyond the quasi-linear ap- 
proach (e.g. for particle transport and acceleration by strong 
turbulent fluctuations ) the re-normal ization equati ons were de - 
veloped bv iBvkov and Toptvginl (119931 ) and IZank eTal\ (I2004T ). 
assuming that the particle propagation regime is diffusive 
Under MHD waves we mean shear Alfven waves and fast mag- 
netosonic waves. Because the phase velocity of the magne- 
tosonic mode is almost always larger than the Alfven velocity 
Va, the magnetosonic wave is often called the "fast" hydromag- 
netic wave. The dynamics of the third MHD mode, the "slow" 
wave, has been shown to be entirely controlle d by the Alfven 
wave cascade by Lith wick and Goldreichl (120011 ); the slow wave 
spectrum is basically the same as the Alfven wave spectrum. 
Landau damping occurs due to the energy exchange between a 
wave with phase velocity vp and particles with velocity v ~ vp , 
which can interact strongly with the wave. Particles with v < vp 
will be accelerated by the wave electric field, while those with 
v > vp will be decelerated, losing energy to the wave. 
The adiabatic deceleration of a positron diffusing in an expand- 
ing shell has a typical time scale r a d ~ |V-U(r,£)| , where 



1. Local e + transport in the ISM and re-acceleration 

In the solar wind - the only natural laboratory for 
direct study of the propagation effects - sub-MeV elec- 
trons and positrons resonantly interact with waves of 
high enough frequencies that contain only a tiny frac- 
tion o f energ y of the turbulent wave cascade (see e.g. 
iRagotL l2005f ). At 1 AU from the Sun 35 the solar wind 
turbulence steepens at a cut-off scale Vx/Q, ce ~10 7 cm, 
where the Alfven velocity Va ~ 2.1 10 5 B^q n~ x l 2 cm 
s _1 and n(cm -3 ) is the local plasma density. This means 
that for pitch-angles /i < Va/v (where v is the positron 
velocity) sub-MeV particles can no longer gyroresonate 
with waves in the inertial range of the solar wind tur- 
bulence. However, if non-resonant magnetosonic (i.e. 
" fast" ) waves of lower frequencies are present in the solar 
wind, they could dominate the propagation of sub-MeV 



U(r, t) is the bulk plasma velocity. The cooling time in expand- 
ing wind, shell or jet (with V • U(r, t) > 0) can be much shorter 
than the Coulomb stopping time in a rarefied plasma. 
1 Astronomical unit (AU)=1.5 10 8 km, the Earth-Sun distance. 
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particles (e.g. iRaeotl 120061 : iToptygfrl Il985h . Depending 
on the level of non-resonant fast-mode wave turbulence 
and the injection distribution of electrons, this can in- 
duce diffusion in energy space, i.e. Fermi re-acceleration 
of su b-MeV electrons in the inner heliosphere (jRagotl . 
l2006h . The conditions for efficient re-acceleration can be 
estimated as follows: 

If positron velocities are above the thermal electron 
velocities, the ratio of the energy diffusion rate D pp /p 2 
and the positron pitch angle scattering rate (both 
in s _1 ) is 

D PP /P 2 D~* fa (Va/v) 2 <C 1, (21) 

(e.g. iPetrosian and Bvkovl 120081 and references therein) . 
The re-acceleration effect is non-negligible when the re- 
acceleration timescale is shorter than the energy loss 
timescale (or the e + propagation timescale). Stated dif- 
ferently, the acceleration rate D pp /p 2 must exceed the 
momentum loss rate Pl/p = [dE/dt]/ \m e v 2 ^), where - 
apart from fastly expanding regions - the energy loss rate 
of sub-MeV particles is dominated by Coulomb losses 
(Eq. ([14]) . Therefore, for re-acceleration to be impor- 
tant, the positron scattering rate by MHD waves 
must satisfy the condition 

n 2 

23 w >0.02 i 8- 1 7" 1 oT" (O, ( 22 ) 

which corresponds to a particle mean free path A(= 
1/3 vD~^) satisfying 

B 2 

A<0.1/3 2 7 ^- (AU). (23) 

nj 

In the solar wind, the electron mean free path A reaches 
a plateau at e nergies ~1 MeV with A about 0.01 — 1 AU 
([RagotL [l999), depending on the parameters chosen for 
the interplanetary medium (IPM). This is comparable to 
typical values of A estimated from observations of elec- 
trons propagating in the inner heliosphere (see Fig. 133]). 
Interestingly, those observations show no dependence on 
electrons' rigidity (or momentum). This can be inter- 
preted as an indication that re-acceleration of low en- 
ergy particles is importan t in the case of the solar wind 
(for f urther discussion see lDrog e, 2000 and Sha lchi et all 
120061) . 

In the ISM, Eq. [23] implies that the effect of re- 
acceleration is expected to be important in the tenuous 
phases filled with strong MHD turbulence. Fig. [34] shows 
the ISM parameter space (number densities vs. mean 
magnetic fields) where MeV positrons may be affected 
by re-acceleration, depending on their mean free path A, 
which depends in its turn on the level of MHD turbu- 
lence. Re-acceleration becomes progressively important 
in regions of low density and/or high magnetic field (such 
as those of the inner Bulge), resulting in longer thermal- 
ization timescales for positrons and increasing the sizes 
of corresponding e + annihilation regions. 
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FIG. 34 Regions of the ISM in the B vs n (magnetic field vs 
number density) plane. Physical conditions for HIM, WIM, 
WNM and Inner Bulge are discussed in Sec. Ill, while IPM 
stands for the interplanetary medium. The two lines corre- 
spond to Eq. [23] for MeV positrons, with the particle mean 
free path A=0.3 AU (solid, appropriate for the interplanetary 
medium from Fig. I33[) . and A=1000 AU (dotted, provided for 
illustration purposes only, since its value is unknown in the 
ISM). In regions of the ISM found to the left of those lines, 
re-acceleration effects may be important for positrons (pro- 
vided there is a sufficient level of small-scale turbulence), as 
suggested from observations of the IPM. 

Unfortunately, our knowledge of the small-scale ISM 
turbulence is limited at present. While the quasi-linear 
theory provides a simple analytical description of charged 
particle transport, recent test particle si mulations re- 
veale d some problems of this approach (e.g. iTautz et~aH 
120061 ). In particular, while non-linear effects were shown 
to be essential for partic l e tran s port in MHD turbulence 
(e.g. lie Roux and WebbL 120071 : iMatthaeus et all 120031 : 
IZank et 01.1 .12004). test-particle simulations of the coef- 
ficient D pp do not exist at present. Therefore, one has 
to rely on simulations of and use that quantity in 
order to estimate D pp from Eq. ([2T]) : in their turn, sim- 
ulations of require assumptions on the polarization 
and spectral properties of MHD waves. 

The aforementioned uncertainties affect any attempt 
to evaluate the propagation of MeV positrons in the 
Galaxy without detailed information on the magnetic 
fluctu ation spectra i n the p ropagation region. In a recent 
work, iHigdon et all (|2QQ9h adopted a particular model 
calculation of the charged particle mean free path in 
the interplanetary medium (in the quasi-linear regime), 
which they extrapolated to the conditions of the ISM (see 
Sec. VI. C for details of their model). However, in view 
of the different conditions - density, magnetic field and 
unknown level of small-scale turbulence - in the Bulge, 
such a calculation appears rather arbitrary (although not 
necessarily wrong) . Their model can be considered as a 
quantitative illustration of a possible scenario, the plau- 
sibility of which remains to be shown. 

Notice that small-scale MHD waves can be injected 
into the ISM in-situ by a variety of kinetic instabili- 
ties. A kinetic instability may be locally triggered, for 
instance, by the streaming (through the ISM) of cosmic 
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rays with bulk v elocity larger th an a few times the lo- 
cal Alfven speed (We ntzej Il974f ). Such streaming insta- 
bility is expected to develop in the intercloud medium 
and may provide particle re-acceleration, compensating 
for the strong ionization losses inside molecular clouds. 
In other terms, low-energy cosmic rays scatter off their 
self-generated wa ves and are, therefo r e, pu shed outside 
molecular clouds (Skil ling and Strong! . Il976h . The waves 
thus generated can help to confine the positrons in the 
warm phases around the molecular clouds. This mecha- 
nis m was adapted to the transport of cosmic-ray electrons 
by iMorfiiil (|l982h . but it operates over a limited range 
of positron momenta and its ability to confine positrons 
within the ionize d phases of th e ISM is questionable 
(jjean et all . l2QQ9h . iHigdon et all (j2QQ9h proposed an al- 
ternative mechanism in which positrons scatter off their 
own self-generated waves. It is possible that kinetic in- 
stabilities play a key role in the observed e + annihilation 
(mostly) in the warm phases of the ISM. 



2. Global positron transport 

Description of positron transport on a small scale 
r g < L < A (where A ~ v is the positron mean free 
path) requires the solution of the Fokker-Planck equa- 
tion. On much larger spatial scales (L> A)) the pitch- 
angle distribution is nearly isotropic. This justifies the 
use of a simplified Fokker-Planck equation, in the so- 
called diffusion approximation; the correspondin g equa- 
tion i s called the diffusion - advec t ion equation (JBykoyJ, 
l2QQlfc iBvkov and Toptvgid . Il993h - If the amplitude of 
velocity fluctuations is u > v\/L, the particle transport 
is dominated by the turbulent advection. 

It is important to note that only compressible large 
scale turbulent motions pro vide efficient re- acceleration 
([Bvkov and Toptvgid . Il993h . In the regime dominated 
by turbulent advection (TA), the re-acceleration rate 
can be approximated as £> nn _T A ~ (u/9L)p 2 , e.g. 
IBvkov and Toptvgid (|l993l 1200 lh . The re- acceleration 
effect is important if D pp -ta/p 2 > Pl/p, where the stop- 
ping rate due to Coulomb collisions Pl/p is discussed 
in Sec. IVI.B.ll Therefore, particle re-acceleration would 
significantly affect the global positron propagation in the 
bulge if the amplitude of compressible turbulence u$o 
(measured in units of 50 km/s) and the energy contain- 
ing scale L50 (measured in units of 50 pc) satisfy the 
condition 36 



U5o/L 50 > 300 



(24) 



A word of caution is in order here. The re- acceleration 
effect, in principle, could boost a substantial amount 



36 In case of important energy losses, other than from Coulomb 
collisions, higher values of u^q/L^q than given in Eq. [24] are re- 
quired. 



of mildly relativistic MeV-positrons to ultra-relativistic 
energies. It might therefore lead to violation of the 
limit on 7-ray continuum emission from in-flight annihila- 
tion of positrons imposed by Beacom and Yiiksel (2006), 
which was discussed in Section IV. In the GC vicinity 
and in the presence of a strong magnetic field, the se- 
vere synchrotron-Compton losses of relativistic positrons 
can make it impossible to accelerate them to energies 
above a few MeV by the Fermi-type mechanism. In the 
absence of synchrotron-Compton losses the differential 
spectral index of the Fermi re-accelerated positron mo- 
mentum distribution N(p) oc p~ a can be approximated 
as a = -1/2 + V / 9/4 + 36(L/i?)/2, where R is the min- 
imal size of the confinement region (if anisotropic); in 
the bulge R is likely below 1 kpc and it is much less in 
the disk. Exact estimates of the positron spectra can be 
done only if one knows the energy (or momentum) depen- 
dence of the positron mean free path in the MeV regime, 
which is governed by yet uncertain small scale magnetic 
fluctuations in the regions of interest. 

The question of whether MHD cascades may extend 
over several orders of magnitude on spatial scales is still 
a matter of active debate. Electron density power spectra 
in the local I SM have been measure d down to scales of 
10 10 cm (e.g. lArmstrong et a/.L[l995[ ). It was established 
that radio scattering in the interstellar plasma implies 
widespread inhomogeneities in density with a power law 
spatial spectrum and suggests a turbulent origin. Obser- 
vations suggest that the small scale plas ma turbulence 
is often highly aniso t ropic and clumpy (iBrisken et all , 
l2010l : iDesai and Fevl l2QQlfc IShebalin et all Il983h ~s 
fluctuations are likely to concentrate in filamentary 
densi ty struc t ures aligned by the local magnetic field 
(e.g. iHigdonl Il984l ). Since the scattering would be- 
come isotropic if the irregularities were uniformly dis- 
tributed over many le ngth scales of the magnetic field, 
iRickett and Colesl (|2004l ) concluded that the plasma tur- 
bulence is distributed intermittently. 

To estimate the positron scattering effect from the elec- 
tron density measurements one has to assume a relation- 
ship between the amplitudes of density and magnetic field 
fluctuations. Fluctuations due to magnetosonic waves 
could contribute to non-resonant scattering of MeV par- 
ticles or - if they survive down to small scales - to reso- 
nant scattering. However, if the observed electron density 
fluctuations are simply due to entropy- type ("isobaric") 
fluctuations, they are inefficient in e + scattering. 

Depending on the ISM phase, the MHD wa ves can suf- 
fer fr om collisionless or collisional damping. Ijean et all 
(2009) concluded that in the neutral atomic and molecu- 
lar phases of the ISM the Alfven and fast magnetosonic 
wave cascades are both cut off by ion- neutral collisions, 
on scales considerably larger than the gyroradius of MeV 
positrons; therefore, MHD waves cannot resonantly in- 
teract with positrons. The situation is different in the 
ionized phases of the ISM, where the Alfven wave cas- 
cade suffers insignificant (collisional) damping down to 
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TABLE XII Collionless and collisional transport processes in 
ISM phases (see Table IV). 









ISM phase 




Transport processes 


CM 


CN 


WN 


WI 


HI 


NR MHD modes: Large scale 
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R MHD modes Alfven 
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R MHD modes Fast MS 
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N 


N(Y) 


N(Y) 


Streaming modes 


Y? 


Y? 


Y? 


Y 


Y 


Collisions 


Y 


Y 


Y 


Y 
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Y (N) means that the process can (cannot) take place. NR 
(R) stands for nonresonant (resonant) MHD motions at 
large scales. Streaming modes are generated through the 
streaming instability produced by the CRs propagating 
away from their sources The streaming processes not yet 
studied in depth are indicated with question marks. 



the thermal proton Coulomb mean free-path 37 . This 
leaves some room for possible resonant interactions of 
the Alfven wave cascade with positrons 38 but the effect 
has not been properly investigated up to now. 

In summary, contrary to the case of the solar wind 
(where nonresonant compressible perturbations may con- 
trol the positron mean free path), the cascade of MHD 

10 - 100 



waves in the ISM (injected at scales Li 



mj 



pc by supernovae) is expected to be damped at scales 
well above r g . Small scale MHD (and whistler) waves 
could be generated by anisotropic distributions of ener- 
getic particles, through streaming instabilities, but the 
waves could b e damp ed by ion-neutral collisions as well 
( Higdon et all l2QQ9h : in any case, the effects have not 
been studied in detail up to now. Streaming instabilities 
could help to confine positrons at the border of the molec- 
ular clouds, thus enhancing the fraction of positrons that 
annihilate in warm phases. In ISM regions of low density 
(and/or high magnetic field) with even a moderate level 
of MHD turbulence (see Fig. [32]) particle re-acceleration 
could substantially suppress the positron stopping by 
Coulomb collisions (even in the absence of any other 
cooling effect, e.g. adiabatic losses), thus providing a 
possibility to extend largely the volume filled with anni- 
hilating positrons. Most of the current models rely on 
quasi-linear wave treatment, but nonlinear effects could 



37 On smaller scales, the cascade enters the collisionless regime, and 
cuts off by linear Landau damping around the proton inertial 
length c/cjpp (where (jU pp = 47rnne 2 /m p ). 

38 Notice that, theoretical considerations show that the energy 
transfer in the Alfven cascade proceeds mostly through the 
interaction of oppositely propagating wave packets. The dis- 
tortion of the wave packets during the interaction produces 
anisotropic fluctuations e longated along the mean magnetic field 
([Lithwick and Go ldreich, 200l|). Using numerical simulation, the 
magnetosonic cascade has been fo und to keep its isotropy a long 
the whole range of spatial scales (|Yan and La zarian. 2004|). In 
case of strong anisotropy of the cascade, the scattering efficiency 
of MeV positrons by small scale Alfven waves could be consid- 
erably reduced while the effect of magnetosonic waves could be 
enhanced. 



substantially change some estimates. Notice that, in gen- 
eral, a quasi-homogeneous distribution of the sources of 
turbulence is assumed, although the intermittency effects 
(due to e.g. large scale ISM shocks) can modify the anal- 



ysis. In Table IXHI we list processes that are potentially 
important for po sitron propagation in the ISM (see also 
I Jean et a/ll2009l ) and which require further investigation. 
Our current understanding of MeV positron propagation 
does not allow one to conclude whether such positrons 
undergo strong diffusion or essentially free propagation. 



C. Implications of e + propagation for 511 keV emission 

The implications of low energy positron tra nsport for 
the G alactic 511 keV emission were raised bv IPrantzosI 
(|2QQ6h . who pointed out that the morphology of the 511 
keV emission does not necessarily reflect the morphol- 
ogy of the underlying e + source distribution. He no- 
ticed that e + from SNIa are released in the hot and rar- 
efied ionized medium, since the scaleheight of SNIa is 
considerably larger than the scaleheight of the disk of 
cool, dense, gas (see Fig. [10]). The e + propagation dis- 
tances are then quite large (Fig. l3T]h allowing e + from 
the disk to annihilate far away from their sources (per- 
haps in the halo, where a low surface brightness emis- 
sion should be expected); this fact may considerably re- 
duce the bulge/disk ratio of 511 keV emission of any 
class of astrophysical e + sources, thus alleviating the 
morphology problem discussed in Sec. IV. D. Models 
of cosmic ray propagation in the Galaxy (like those de- 
scribed in Sec. IIV.B.2P usually adopt an iso tropic diffu- 
sion coefficient. iGebauer and de Boerl ( 20091 ) introduced 
recently an anisotropic diffusion coefficient in the GAL- 
PROP code to simulate the advection of cosmic rays by 
a Galactic wind. They find that, by adopting observa- 
tionally derived wind velocities, their scheme naturally 
produces large escape fractions (>50%) of positrons from 
the disk. 

Furthermore, IPrantzosI ( 2006[ ) suggested that if the 
halo magnetic field of the Milky Way has a strong 
polo i dal component, as suggested by s everal authors (e.g. 
iHanl 120041 : iProuza and gmidal 120031 and Fig. [35j top) 
then some positrons escaping the disk may be channeled 
into the bulge and annihilate there, enhancing even more 
the bulge/disk e + annihilation ratio; he noticed that, in 
that case, positrons from SNIa may suffice to explain 
quantitatively both the total observed e + annihilation rate 



(~2 10 43 e + s _1 ) and the corresponding bulge/disk ratio, 
provided that the escaping e + fraction from SNIa is ~4% . 
As discussed in Sec. IIII.Dl it is rather unlikely - although 
it cannot yet be ruled out - that the poloidal compo- 
nent of the regular Galactic magnetic field is close to a 
dipole. Observations of external spirals suggest rather an 
X-shaped halo field (e.g. Fig.[35j bottom), in which case it 
would be difficult for disk positrons to find their way into 
the bulge. Still, the issue is of considerable interest and 
urgently calls for a better assessment of the poorly known 
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FIG. 35 Possible configurations of the large scale magnetic 
field of the Milky Way. Top: As derived from Faraday polar- 
ization measurements of the MW, according to lHar 
Bottom: Scetch of the observable components of the large 
scale magnetic field of the disk galaxy NGC253 (which shows, 
however, signs of starburst activity, unlike the Milky Way); 
the halo magnetic field is even and pointin g; outwards, wherea s 
the dotted parts are not observed (from iHeesen et al.l\200$ ). 



gl obal configuration o f the Galactic magnetic field. 

iHigdon et all (2009) suggested that positron propaga- 
tion may be the key for understanding not only the spa- 
tial morphology of the 511 keV emission, but also its 
spectral properties. They made the bold assumption that 
radioactivity (from 26 Al, 44 Ti and, mostly, from 56 Co, 
see Sec. IV. A) is the sole e + source in the Galaxy. They 
considered (i) a fairly detailed description of the various 
phases of the ISM (Sec. III.B) and (ii) a particular phe- 
nomenological model of collisionless scattering of MeV 
positrons by turbulent fluctuations of the ISM, that was 
used to describe energetic particle (electron or proton) 
propagation in the interplanetary medium. Depending 
on the nature of the medium, positrons are assumed to 
propagate either by diffusion along magnetic flux tubes 
(in ionized media, where turbulence cascades down to 
the gyroradius r g ) or by streaming with mean velocity 
v = j3c (in neutral media, where turbulence is quenched 
by ion- neutral collisions at scales >> r g ). 



Putting together the aforementioned ingredients, 
IHigdon et all (|2QQ9h proceed then in an impressive cal- 
culation of positron production rates from radioactiv- 
ity along the Galaxy, typical distances of e + propaga- 
tion in the corresponding ISM phases and probabilities 
that positrons will finally annihilate in one or another 
of those phases. In the end, they find excellent agree- 
ment with each and every observable of the 511 keV 
emission available so far (Sec. II. D): spatial morphol- 
ogy, i.e. bulge/disk annihilation ratio, spectral futures 
(including a narrow and a broad 511 keV line) and even 
the claimed asymmetry between fluxes from negative and 
positive longitudes; they explain the latter as due to the 
corresponding overall asymmetry of the spiral arm pat- 
tern of the Milky Way disk, as viewed from the Sun. 

The work of IHigdon eTaR (2009) constitutes the first 
study of Galactic e + production, propagation and anni- 
hilation in a "global" framework, trying to include all (or 
most of) the various aspects of this complex topic and to 
account for all the available observational data. How- 
ever, its extremely precise "predictions" for the various 
properties of the 511 keV emission (which fit extremely 
well - to better than 10% - each and every observable), 
concealed the various uncertainties of the problem. For 
instance, it is assumed that most Galactic positrons re- 
sult from 56 Co produced in SNIa: but, as discussed in 
Sec. IV. A. 4, recent observations suggest that the e + es- 
cape fraction from such objects is very small, at least at 
late times. Only an early e + escape (due, perhaps, to 
3D effects, not yet theoretically treated and unobserv- 
able at present) could make SNIa plausible e + sources 
again. And even if that were the case, the poorly known 
SFR of the bulge does not allow one to estimate the SNIa 
rate of that region to better than a factor of two (see Ta- 
ble [Vj). The same criticism applies to their treatment of 
e + propagation. Their prescription for deriving the e + 
mean free path is, strictly speaking, valid in the frame- 
work of interplanetary plasma, while the properties of 
the ISM and magnetic fields in the inner Galactic re- 
gions are too poorly understood to allow for any strong 
conclusions on the nature of turbulence there. For in- 
stance, reacceleration - which extends the e + annihila- 
tion zone - is expected to be particularly efficient in the 
ionized, low density, medium of the bulge (Sec. VI.B.l). 
As for the claim that the observed spatial asymmetry of 
the 511 keV emission is simply due to the corresponding 
asymmetry of the spiral arms as viewed by the Sun, it 
is apparently in contradiction with the fact that no such 
strong asymmetry is observed in the 1.8 MeV emission 
of 26 Al (a product of massive stars), which should be a 
direct tracer of spiral arms, since 26 Al nuclei travel much 
slower than positrons. 

One of the most interesting results of the work of 
IHigdon et all ( 20091 ) is, perhaps, their prediction of a spa- 
tial differentiation between the broad and narrow 511 keV 
line components, revealed by the spectral analysis of SPI 
data (Sec. II. C). The former should occur in the "middle 
bulge" (the region of radius 0.5 <r<1.5 kpc, according to 
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the authors' terminology), through in-flight positronium 
formation, while the latter in the "inner bulge" (r<0.5 
kpc) through positronium formation at thermal energies. 
An analogous spatial differentiation is predicted to ex- 
ist between the bulge and disk 511 keV line components. 
It is unlikely that such predictions will be checked with 
INTEGRAL/SPI observations, but they certainly hold 

important clues to positro n propagation issue s. 

In a companion p aper, H?i ngenfelt er et al\ ( 20091 ) ar- 
gued that, since the Iffigdon et all (|2QQ9) model (with 
propagation of radioactively produced positrons from 
SNIa) can fully account for the observations, dark matter 
should be excluded as a major positron source. However, 
such a conclusion is premature, since it is not yet clear 
whether positrons from SNIa escape at all. It is true that, 
since the confirmation of the disk 511 keV emission, DM 
has lost a lot of its "appeal" as positron source, but it 
cannot yet be excluded as such, at least for the bulge: in- 
deed, in that case e + propagation merely smears out the 
spatial profile of the 511 keV emission, and even decaying 
DM cannot be excluded then (contrary to the arguments 
of Sec. IV. D, which neglect positron propagation). 



VII. SUMMARY AND PERSPECTIVES 

The Galactic 511 keV emission from e + annihilation 
is the first 7-ray line detected from outside the solar 
system. Its unambiguous identification came soon after 
its detection, with high resolution Ge-detectors aboard 
balloon experiments. However, its spatial morphology 
remained elusive for almost three decades after its first 
detection. Only long running experiments aboard satel- 
lites could tackle this issue, in view of the importance 
of the treatment of the background at those energies 
(mostly created inside the detectors by cosmic ray inter- 
actions). Observations in the 1990s by OSSE/CGRO of- 
fered the first hints for an abnormally high bulge / disk ra- 
tio (compared to the situation at any other wavelength). 
That property was firmly established only after obser- 
vations in the 2000s with SPI/INTEGRAL, which fur- 
ther detected for the fist time an unambiguous disk emis- 
sion. The latter appears to be asymmetric, according to 
Weidens pointner et al. I (l2008ah . with emission from neg- 
ative longitudes being 80% brighter than from positive 



ones; ho wever, that claim is not supp orted by a different 
analysis ( iBouchet et 
confirmed. 



um is no t supp orted by ; 
3, 120081 [2010) and has 



yet to be 



According to the latest ima ging analysis of SPI data 
([Weidenspointner et all [2008a) the total Galactic e + an- 
nihilation rate is at least N e + ~2 10 43 s _1 , with a lumi- 
nosity bulge/disk ratio B/D=1.4. This model (Table I), 
is further refined by considering a narrow (FWHM = 
3°) and a broad (FWHM = 11°) bulge, the former con- 
tributing to ~35% of the total bulge emission. However, 
the data analysis also allows for other morphologies, in- 
volving extended regions of low surface brightness but 
high total emissivity, e.g. a "halo" of total N e + ~ 3 



10 43 s -1 and a thin disk of N e + -5 10 42 s~\ leading to 
a high B/D~6 (Table I). Obviously, the poorly known 
configuration of the Galactic 511 keV emission precludes 
at present the formulation of definitive statements about 
the origin of annihilating positrons. 

Information on the origin of those positrons is also 
obtained via the spectral analysis of the 511 keV emis- 
sion: the observed flux at ~MeV energies from the in- 
ner Galaxy constrains the initial energy of the positrons 
to less than a few MeV (otherwise the emission from 
in-flight annihilation would exceed the observed flux). 
Moreover, the spectral analysis provides important in- 
formation on the physical properties of the e + annihila- 
tion sites. The large positronium fraction fp s ~94-97 
% implies that positrons annihilate mostly at low ener- 
gies, since direct annihilation cross-sections are impor- 
tant only at high energies (Sec. V.D). The overall spec- 
tral shape suggests that annihilation occurs mostly in 
warm (T~8 000 K) media, at about equal amounts in 
neutral and ionized phases but it cannot be excluded that 
less than 23% of annihilation occurs in the Cold neutral 
medium (T ~ 80 K; Sec. V. E); annihilation in the neu- 
tral media may account for the presence of a broad 511 
keV line component (FWHM ~5 keV) and the annihi- 
lation in the warm ionized medium for the narrow one 
(FWHM -1 keV). 

Among the various astrophysical sources of positrons 
proposed so far, the only one known with certainty to 
release e + in the ISM is /3 + radioactivity of 26 Al; the ob- 
served intensity of its characteristic 1.8 MeV emission in 
the Galaxy corresponds to ~3-4 10 42 e + s _1 (Sec. II. C. 2). 
A similar amount is expected from the decay of 44 Ti, on 
the grounds of nucleosynthesis arguments (Sec. IV. A. 2). 
Both radionuclides are produced mostly in massive stars 
and their positrons should be released along the Galac- 
tic plane, as traced by the 1.8 MeV emission; they could 
thus account for the observed disk 511 keV emission. 

Radioactivity of 56 Co from SNIa was traditionally con- 
sidered to be the major e + producer in the Galaxy. Both 
the typical 56 Ni yield of a SNIa and the Galactic SNIa 
rate are rather well constrained, resulting in 5 10 44 e + s _1 
produced inside SNIa. If only i esc ~4% of them escape 
the supernova to annihilate in the ISM, the observed total 
e + annihilation rate can be readily explained. However, 
observations of two SNIa, interpreted in the framework 
of 1-D (stratified) models, suggest that the positron es- 
cape fraction is negligible at late times. On the other 
hand, both observations of early spectra and 3-D models 
of SNIa suggest that a sizeable fraction of 56 Ni is found 
at high velocity (close to the surface), making - perhaps 
- easier the escape of 56 Co positrons (Sec. IV. A. 4). In 
our opinion, SNIa remain a serious candidate, with a po- 
tential Galactic yield of 2 10 43 e + s _1 . But the expected 
spatial distribution of SNIa in the Galaxy corresponds to 
a much smaller B/D ratio than that of the observed 511 
keV profile. 

Most of the other astrophysical candidates can be con- 
strained to be only minor e + sources, on the basis of 
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either weak e + yields (novae, Galactic cosmic rays), high 
e + energy (compact objects, like pulsars or magnetars), 
spatial morphology of sources (hypernovae, gamma ray 
bursts) or a combination of those features (e.g. cosmic 
rays), as discussed in Sec. IV. E. Only two astrophysical 
candidates remain as potentially important contributors: 
LMXRBs (or the microquasar variant of that class of 
sources, Sec. IV. B. 3) and the supermassive black hole at 
the Galactic center (Sec. IV. B. 4). It should be stressed 
that there is no evidence that either of those sources pro- 
duces positrons and the e + yields evaluated by various 
authors are close to upper limits rather than typical val- 
ues (Sec. IV.D.l). Furthermore, because of the current 
low activity of the central MBH (much lower than that of 
LMXRBs) it has to be assumed that the source was much 
more active in the past, thus dropping the assumption of 
"steady state" between e + production and annihilation, 
which is likely in all other cases. 

The observed spatial distribution of LMXRBs is sim- 
ilar to the theoretically derived one for SNIa (Sec. 
IV. D. 3), as expected, since both classes of sources have 
old and young stellar components; however, none of them 
has the large B/D ratio of the observed 511 keV emission. 
The only hint that LMXRBs may contribute to, at least, 
the disk 511 keV emission stems from the asymmetric dis- 
tribution of the hard LMXRBs in the 3rd IBIS catalogue, 
of similar magnitude to the detected 511 keV emission; 
however, the former asymmetry is not confirmed (the 
4th IBIS catalogue shows no strong evidence for such an 
asymmetry) and even if the 511 keV asymmetry is con- 
firmed by future analysis, there is considerable debate on 
whether such a similarity has a causal origin or is just 
fortuitous (Sec. IV.D.3). 

Dark matter (DM) has been proposed as an alterna- 
tive e + source, at least for the bulge 511 keV emission; in 
principle, it could complement disk emission originating 
from radioactivity of 26 Al and 44 Ti or 56 Co. Observations 
of the MeV continuum from the inner Galaxy constrain 
the large phase space of DM properties. The mass of 
annihilating or decaying DM particles should be smaller 
than a few MeV, otherwise their in-flight annihilation 
would overproduce the MeV continuum (Sec. II.C.l and 
IV.B.l). Scalar light DM particles with fermionic inter- 
actions still appear as a possible candidate; alternatively, 
the collisional de-excitation of heavy (100 GeV) DM par- 
ticles could provide the required positrons, provided the 
energy separation between their excited levels is in the 
MeV range (Sec. IV. C). On the other hand, the ob- 
served spatial profile of the 511 keV emission constrains 
the production mode of DM positrons, if it is assumed 
that they annihilate close to their production region: only 
" cuspy" profiles are allowed in the case of annihilating or 
de-exciting DM particles (for which p 1 cx p 2 DM ), while 
decaying DM particles (for which p 1 cx Pdm) are ex- 
cluded (Sec. IV.D.3); the problem is that observations 
of external galaxies suggest rather flat, not cuspy, DM 
profiles (Sec III.E). 

Positrons produced in the hot, tenuous plasma filling 



the bulge (either from SNIa, LMXRBs or DM), have to 
travel long distances before slowing down and annihilat- 
ing. This is corroborated by the spectral analysis, which 
suggests that positrons annihilate in warm gas: such gas 
is filling mostly the inner bulge. Positron propagation 
appears then unavoidable, undermining the assumption 
that the e + production and annihilation profiles are cor- 
related, at least in the bulge. A similar situation should 
hold for positrons produced away from the plane of the 
disk (i.e. from SNIa or LMXRBs), which is also domi- 
nated by hot, tenuous gas. The situation is less clear for 
positrons produced by massive star radioactivity, in the 
plane of the disk and inside spiral arms: although some 
of them may fill hot bubbles and cavities created by the 
SN explosions and ultimately escape from the disk, an- 
other fraction may annihilate in closeby dense molecular 
clouds. Propagation of MeV positrons in the ISM may 
then hold the key to understanding the 511 keV emis- 
sion. It depends on the physical properties of the ISM 
(density, ionization, Sec. V) but also on the properties 
of turbulence and magnetic field configuration (Sec. VI). 
Preliminary attempts to evaluate the extent of positron 
propagation (Sec. VI. A) and their implications for the 
Galactic 511 keV emission (Sec. VI. D) are promising 
in that respect, but the situation is far from clear at 
present: the entanglement between the various uncer- 
tainties (concerning e + sources, e + propagation and an- 
nihilation sites) does not allow any strong conclusions to 
be drawn. 

More than 30 years after its discovery, the origin of the 
first extra-solar 7-ray line remains unknown. Progress in 
the field will require advances in several directions: 

(i) Observations of 511 keV emission: what is the 
true spatial distribution of the emission? how far do 
the spheroid and disk extend ? are there yet undetected 
regions of low surface brightness? is the disk emission 
asymmetric indeed? how do the 1.8 MeV and 511 keV 
disk emissions compare to each other? A much deeper 
exposure of the Galaxy and a better understanding of the 
backgrounds will be required to tackle those issues. Even 
if INTEGRAL'S mission is extended to 2012, it may not 
provide the answers; unfortunately, no other mission of 
similar scope is on the horizon. 

(ii) Physics of e + sources: what is the e + escaping 
fraction in SNIa ? what is the SNIa rate in the inner (star 
forming) and in the outer (inactive) bulge? what are the 
e + yields, activity timescales, and spatial distribution in 
the inner bulge of LMXRBs or microquasars? how can 
the past level of activity of the central supermassive black 
hole be reliably inferred? 

( Hi ) Positron propagation: what is the large scale con- 
figuration of the Galactic magnetic field? what are the 
properties of interstellar plasma turbulence and how do 
they affect the positron transport? what are the domi- 
nant propagation modes of positrons and what is the role 
of re-acceleration? 

The many facets of the Galactic 511 keV emission make 
this problem one of the most intriguing problems in high 
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energy astrophysics today and for many years to come. 
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